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PREFACE  TO  REVISION  A OF  THE  "PAC"  USER’S  MANUAL 


The  original  publication  of  the  PAC  User’s  Manual  issued  in 
August  of  198  0 is  now  somewhat  dated.  In  the  time  that  has  elapsed 
since  the  original  manual  was  written,  several  modifications  have 
been  made.  The  reason  for  the  modifications  have  been  primarily  two- 
fold. First  some  program  "bugs"  were  isolated  in  the  program  and 
reguired  fixing.  Second,  certain  enhancements  have  been  made  to 
the  program  that  should  be  documented. 

In  the  list  that  follows  these  modifications  will  be  described. 

a)  The  parameters  "DELTTB"  and  "LT"  have  been  removed  from  the 
program  input  file.  These  parameters  are  now  solved  for  in 
the  program  in  terms  of  other  input  parameters. 

b)  The  parameters  "XSTOP",  "STEP",  "PINTl"  and  "PINT2"  which  were 
previosly  assigned  values  internally  in  the  program  have  now 
been  placed  in  the  input  file  for  convenience.  These  values 
are  the  time  after  T-0  at  which  the  user  wishes  computation 

to  cease,  the  integration  interval  step  size,  and  the  print 
interval  for  the  "PRINTl"  and  PRINT2"  subroutine  block  out- 
put respectively. 

c)  The  computation  of  HIC  (Head  Injury  Crition)  has  been  added. 

d)  The  output  has  been  modified  to  be  more  pertinent  to  the 
average  user  by,  first,  limiting  output  to  values  of  more 
general  interest  and  by,  second,  providing  an  option  to  ob- 
tain only  an  abbreviated  output  of  the  critical  parameters  if 
one  so  desires. 
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e)  The  neck  torque  algorithm  has  been  modified  to  include  a 
resisting  torque  during  the  unloading  phase. 

f)  The  "SPRING"  subroutine  has  been  modified  to  prevent  certain 
subscripts  from  becoming  too  large. 

g)  Errors  found  in  the  expressions  for  head  and  chest  g’s 
have  been  fixed. 

h)  The  meaning  of  the  input  values  "Xl"  and  "Yl"  have  been 
changed  from  the  coordinates  of  the  airbag  center  to  the 
inf la tor  center.  This  was  done  to  prevent  the  user  from 
from  having  to  change  these  values  in  the  input  file  each 
time  a new  airbag  shape  was  simulated.  The  program  then 
makes  the  calculation  to  move  these  coordinates  to  the  air- 
bag center  without  user  involvement. 

i)  The  number  of  points  that  may  be  specified  for  the  crash 
pulse  has  been  increased  from  30  to  50. 

j)  The  neck  tensile  force  and  moment  computation  was  removed 
from  the  program  due  to  lack  of  data  to  which  these  com- 
putations be  be  checked  for  accuracy. 

k)  The  knee  trajectory  equations  were  modified  to  allow  the 
tibia  angle  to  increase  beyond  ninety  degrees  (the  vertical 
position)  without  the  algorithm  "blowing  up". 

In  addition  to  these  modifications,  a new  validation  of  the  program 
has  been  included  in  the  manual  as  an  example  case.  The  manual  also 
has  been  extensively  rewritten  so  the  latest  information  is  included. 

Michael  Fitzpatrick 
September,  1982 
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1.0  Introduction 


This  manual  is  written  to  give  the  potential  user  of  the  "PAC" 
computer  program  the  specific  information  he  will  need  to: 

a)  Set  up  the  data  input  file 

b)  Run  the  program 

c)  Interpret  the  results 

Prior  to  discussing  these  main  items  however,  let  us  present 
some  background  information  on  the  program. 

PAC  is  an  acronym  for  "Passenger  Air  Cushion".  As  this  title  in- 
dicates, the  program  was  written  to  describe  the  interaction  be- 
tween the  passenger  of  the  vehicle  and  his  airbag  in  a crash  sit- 
uation. Other  programs  have  been  written  in  the  past  to  describe 
such  an  interaction,  but  none  incorporate  the  combination  of  use- 
ful features  of  this  program  which  are  described  below. 

a)  Simulates  the  entire  deployment  sequence  of  the  unfurling 
airbag  from  its  stowed  condition  through  the  rebound  of  the 
passenger  away  from  the  fully  deployed  bag  at  the  conclusion 
of  the  crash  event. 

b)  Has  the  flexibility  of  simulating  virtually  any  normal,  pass 
enger  airbag  shape  from  a circular  cylinder  airbag  with  hemi 
spherical  ends  to  an  ellipsoidal  cylinder  airbag  with  ellip- 
soidal ends. 

c)  Has  the  versatility  of  specifying  a variable  airbag  deploy- 
ment angle  and  a variable  airbag  fabric  weight. 

d)  Has  the  flexibility  of  the  user  being  able  to  specify  the 
up-down  and  fore-aft  location  of  the  gas  generator  relative 
to  the  seated  passenger. 
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e)  Is  simple  ands  therefore,  inexpensive  enough  to  run  on  a 
small  in-house  computer  with  no  "library  routines"  from  an 
external  source  required  for  execution.  Therefore  the  pro- 
gram is  "self  contained"  as  well  as  inexpensive  and  simple 
to  operate  so  it  can  be  efficiently  used  as  a design  tool. 

f)  As  a design  tool,  the  program  is  oriented  to  the  user 
requirements  of  a typical  restraint  systems  engineer  with 
both  the  formulation  and  the  input  and  output  in  units  comm- 
only used  and/or  measured  in  a typical  test  situation. 

g)  As  a design  tool,  the  program  is  also  oriented  toward  the 
test  hardware  actually  encountered  in  most  situations.  For 
example,  past  computer  programs  might  model  the  passenger 
anthropometric  properties  very  well  but  might  neglect  the 
bag  shape  actually  used  and/or  the  bag  deployment  forces  that 
can  be  all  important  in  the  out-of-position  child  situation. 

h)  Finally,  the  program  is  balanced  so  that  the  various  comp- 
onents that  comprise  a normal  restraint  system  are  modeled 
to  be  of  approximately  equal  detail,  complexity  and 
accuracy. 

With  these  features  in  mind,  we  will  now  discuss  the  PAC  pro- 
gram in  some  detail. 
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2 . 0 Program  Description 


PAC  is  a two-dimensional,  lumped  mass  computer  model  of  a veh- 
icle passenger  interacting  with  an  initially  deploying  or 
already  deployed  airbag  of  arbitrary  shape  in  a crash  situation. 
The  model  includes  six  primary  masses  - four  to  degribe  the  pass- 
enger and  two  to  describe  the  airbag.  The  masses  that  arc?  used  to 
describe  the  passenger  are  the  head  mass,  the  main  torso  mas's,  the 
mass  of  the  sternum  (especially  useful  in  those  cases  where  "bag'-., 
slap"  is  of  interest),  and  the  lower  body  mass  (legs  and  hips). 

The  airbag  masses  are  composed  of  the  mass  that  impacts  the  pass- 
enger and  the  mass  that  surrounds  the  passenger  and  is  relatively 
free  to  expand.  These  two  airbag  masses  are  known  as  the  restrain- 
ed airbag  mass  and  the  unrestrained  airbag  bag  mass  respectively. 

The  passenger  airbag  is  simulated  by  an  ellipsoidal  cylinder  with 
ends  of  arbitrary  curvature  into  which  a programmed  amount  of  gas 
flows.  By  adjusting  the  airbag  vent  size,  a selected  amount  of  gas 
can  be  vented  during  the  simulated  crash  in  order  to  attenuate 
peak  chest  g’s  and  rebound  effects  thereby  aiding  in  the  selection 
of  an  optimum  airbag  design  for  a given  impact  situation. 

In  addition  to  the  basic  model  characteristics,  the  model  has 
specific  features  that  enable  the  user  to  ascertain  how  well  a 
basic  airbag  design  will  protect  the  out-of-position  child 
which  may  be  sitting  very  close  to  the  airbag  at  the  instant 
of  airbag  actuation.  In  order  to  make  this  determination,  the  com- 
plete airbag  deployment  process  had  to  be  modeled.  In  order  to 
lend  as  much  flexibility  to  this  feature  as  possible,  we  have 
made  it  possible  for  the  user  to  specify  in  detail  the  airbag  de- 
ployment geometry.  Thus,  the  airbag  deployment  angle  and  the 
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up-down  and  fore-aft  location  of  the  gas  generator  relative  to 
the  passenger  are  input  variables  specified  by  the  user. 

Specific  output  tailored  to  the  out-of-position  child  are  the  bag- 
slap  forces  and  the  chest  and  sternal  deflections,  velocities 
and  accelerations  induced  by  the  deploying  airbag. 

In  addition  to  these  basic  model  characteristics,  neck  rotational 
resistance,  seat  friction,  force-deflection  and  damping  properties 
of  the  chest  and  sternum,  and  knee  restraint  force-deflection  prop- 
erties are  additional  variables  that  need  to  be  specified  prior  to 
running  the  program. 

2 . 1 Mathematical  Formulation 

The  mathematical  formulation  of  the  eguations  of  motion  follows 
the  classical  Lagrangian  derivation  (Appendix  A)  with  body  pivot 
points  at  A and  B of  Figure  1 . 

The  lower  body  mass,  consisting  of  hips  and  legs,  is  constrained 
to  move  horizontally.  The  torso  mass  and  the  head  mass  are  free 
to  both  translate  and  rotate  about  the  H-point  and  neck  pivots. 

The  sternal  mass  is  constrained  to  move  in  a direction  perpend- 
icular to  the  chest. 

PAC  uses  a fixed  time  step  integration  routine  to  solve  the  diff- 
erential eguations  of  motion  numerically.  The  integration  routine 
chosen  was  the  Adams-Moulton  predictor-corrector  method  with  the 
fourth  order  Runge-Kutta  method  employed  to  determine  the  first 
four  solution  points. 

PAC  has  been  written  in  Fortran  IV  and  set  up  to  operate  in  an 
interactive  time  share  mode  for  maximum  utility  to  the  design  eng- 
ineer. The  program  is  self  contained  in  that  no  external  routines 
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Schematic  of  "PAC"  Computer  Model 


Figure  1 . 
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are  required  for  execution.  The  program  is  also  set  up  in  modular 
format  to  facilitate  the  addition  of  other  subroutines  at  a future 
time . 

Data  input  is  read  by  PAC  from  a previously  created  file  which  is 
then  displayed  and  described  as  a portion  of  the  total  output  im- 
mediately before  the  main  program  output. 

As  previously  mentioned,  PAC  has  been  programmed  in  modular  format 
with  several  subroutines.  This  enables  the  program  to  exhibit 
flexibility  in  two  important  ways.  First,  additional  algorithms 
may  be  added  if  desired.  Second,  it  makes  for  an  easy  way  to 
provide  for  tabular  data  input.  PAC  provides  for  this  in  two  ways. 

For  data  in  which  the  particular  value  of  the  dependent  variable 
is  a function  only  of  the  value  for  the  independent  variable,  a 
simple  table  look  up  and  interpolation  subroutine, "LOOKUP"  is  pro- 
vided. Inflator  gas  flow  versus  time,  vehicle  g's  versus  time, 
neck  torque  versus  angle,  and  chest  and  sternal  force  versus  def- 
lection are  examples  of  this  method  of  data  operation. 

However,  in  those  cases  where  the  dependent  variable  is  a function 
not  only  of  the  independent  variable,  but  also  depends  on  whether 
the  independent  variable  is  increasing  or  decreasing  , a different 
subroutine,  "SPRING",  that  allows  for  complex  plastic  behavior  is 
used.  In  this  case  one  must  not  only  specify  the  values  of  the 
dependent  variable  for  different  values  of  the  independent  var- 
iable, but  must  also  specify  the  "unload  slopes"  for  those  con- 
ditions in  which  the  member  is  undergoing  unloading  during  a less- 
ening of  the  degree  of  deformation.  Knee  restraint  force  as  a fun- 
ction of  crush  and  seat  friction  force  as  a function  of  stroke 
are  handled  by  this  subroutine. 
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2 . 2 Airbag 


The  airbag  shape  chosen  for  use  in  the  PAC  model  has  been  chosen 
to  be  of  as  general  a shape  as  possible  without  introducing  an 
airbag  algorithm  that  would  require  unrealistically  high  amounts 
of  run  time  on  the  computer  or  be  incompatible  with  the  degree  of 
complexity  of  the  remainder  of  the  program.  Figure  2 shows  the 
airbag  shape  chosen  for  the  PAC  model.  As  can  be  seen  from  the 
figure,  three  ellipsoidal  axes  dimensions  (a^,b^,  and  c^ ) plus 
the  half  lengh  of  the  cylindrical  portion  (A  ) are  necessary  to 
describe  the  general  airbag  shape. 

In  addition  to  these  variables,  three  more  are  required  to  describe 

the  shape  of  the  bag  in  area  where  the  bag  encloses  the  inflator. 

These  variables  are  the  inflator  diameter  (D.  ,_),  the  distance  the 

inf 

"manifold  sock"  extends  forward  of  the  dash  (D  ),  and  the  width 

ms 

of  the  manifold  sock  (wsock)  "which  is  equivalent  to  the  inflator  length. 

Initially  the  airbag  is  assumed  to  be  rolled  up  and  have  the  shape 
of  a very  small  ellipsoidal  cylinder  with  an  axis  length  c^  equal 
to  After  bag  inflation  begins,  the  bag  is  assumed  to 

inflate  along  all  three  airbag  axes  a proportionate  amount  so  that 
a constant  ratio  of  axis  lengths  exist  until  full  airbag  inflation 
is  obtained.  The  direction  the  bag  deploys  is  specified  by  assigning 
the  value  desired  to  0^  in  the  input  file. 

In  order  to  properly  describe  the  deployment  sequence,  it  is  nec- 
essary to  keep  track  of  two  separate  airbag  masses;  the  airbag  mass 
that  impinges  upon  the  torso  of  the  passenger  which  is  called  the 
"restrained"  airbag  mass,  and  the  airbag  mass  that  is  free  to  expand 
until  full  deployment  which  is  called  the  "unrestrained"  airbag  mass. 
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"PAC”  Airbag  Schematic 


Figure  2 


Appendix  B contains  the  derivations  of  key  equations  that  were  nec- 
essary to  describe  the  airbag  deployment  sequence  as  well  as  other 
airbag  related  items  such  as  derivation  of  wraparound  forces,  press- 
ure forces,  contact  areas,  airbag  volume,  etc. 

2 . 3 Knee  Restraint  and  Seat  Friction 


As  previously  mentioned,  the  program  is  set  up  to  accept  tabular 
input  for  the  force  versus  crush  properties  of  the  knee  restraint 
and  the  force  versus  displacement  properties  of  seat  friction.  It 
is  these  values  that  primarily  determine  the  kinematics  of  the  lower 
body.  Appendix  C contains  the  derivation  of  knee  restraint  equations. 

The  user  specifies  in  tabular  format  in  a data  file  what  these  prop- 
erties are  to  be.  In  addition,  the  user  specifies  the  "unload  slope" 
so  the  program  can  compute  the  path  the  unloading  sequence  is  to 
take  during  rebound  away  from  the  knee  bolster  or  movement  rearward 
over  the  passenger  seat.  Specifics  on  how  this  and  all  other  data 
input  is  handled  is  discussed  in  Section  3.1. 

2 . 4 The  Passenger 

The  passenger  is  modeled  by  four  masses,  the  head,  the  main  torso 
mass,  the  sternal  mass,  and  the  lower  body  mass.  Pivot  points  exist 
at  points  A and  B as  shown  in  Figure  3.  This  figure  also  describes 
the  overall  passenger  geometry  and  location  of  the  passenger  with 
respect  to  the  vehicle  interior. 

Specific  details  of  providing  the  passenger  related  input  into  the 
data  file  is  discussed  in  Section  3.1. 

The  torso  and  sternal  masses  are  assumed  to  be  connected  in  series 
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Inf la tor 
Center 


Figure  3. 
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as  shown  below.  The  spring  forces  as  a function  of  the  deflection 
of  the  spring  are  input  in  tabular  format.  The  damping  forces  are 
assumed  to  be  a function  of  velocity  and  are  computed  by  multiply- 
ing the  relative  velocity  between  the  bag  mass  and  the  sternal  mass 
and  the  sternal  mass  and  the  main  torso  mass  by  an  appropriate  fac- 
tor included  as  part  of  the  data  input.  These  factors,  one  for  each 
mass,  are  known  as  the  damping  coefficients. 


A further  comment  is  necessary  in  regard  to  the  resisting  torgue 
generated  by  neck  muscular  resistance  and  anatomical  interferences 
due  to  relative  displacement  between  the  head  and  the  torso  since 
these  forces  are  only  applied  if  certain  conditions  are  met. 

In  those  cases  where  the  head  is  returning  to  be  more  nearly  in 
line  with  the  passenger's  torso  (unloading),  the  neck  torgue 
resisting  head/neck  relative  movement  is  reduced  to  0.5  times  the 
"loading"  value  for  hyperextension  and  0.4  times  this  value  for 
hyper f 1 exi on . This  has  been  done  to  make  the  neck  unloading  con- 
form to  recommended  performance  parameters  for  mechanical  necks. 

In  addition  to  the  neck  resisting  torgue  based  upon  angular  dis- 
placement of  the  head  relative  to  the  torso,  we  have  also  included 


— Chest  - 


- — Sternum 


Chest  and  Sternum  Dampers 


Chest  and  Sternum  Springs 
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a damping  coefficient  for  the  head  based  based  upon  the  angular 
velocity  of  the  head  relative  to  the  torso.  This  value,  known  as 
DCN  in  the  input  data  file,  is  multiplied  by  this  angular  velocity 
and  added  to  the  neck  torque  calculated  as  a function  of  the  ang- 
ular displacement  of  the  head  relative  to  the  torso  for  the  total 
neck  torque  resisting  head/torso  relative  motion. 

For  an  overall  listing  of  the  PAC  program,  please  refer  to 
Appendix  D. 
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3 . 0 Example 


The  best  way  to  illustrate  the  use  of  the  PAC  model  is  by  an 
example.  For  our  example  we  chose  to  simulate  a recent  sled 
test  supervised  by  Fitzpatrick  Engineering  and  conducted  at 
Transportation  Research  Center  of  East  Liberty,  Ohio.  This 
test  was  conducted  as  part  of  Contract  No.  DTNH22-81 -C-071 32, 

"A  Systems  Analysis  Approach  into  Integrating  Airbags  into  A 
Production  Ready  Small  Car".  In  this  program,  Fitzpatrick  Eng- 
ineering designed  air  cushion  restraint  systems  for  the  driver 
and  passenger  of  the  DeLorean  sports  car.  By  simulating  an 
actual  sled  test  we  will  not  only  be  able  to  illustrate  the  use 
of  PAC  by  example,  but  we  will  also  be  able  to  compare  test 
results  with  computer  results  giving  some  indiction  as  to  the 
degree  to  which  the  program  is  able  to  predict  what  will  happen 
in  an  actual  crash  situation. 

The  test  case  chosen  was  TRC  sled  run  no.  7 of  the  above  ref- 
erenced contract.  This  test  was  a frontal  test  at  34  mph  with 
an  airbag  sensing  time  of  15  msec.  The  Part  572,  50t*1  percentile 
male  ATD  was  seated  normally  in  the  DeLorean  passenger  seat. 

The  gas  flow  rate  into  the  airbag,  appropriate  airbag  dimensions, 
airbag  deployment  geometry  and  other  information  needed  to 
specify  the  information  required  for  computer  input  were  obtain- 
ed from  the  appropriate  source  and  prepared  for  computer  input. 

3 . 1 Creation  of  the  Data  Input  File 

The  first  thing  one  must  do  in  preparation  for  making  a computer 
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run  is  to  set  up  the  data  input  file.  This  file  is  a listing  in 
a pre-ordained  format  of  the  information  the  program  reguires  in 
order  to  run  the  given  case.  Pages  15  through  18  list  the  para- 
meters for  the  data  file  for  this  particular  case  as  well  as  show- 
ing how  to  set  up  the  data  file  for  any  other  case.  We  will  now 
discuss  how  the  information  presented  on  pages  15  through  18  relates 
to  the  data  file. 

On  page  15  the  first  column  shows  the  location  of  a particular 
piece  of  data  in  the  input  file.  For  example  under  "LOCATION",  the 
"LINE  NO ."  corresponds  to  the  line  number  in  the  data  file.  Line 
number  "1"  would  correspond  to  the  first  line  of  the  data  file, 
line  number  "2"  the  second  line,  etc. 

The  second  number  under  the  "LOCATION"  heading  is  the  "LINE  LOC . " 
which  is  the  location  in  the  line  of  a particular  piece  of  data. 

For  example,  a "2"  in  this  column  indicates  the  second  piece  of  data 
in  that  particular  line. 

In  the  second  main  column  we  encounter  a column  heading, "NAME/UNI TS" 
The  "NAME"  is  the  name  of  the  particular  piece  of  data  referred  to 
in  the  first  column.  Further,  this  name  is  as  programmed  in  the  PAC 
program  and  is,  therefore,  alpha-numeric  in  character.  The  "UNITS" 
part  of  the  second  column  gives  the  units  in  which  the  particular 
data  piece  must  be  entered  in  the  file. 

In  the  third  main  column  a short  description  of  the  data  in  the  part 
icular  file  location  is  given. 

In  column  four,  the  actual  value  used  in  the  file  for  this  particula 
example  is  given. 
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15 


* 

Notes  The  second  value  entered  on  this  line  must  be  the 
sensing  time. 
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17 
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In  summary,  consider  the  first  line  on  page  15.  In  this  line  the 
first  column  tells  us  that  the  information  being  described  is  en- 
tered on  the  fi rst  line  of  the  data  file  in  the  f i rst  position  in 
the  line.  The  second  main  column  tells  us  that  the  piece  of  data 
that  appears  in  this  first  line  in  the  first  position  is  known  as 
"Y(4)M  in  the  program  and  has  units  of  "MPH"  in  the  file.  The 
third  column  tells  us  that  this  variable  is  the  "Vehicle  Impact 
Velocity"  while  the  fourth  and  final  column  tells  us  that  its  value 
in  this  particular  example  is  "34.0"  mph. 

The  data  file  that  results  from  the  information  given  on  pages  15 
through  18  is  shown  by  Figure  4.  We  have  chosen  to  call  this  file 
"TRCST7P" . 


3 . 2 Running  the  Program 


Once  the  input  file  has  been  created  and  saved,  the  program  is 
ready  to  be  run.  At  this  point  the  user  accesses  PAC  and  tells  the 
computer  to  run  the  program.  The  computer  will  respond  by  asking 
the  user  to  name  the  data  input  file.  In  this  case  we  would  respond 
by  entering  "TRCST7P"  from  the  keyboard,  as  shown  by  the  sample  run 
shown  as  Appendix  E.  Next  the  user  is  prompted  for  a choice  on  whether 
he  wishes  the  full  or  abbreviated  list  of  output.  Assuming,  a full 
list  is  desired,  the  user  types  ”1"  . (if  the  abbreviated  list  were 
selected  by  typing  "2",  the  output  would  be  as  shown  in  the  sec- 
ond list  of  output  in  Appendix  E) . Once  these  answers  are  given, 
the  computer  will  print  out  the  input  data  followed  by  the  selected 
output.  We  will  now  discuss  this  output. 

Altogether  there  are  ten  separate  blocks  of  output,  most  with  seven 
pieces  of  information  presented.  In  each  block  of  output,  the  first 
piece  of  information  is  always  the  same  and  is  the  elapsed  time  from 
the  beginining  of  impact  in  msec. 
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The  first  block  of  output  shown  in  Appendix  E consists  of  the 
following  items  as  we  read  from  left  to  right  across  the  page* 
elapsed  time  (TIME),  vehicle  g's  (VEH  G'S),  vehicle  velocity 
(VEH  VEL) , vehicle  crush  (VEH  DI SP ) , amount  of  bag  penetration 
by  the  chest  (CHEST  BP),  chest  wraparound  force  (CWA  FORCE),  and 
chest  pressure  force  (CPR  FORCE) . The  units  for  these  variables 
are  indicated  immediately  below  the  headings. 

The  second  data  block  contains:  elapsed  time  (TIME),  overground 
displacement  of  the  H-point  (H-P  DISP),  the  overground  velocity 
of  the  H-point  (H-P  VEL) , the  force  imparted  to  the  passenger 
through  seat  friction  in  sliding  over  the  seat  (SEAT  FR.),  the 
force  imparted  to  each  femur  by  the  knee  restraint  (FEM  FORCE), 
the  angle  of  the  femur  from  horizontal  (see  Figure  1)  (KNEE  ANG) , 
and  the  angle  of  the  tibia  from  horizontal  (see  Figure  1).  This 
parameter  is  shown  as  "TIB  ANG"  in  the  output. 

The  third  data  block  contains:  elapsed  time  (TIME),  the  overground 
horizontal  displacement  of  the  torso  center-of-gravi ty  (TORSO  DISP), 
the  angle  of  torso  inclination  with  respect  to  a vertical  line  as 
shown  in  Figure  3 (TORSO  ANG),  the  angular  velocity  of  the  torso 
(TORSO  VEL),  the  angular  acceleration  of  the  torso  (TORSO  ACC), 
the  displacement  of  the  torso  in  a horizontal  direction  with 
respect  to  the  vehicle  compartment  (TORSO  R.D.),  and  the  velocity 
of  the  torso  in  a horizontal  direction  with  respect  to  the  vehicle 
compartment  (TORSO  R.V.). 

The  fourth  data  block  contains  the  exact  egui valent  to  the  third 
data  block  with  the  exception  that  the  data  is  for  the  head  in- 
stead of  the  torso  and  the  last  item  in  the  data  block  is  the  head 
angle  relative  to  the  torso  .(HEAD  R.  ANG)  as  shown  in  Figure  1 . 
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The  fifth  data  block  contains:  elapsed  time  (TIME),  the  restrained 
bag  (part  of  bag  impinging  on  chest)  acceleration  (R  BAG  ACC),  the 
restrained  bag  velocity  with  respect  to  the  ground  in  the  horiz- 
ontal direction  (RBV  WR  GND),  the  restrained  bag  velocity  with 
respect  to  the  chest  in  a direction  normal  to  the  chest  (RBV  WR  CST) , 
the  restrained  bag  velocity  with  respect  to  the  dashboard  in  the 
direction  of  airbag  deployment  (RBV  WR  DSH) , the  restrained  bag 
displacement  relative  to  the  ground  in  a horizontal  direction 
( RBD  WR  GND),  and  the  restrained  bag  displacement  with  respect 
to  the  dashboard  in  the  direction  of  airbag  deployment  (RBD  WR  DSH). 

The  information  contained  in  the  sixth  data  block  is  exactly  equiv- 
alent to  that  presented  in  the  fifth  data  block  with  the  except  Lon 
that  the  information  pertains  to  the  unrestrained  portion  of  the 
airbag  (that  part  not  in  contact  with  the  passenger). 

The  seventh  data  block  contains:  elapsed  time  (TIME),  the  force 
applied  to  the  chest  due  to  "bagslap"  (CST  F BSP),  the  force  applied 
to  the  sternum  due  to  bagslap  (STN  F BSP),  the  velocity  of  the 
sternum  with  respect  to  the  chest  in  a direction  normal  to  the 
chest  surface  (STV  WR  CST),  the  deflection  of  the  sternum  by  the 
impacting  bag  roll  (RLD  WR  STN),  the  deflection  of  the  sternum 
with  respect  to  the  chest  (STD  WR  CST),  and  the  distance  from  the 
aft  edge  of  the  inf la tor  to  the  torso  in  the  bag  deployment  dir- 
ection (DTORSO) . 

The  eighth  data  block  contains:  elapsed  time  (TIME),  the  amount  of 
bag  penetration  by  the  head  (HEAD  BP.),  the  airbag  volume  (BAG  VOL.), 
the  airbag  pressure  (BAG  PRESS.),  the  head  wraparound  force  (IIW/A 
Force),  the  head  pressure  force  (HP  FORCE),  and  the  volume  of  the 
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airbag  intercepted  by  the  chest  and  head  (INT.  VOL) . 

The  ninth  data  block  contains:  elapsed  time  (TIME),  the  chest  A-P 
(anterior-posterior)  g's  (CHEST  AP),  the  chest  S-I  (superior-inf- 
erior) g* s (CHEST  Si),  the  head  A-P  g's  (HEAD  AP),  and  the  head 
S-I  g’ s (HEAD  SI ) . 


The  tenth  and  final  data  block  presents  the  following  data:  elapsed 
time  (TIME),  the  diameter  of  the  bag  roll  (ROLL  DIA),  the  distance 
from  the  H-point  to  the  point  where  the  bag  roll  impacts  the  chest 
(ROLL  RAD),  the  X coordinate  of  the  airbag  center  from  the  chosen 
reference  point  (XC  B CTR) , the  Y coordinate  of  the  airbag  center 
from  the  chosen  reference  point  (YC  B CTR),  the  weight  of  the  rest- 
rained portion  of  the  airbag  (WRB),  and  the  weight  of  the  unrest- 
rained portion  of  the  airbag  (WURB) . 

At  the  end  of  the  output,  the  user  is  asked  whether  he  wishes  a HIC 
computation.  If  so,  the  user  types  "1"  and  the  HIC  is  computed  along 
with  the  beginining  and  ending  times  of  the  computation. 

Sign  Convention 

Now  that  the  input  file  and  the  output  values  have  been  described 
in  some  detail,  it  is  necessary  to  discuss  the  sign  convention 
which  has  been  used  in  setting  up  these  groupings  of  data.  The 
input  file  will  be  discussed  first  followed  by  the  output  file. 

Input 


The  sign  convention  for  the  data  in  the  input  file  in  Figure  4 
is  to  be  generally  positive  except  for  the  following  cases 
where  specific  signs  must  be  given  to  avoid  confusion. 

Y( 6 ) ; Head  Angle;  Line  1,  Location  2;  Positive  as  shown  in  Fig.l. 
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Y(7);  Torso  Angle;  Line  1,  Location  3;  Positive  as  shown  in  Fig.  1. 


FNECK(1,K);  Neck  Torque  Relative  Angle  (Y(6)-Y(7));  Line  14,  All; 
Positive  when  head  pitched  forward  relative  to  chest  as  shown  in 
Fig . 1 . 

FNECK ( 2 , K) ; Neck  Torque;  Line  15,  ALL;  Positive  when  in  direction 
of  positive  Y(6) . 

CF(1,K);  Chest  Deflection  due  to  bagslap;  Line  22,  All;  Positive 
when  chest  compressed. 

CF(2,K);  Chest  Force  due  to  bagslap;  Line  23,  All;  Positive  when 
producing  chest  compression. 

This  completes  the  discussion  on  sign  convention  for  the  input  file. 
We  will  now  discuss  the  sign  convention  used  in  presenting  the  out- 
put values  presented  in  Appendix  E. 

Output 

Displacements,  velocities  and  accelerations  are  considered  pos- 
itive when  they  are  in  either  the  upward  direction  in  the  ver- 
tical mode  or  in  the  direction  of  original  vehicle  travel  in 
the  horizontal  mode. 

When  the  direction  of  displacement,  velocity  or  acceleration  is 
neither  horizontal  nor  vertical  (such  as  the  restrained  bag 
diameter  with  respect  to  the  dash,  RBD  WR  DSH,  which  is  in  the 
bag  deployment  direction)  the  sign  given  in  the  output  is  the 
sign  that  corresponds  to  the  direction  of  its  horizontal  comp- 
onent. Thus  a negative  value  for  "RBD  WR  DSH"  mean  that  the 


24 


restrained  bag  displacement  with  respect  to  the  dash  is  in  a 
direction  opposite  to  the  direction  of  original  vehicle  travel; 
i.e.  toward  the  passenger. 

In  the  case  of  forces  and  moments  applied  to  the  passenger  in 
which  the  direction  of  application  is  obvious,  the  values  listed 
in  the  output  are  listed  as  positive  values.  For  example,  the 
following  forces;  the  CWA  FORCE,  the  CPR  FORCE,  the  HW/A  FORCE, 
the  H P FORCE,  FSTBS , FCBS,  the  FEM  FORCE  and  the  SEAT  FR . force 
are  listed  in  the  output  as  positive  values  and  are  understood 
to  be  in  a direction  opposing  further  airbag  penetration,  knee 
restraint  penetration  and  sliding  across  the  seat  when  applied 
to  the  passenger. 


25 


3 . 3 Comparing  Computer  Results  with  Test  Results 


Test  conditions  for  TRC  Sled  Test  No.  7 were  duplicated  on  the 
computer  as  accurately  as  possible  and  then  a simulation  of  this 
test  was  made  to  verify  the  ability  of  PAC  to  reproduce  actual 
test  results.  Figure  5 shows  the  sled  test  "crash"  pulse  and  the 
data  points  used  for  input.  Other  data  needed  for  input  were  ob- 
tained from  pre-test  measurements,  "known"  dummy  properties  and 
gas  flow  information  provided  by  Thiokol,  the  gas  generator 
manufacturer . 

Figures  6 through  10  show  the  comparison  of  computer  results  and 
test  results  (Appendix  E shows  the  actual  computer  run) . As  is 
evident  from  the  figures,  very  good  correlation  was  obtained. 

Figure  11  shows  the  passenger  trajectory  as  computed  and  then 
plotted  by  computer.  Comparison  of  this  trajectory  with  the 
results  actually  observed  in  the  test  show  very  similar  dummy 
movement.  The  only  difference  noticed  was  that  in  the  test  the 
dummy  H-point  appeared  to  be  approximately  two  inches  lower  than 
was  measured  before  the  test.  This  was  perhaps  due  to  seat  com- 
pression which  is  not  computed  in  the  PAC  simulation.  Thus  the 
head  is  shown  very  close  to  the  windshield  in  Figure  11  whereas 
in  actuality  is  was  not  quite  so  close. 

The  table  on  the  following  page  compares  the  injury  measures 
actually  measured  in  the  test  with  those  predicted  by  the  com- 
puter . 


26 


SLED  TEST  NO.  7 


Injury  Measure 


Test  Results 


Computer  Results 


HIC 


606 


619 


Pk.  Res.  Chest  G’s 
(-3  msec) 


50 


50 


Pk.  Femur  Loads 

(Lb:  Righr,  Left)  1070,  1140 


1097,  1097 
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REST , MU  DEL  VALID.  TEST 


r-  a 


co 

CO 

CO 


t-ij 


LiZi 

CO 

<=> 

>70 


'X: 


=*•  co 

c=i  rr> 
c-.j  i — 

CO  Ll_ 


D®  ' 


s 

■S3 

ur" 

CO 


28 


FLOOR  PAN  TUNNEL  ACCELERATION 


flEbT . MuDEL  v'fiLID.  TEoT 
a 2 ®7 4 FILTER  5 fiLFF 

HEDRG2  MIN,  MAX  VRLUEB 


r=> 


LO 


CL- 

UJ 

i_5) 


CO 

Ul 

OZ 


i T 


LLJ 


! ) 
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S5 1075050  , TRC862  31-KHh-a^  1 

REST.  MODEL  VALID.  TEST 
82374  FI  HER  : BLFF 

CS1RG2  MIN,  HR X VALUE 


sa 

<S» 

lO 


30 


C H E S T fl  C C ELBE  R T ] 0 N Fi  E 3 U L T R N T P R S 5 E N G E R 


S5 1075050  , T Pi C 06 2 3 1 - MRR A„.  iu:  uu 

REST.  MODEL  VRL.ID.  TEST 

82074  FILTER  = BLPF  1000/  3170/  -40 

LFMF2  MIN » MAX  VALUES  = *377.18  & 23.83,  1080.89 


<=> 

to 
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LEFT  FEMUR  FORCE  PASSENGER 


S51075O50  , TRC862  3i -NfiR-82  13:10:10 

REST . MODEL  VALID.  TEST 

82074  FI  HER  = BLPF  1000/  3170/  -40 

RFNF2  HIM,  M Fl X VALUES  = -734.29  e 23.75,  1142 


S3 

(=> 

in 
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RIGHT  FEMUR  FORCE  PASSENGER 


REST.  MODEL  VRLID.  TEST 

32074  FILTER  = BLPF  300/  950/  -40 

flBF‘2  MIN.  MUX  VALUES  = * 1.313  » 37.25.  15.05  e 18,50 


<s> 

IS* 


lO 


LU 
r_ £-i 

LU 

CO 

CO 


LiJ 

Q_ 

i 

CO 

CO 

LU 

rt 


f - j 

cn 

CQ 
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COMPUTER  SIMULATION 


OF : 


DELOBEAN  SLED  TEST  HO.  7 - 50TH  PERCENTILE  PASSENGER  - 34  MPH 

HIC  : 619 

PR.  CHEST  G'S  = so  g 
PE.  FEMUR  LOAD  : 1097  lb 


Figure  11 . 


3 . 4 Helpful  Hints 


In  the  following  list  we  will  leave  the  user  with  a few  "helpful 
hints"  which  should  assist  the  user  of  the  PAC  model  in  obtaining 
the  most  accurate  results  possible. 

1.  The  integration  interval  (STEP)  should  not  be  greater  than 
.001  sec  for  adult  passengers  nor  greater  than  .0005  sec 
for  child  passengers.  Further,  for  sternal  masses  less  than 
1.5  lbs,  the  integration  interval  should  not  exceed  .0002 
sec . 

2.  The  chest  should  not  be  placed  closer  than  one  inch  to 
the  dash. 

3.  As  presently  set  up,  the  H-point  is  assumed  to  be  the  lower 
body  pivot  point.  It  is  assumed  that  the  airbag  contact  line 
passes  through  this  point  and  the  neck  pivot.  For  greatest 
simulation  accuracy,  the  effective  pivot  points  should  be 
moved  foreward  slightly  (approximately  one-half  the  torso 
thickness)  so  that  the  line  corresponds  with  the  front 
surface  of  the  torso  and  not  the  torso  centerline. 

4.  For  airbags  which  are  unsymmetrical  about  a horizontal  line, 
the  user  should  account  for  this  in  the  simulation  by  raising 
or  lowering  the  effective  mounting  point  of  the  airbag  an 
appropriate  amount. 

5.  For  airbags  which  deviate  significantly  from  the  ellipsoidal 
cylinder  shape  assumed  in  the  PAC  derivation,  the  user  should 
approximate  the  actual  airbag  by  an  ellipsoidal  cylinder  of 
most  similar  shape. 
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6.  As  presently  derived*  the  PAC  model  assumes  a seated  pass- 
enger. If,  for  example,  one  wishes  to  simulate  a standing 
child,  the  "RT"  dimension  must  be  modified  to  be  the  distance 
from  the  H-point  to  the  "effective”  C.G.  of  the  upper  and 
lower  body.  This  roughly  accounts  for  the  mass  of  the  legs 
distributed  below  the  H-point  which  is  significant  for  the 
standing  child.  Past  experience  has  shown  this  value  for 
R to  be  much  smaller  than  the  normal  value  for  the  seated 
chi Id . 
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APPENDIX  A 


Derivation  of  the  Equations  of  Motion  for  the  Passenger 


APPENDIX  A 


DERIVATION  OF  THE  EQUATIONS  OF  MOTION 

The  derivation  of  the  equations  of  motion  will  be  formulated  utilizing  Lagran- 
gian  techniques  based  upon  the  geometrical  representation  in  Figure  A-l. 

Writing  an  expression  for  the  total  kinetic  energy  of  the  occupant , we  have: 

(1)  T = 15  ^ + Mj.  (X^  + Yp  + M^] 

Note  that  Y,  5 0,  as  no  movement  normal  to  the  X-direction  is  allowed  for  the 
hip-leg  mass. 

Mpj  = Head  mass 

Mp  = Torso  mass 

Mp  = Hip- leg  mass 

Xp  = Horizontal  translation  of  the  hip -leg  mass  with  respect  to  inertial 
reference  point  -which  is  positive  when  it  is  in  direction  shown. 

Xp  and  Xpj  are  similarly  defined 

Ypj  = Vertical  distance  from  H-point  o the  center  of  gravity  of  the  head 
Yp,  = Vertical  distance  from  H-point  to  the  center  of  gravity  of  the  torso 
Successive  dots  indicate  velocity  and  acceleration,  respectively. 

Writing  the  transformation  equations,  we  have: 


(2 ) 

V 

*L  + rT  SIN9T 

(3) 

hT 

ii 

rT  COSQp, 

(4) 

n 

+ rN  SIN0I  + rH  SIN9H 

(5) 

V 

rN  COS0T  + rH  C0S6H 

(6) 

ii 

*L  * rT  C0S6T  9T 

(7) 

V 

-r_  sdm  e_ 
T T T 

A 
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MATHEMATICAL  MODEL 


Figure  A-l.  Passenger/Airbag  interaction 

Schematic 
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(S)  Xn  - \ rN  coseT  e?  ♦ rH  coseH  eH 

(9)  Yh  - -rN  SIN6TeT  - rH  SIN0H 
where : 

rT  = Distance  from  hip  H -point  to  torso  center  of  gravity  (The  H- point 
is  assumed  to  be  coincident  with  the  hip-leg  center  of  gravity.) 

r^  = Distance  from  H-point  to  neck  pivot  point 

rH  = Distance  from  the  neck  pivot  point  to  the  center  of  gravity  of  the 
H head 

0^  and  0y  are  as  defined  in  Figure  A-l„ 

Substituting  Equations  6 through  9 into  Equation  1,  we  have: 

(10)  I - ■=  {mjX*  + Kj.  [x*  + 2 C0S6T  0T  ♦ r*  8*]  * Mh  [x*  * 

2 XL  (rN  C0SeT  8t  ♦ rH  C0S9H  0^  ♦ 2 rN  rH  (COS0T  COS0H  0?  ^ ♦ 
SIN9t  S!N0h  9t  0h)  ♦ r>  6|*r*0*j} 

The  potential  energy  portion  of  the  Lagrangian  is: 

(11)  VT  - Mj.  g rT  COS0T 

(12)  VH  = ^ * <rH  C0S6H  * rN  W 

Note:  The  applied  forces  and  moments  will  be  treated  separately  later  on. 

Writing  the  Lagrangian,  we  have: 

(13)  ' L = T-V  = T - (VT  + VH) , 

where  the  values  to  be  substituted  into  this  equation  are  given  by  Equations  10 
11  and  12. 
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The  basic  equation  in  Lagrangiai  mechanics  is: 


(14) 


d /9L_\  _ 3L  _ 
dt  \ 3q^  / qi 


where : 

q^  = generalized  displacement  of  the  i^1  mass 
= generalized  velocity  of  the  i1^1  mass 
Fqi  = generalized  force  acting  on  the  i^1  mass 

Taking  the  required  derivatives  from  Equation  13  for  substitution  into  Equa- 
te14,  we  obtain: 


(15)  3L 


= (^  ♦ ^ ♦ y ^ ♦ M,.  rT 


coseT  eT  + crN  coseTex  + iH  coseH  y 


(10)  d_  / 3L_  \ =(ML  ♦ Mj.  + Xj_  - MjXj,  SIN6t  0t  - (rN  SIN8T  6*  + 


rH  SIN9h  6*)  ♦ ^ rT  C0S9T  ^ ^ (rN  COS0T  ^ ♦ 

rH  COS0H  y 


(17) 

3L 

3Xl 

= 0 

(18) 

3L 

a a. 

1! 

4s 

0?- 

hh 

8 

8 

H 

+ 

CD  • 

tL 

+ 

"h 

x.  rXI  cose^  + r r (cose^  cose  p 

L N i No  i 

T 

* SIN0T  SIN0H  y ♦ 0T' 
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(19)  5F  (#)  ■ "V  (*L  rT  C0SeT  - *L  rT  SIN9T  ®T  + rT  eT^ 


(20)  3L_ 

30„ 


* ^ |_  *L  rN  COS0T  ' *L  rN  SIN0T  9X  • rN  rH  (SIN9I  9T  C0S9H  9H 

+ eose„  siN9„  e2  - cose„  cose,,  e„  - cose„  L • siNe„  e„ 

T H n I Hn  ii  n n 


- SIN9t  CO S0U  0"  - SIN0  SIN0  0„)  + r*  0T 
T HH  T HH  NT 


= -Hj.  \ rT  SIN0t  9t  - ^ rN  (SIN0T  ©T  ^ + rH  SIN9T  COS0H 


0T  0H  - rH  COS0T  SIN0h  0T  0H)  + Mj.  g SIN0T  ^grN  SIN0T 


(21)  3L 


30 


H 


"h  rH  c°S0h 


rM  ru  (COS0™  COS0U  0„  + SING  SIN0  0 ) 
NH  T HT  T HT 


+ r2  0 
H H 


(22)  d_  / 3L 
dt 


(23)  3L 
30 


H 


/ 3L  \ • « • 

( 30“ ) = ^ \ rH  C0S9H  ’ \ rH  SIN0H  9H  " rN  rH  L°IN6T  COS0H  0T 
H L 

+ coseT  siNeH  eT  eH  - coseT  coseH  eT  - coseT  siNeH  e2 
- SIN0T  coseH  eT0H  - SIN0T  siNeH  eT)  + r2  eH 

^ rH  [SIN9H  0H  \ + rN  (COS0T  SIN0H  9I  9H  - SIN0T  COS0H  * 
9T  0h)]*  “h  8 rH  SIN9H 


Substituting  Equations  16  and  17  into  Equation  14,  we  have: 

(24)  * Mj,  * Mjj)  Xl  - M,.  rT  SIN0T  02  - ^ (rN  SIN0.J.02  + rH  SIN0H  02) 

* *r  rT  COS0T  eT  ♦ Mn  (rN  coseT  eT  * rH  coseH  5 ) = F;a 


which  is  the  equation  of  motion  for  mass  M^. 
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Substituting  Equations  19  and  20  into  Equation  14,  we  have: 


“r  (XL  rT  C0SaT  - \ rT  SIN0T  9T  ♦ r/  y ^ ^ rN  C0SST 


- *L  rN  SINeT  ®T  ' rN  rH  (SIN8T  ®T  C0S6H  ®H  + C0S6T  SIN6H  6 

- cos8t  cos6h  eH  - cose  eT  siNeH  e - siNe  coseH  e* 


H 


- SIN0t  SIN9h  0 ) + r>  0T 


+ Mr*Lri  SIN0  eT  ♦ ^ rN  ( 


SIN0t  0t  ^ + rH  SIN9t  COS0h  ©t  ©h  - rH  COS0T  SIN0H  0^  0H) 


Hj-g^  SIN0t  - g rN  SIN0T  = F 


0T 


Rewriting  the  above  yields: 


(25)  Mj.  (XL  rT  COS0T  + r^.  0T)  + ^ rN  COS0T  - rN  rR  (C0S©T  SIN0H  02 

- COS0T  coseH  0H  - SIN0T  coseH  e2  - sin©t  siNeH  0 ) 


+ r2  0 


N T 


• Hf  g rT  SIN0t  -^grN  SIN0T  = F©T 


which  is  the  equation  of  motion  of  the  torso  mass. 


Substituting  Equations  22  and  23  into  Equation  14,  we  have: 

^ [*L  rH  C0S6H  ' h rH  SIN0H  ®H  ' rN  rH  (SIN9T  C0S6H  9X 

* coseT  SIN0H  eT  eH  - coseT  cos6H  eT  - coseT  sinsh  e* 

- SIN0t  COS0H  0T  0H  - SIN0t  SIN0h  V ♦ r*  eHJ  ♦ ^ rH  [ 

SIN0h  0h  ^ + rN  (COS0T  SIN0h  0?  0R  - SIN0t  COS0h  ^ 0^ 


g r„  SIN0„  = F, 


H 0H 
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Rewriting  the  preceding  yields: 


(26)  M rH  COSeH  - rN  rH  (SIN8T  COS0H  9*  - COS0T  COS0H  3T 


- COS0T  SIN0h  0*  - SIN0t  SIN0h  0t)  + r*  0H 


“h  g rH  SINeH  = F0H  ’ 


which  is  the  equation  of  motion  for  the  head  mass , 


Writing  Equation  24  in  terms  of  XL,  we  Have: 


(27) 


*L 


1 - FXL  + rT  + ^ V Sm0T  9T  + % rH  eH  SIN3 


^ rT  + W 9I  COS0T  - "V  rH  0H  COS0H 


Writing  Equation  25  in  terms  of  8^,  we  have: 

(28)  “if 

eT  - ^ 4 ♦ % rfs  (F0T  - <*r  rT  * «H  V *L  CCS0T 


Mh  rN  rH  [ 0H  (COS0H  COS0T,  ♦ SIN0h  SIN0t) 


♦ e*  (-siN6H  co.seT  + coseH  sinqt) 


+ Hj.  g rT  SIN6t  + Mn  g rN  SIN0T 


A - 7 


Writing  Equation  26  in  terns  of  6^,  we  have: 


(29) 


9h  = 


0H 


- \ 


COS0 


vs 


» •!»  r, 


‘H 


(C0S6U  COS0,,,  + SIN0„  SIN0  ) 0 


H 


H 


+ (SIN0h  COS0t  - COS0H  SIN0t)  02  1 + SIN0h 

rH 


where: 


F0H  ■ FH  rH  + T ““ 

F0T  = FH  C0S(eH  - V rN  + Fx  rF  - T 
FXL  = FH  COS0H  * FI  COS0T  + FL  • 
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APPENDIX  B 


Derivation  of  the  Airbag  Algorithm 
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CALCULATION  SHEET 

Page  ‘ ..  Or 
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Calculation  Serial  No 
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1 t 

13)  X'-  yt-X,  +y.'V  u/smEM 
Cosed 

*4jtj'  = y7  - y,  1-  y" -X'S\U^d 

CO5  0J 
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'THE  ReSUlTMT  HoT\oO  OF  “THE.  'STEE.^VJM  ReTAT TOTrIt  CHH-TT 

IS  6)>w£  kL  Qy  ; 
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^cm  or.r,  K] 

|4jowtrJ^  TAqj  £ UAox  '/Je:  H^y  cav-culatx  hVoi  , TV^e  \ju^s  o>*^ 

Me  TTntt  tl<t  e pT  j 


^)  riw  - ^ H 


A.  O 1 DTVl 


UJHE.  f ZL.  Vy  - y|e 

We>  - Ruoy  VM  . 


'Vox  ■ VVH  HApi 

\05)  V^I  - ^voj  + ~^voi  ~ Kaot  VJi\  4 "TAoj. 

VaJh  v c v\  \s  the  e- (pM  .'LXes\pr-E>  So  tha-t  e cp^ . 76  Caa  se  evaluated  . 


Prepared  By 


Checked  By 


(B-23 ) 


Date 


Date 


APPENDIX  C 


Derivation  of  the  Knee  Restraint  Equations 


CALCULATION  SHEET 

Page  - 1 Or 

Project 

Job  Order  No 

Calculation  Serial  No 

JOng  e Ti a 3~gcrrp sy 

Wheaj  the.  hU  Pt  i-l a ^ 
^\/ed  Wp he  than] 

^ ■yp  j £cs  (_  Relative 

1 ’o  TH  El  (Tom  y a, t~t  m g"  m t y 
THE.  roELoWiM  (j  ETpu'S 
ATPlV  , 


AT  T = £- 


H H “ A ! + (w-lj7-')‘'=  Rp  E (('!•  Of  Cltcu^/EHTPO.  JKT  l]  - FT- 

V-Mem  WE  oo'-ve.  THESt  TWO  g cp^J 's  SiiHWCTAMEojwuy  \we  Solve  Foa  the. 
P.Wee  P^var  Po,nT  Pp  THE.  Y-Y  Coofr^MA  TE  Sy5TE M . 

^o<— v/it«s  The  pa  ^ihl ltawEoo^ E_y  yie^D^  ’ 

?j  kxi  + ynlflnSH'  = ^-t-TTvt  w 


THis  EL<$xi.  VJilw  0e  solved  ^UMEpucAui-y  (?y  THf  uew-rcn~$A  so*J  ^femoo  5Pc. 

4)  ij  * M ^Substituting,  \mtp  4-  , we-  obtain  is) 

« / t 

AMouiAf,  y.t  u^  tug  C, ah  Solve  PbfL  ©p.  S 


X- 


- Tan! 


A ^ 


X 


) 


Prepared  By 


Checked  By 


(C  — 1 ) 


Date 


Date 


L 

CALCULATION  SHEET 

Page  2-  Or 

; Project 

Job  Order  No 

Calculation  Serial  No 

, VJ£_  £anL  Pbn.  "TH  £ £©&C^£. 

* 


6)  Ff  - W CD«  0p 


Prepared  By 

(C-2) 

Date 

Checked  By 

Date 

Appendix  D 


PAC  Program  Listing 


1 


PRC 


0 8- 30  •••••£:£ 


1 OOC  ♦ ♦♦  PRC  ♦♦♦ 

1 1 OC 
1 SOC 

1 3 OC  THIS:  PROGERM  PREDICTS  THE  PRSSENGER  K I NEMRT I CS  IN  R CPRSH  S ITURT I ON 

1 4 OC  IN  WHICH  THE  PRSSENGER  IS  PESTER I NED  BY  RN  RIEBHG  RND  KNEE  PESTER  I NT 

1 5 OC  THE  RIPERG  IS  MOUNTED  TO  THE  DRSH  IN  R SPECIFIED  POSITION  RND 

1 6 OC  DEPLOYS  TOURED  THE  PRSSENGER  RT  R DEPLOYMENT  RNGL.E  SPECIFIED 

1 7 OC  BY  THE  USEE.  THE  COMPLETE  DELOYMENT  PROCESS  OF  THE  BRG  IS  MODELED. 

1 soc 

1 9 OC  THE  KNEE  PESTER I NT  RBSDPBS  THE  KINETIC  ENERGY  OF  THE  LOWER  BODY. 


£ 1 0 C THE  PR S S E N G E P I S M 0 D E L F D B Y FOU P M R S S E S -R  HERD  M R S S - TWO  CHE S T 

28 OC  MRSSE'S  RND  R LOWER  BODY  MESS.  THE  PRSSENGER  IS  CONSTPRINED  TO  HRVE 

23  OC  F'LRNRR  MOTION  SO  THRT  THE  PRO  GERM  IS  STRICTLY  EPPLICRBLE  ONLY 

£4 OC  TO  FRONT HL  CPRSH  S ITURT IONS. 

80  OC 

26  OC  PESTER  I NT  PEEFORMRNCE  CRN  BE  OPTIMIZED  BY  RPF'POPR  I RT  E C HR  NGE  S' 

27 OC  IN  THE  DESIGN  PRRRMETEPS . 

28  OC 

28 OC  TYPICRL.  DESIGN  PRRRMETEPS  THRT  CRN  BE  E VRL UR TED  RPE  BRG  SIZE,  BRG 

3 OOC  SHRF'E-  I NFL.RT  I ON  CHRPRCTER I ST  I CS  ■>  VENT  RPERi  INFLRTOP  LOCRTION  IN 

3 3 0 C BOTH  THE  U P— DO W N R N D FOPE-RFT  DIRE C T 1 0 N S * B R G D E P L 0 Y M ENT  R N G L E - 

32 OC  RIPBRG  FRBRIC  HEIGHT,  LEG  RNGLEG  CHEST  MRSS  RND  DRMPING  PROPERTIES' 

3 3 OC  RS  WELL  RS  OTHER  SYSTEM  PRRRMETEPS „ 

34  OC 

35 OC  THIS  PROGRRM  IS  SELF  CONTRINED  IN  THRT  NO  EXTEPNRL.  FUNCTIONS  DR 

36 OC  SUBROUTINES  RPE  REQUIRED. 

37  OC 


38  OC 

3 8 OC 

4 0 OC 

41  OC 

42  0 C 

43  OC 

44  OC 

45  OC 

46  0 
4 7 0 
48  0 

4 9 0 
5>  i 

5 0 0 
51  0 

58  0 

53  0 

54  0 
550 
56  0 
570 
580 

59  0 


RUTHOR:  MICHREL  F I TZF'RTR  I CL 

P I TZF'RTR  I CL  ENG  I NEEP I NG 
WhF'SRU  ' I ND I hNh  4658 : t 
TEL.  1 8 1 9 1 - 8 6 7 — 4 4 7 

ORIGINALLY  WRITTEN  NOV.  83'  1981.  THIS  IS  REVISION  RS  RUG.  £8-  1 988 


FILENRME  INFILE 
C 0 M M 0 N E L R G - N 0 U T P IJ  T 
PER  L N E C K C H L' ' K R N ' L.  T - L F 

COMMON-- OUT,--  NF  D ? T *.  1 r-  5 1 ' u 1 6 ' 1 r 5,1  ' 1 1 6 ’ 1 i-  5 1 , ,;-:!8  ■:  6- 1>  175.'  ' X..;;  >:  6 1 1 75."1  ' ■ r 1 6 , i 

KY8 '6'  1 75'1  !»  Y9  1 6'  1 75 '•  •<  Y1  0 1 6'  1 75’>  ’ XI 1 C6 » 1 75-1 

COMMON  - OUT  1 ••  NOUT ' TOUT  1 75">  ' X4  C6  '175''  X5 '6 , 1 75  > ' X6  -:'6  - 1 75  > 

C 0 M M 0 N TIME S TER,  X STOP,  DELTRT n R I N T 1 ■>  P I N T 3 
COMMON 8H I C.-  TH  I C 175  --'  HRGS  C 1 75"-  ' CPGS  < 1 75> 

COMMON, -NR  ME-"  I NF I L E 
1 05  0 PRINT  "INPUT  PILE  NR  ME " 

INPUT  'INFILE 

PRINT' "ENTER  1 IF  YOU  WRNT  FULL  LIST  OF  OUTPUTS  ENTER  2 IP  YOU 
Z-.  WRNT  RBBREVIRTED  LIST." 

INPUT, NOUTPUT 


D 


1 


PRC 


6 0 0 

NPD  = 0 

fc.  1 0 

NOUT=0 

68  0 

CRLL  SOLVE' 8"' 

b 3 0 

IF  t'NPD.  GT.  175  ' 

NPD= 1 75 

6 4 0 

IF  - NOUT. RT.  175 

) NOUT=l  75 

65  0 

1 1 8 0 

FOPMRT  • 1 H-'i 

66  0 

1 1 85 

FOPMRT - F7. 4 , £F6. 1 , 8F7. 8 

67  0 

1 1 3 0 

FOPMRT  •;  1 ',■■■  i 7F1  1 

;i 

6 8 0 

114  0 

FOPMRT CV) 

6. 9 0 

1 1 5 0 

p n p m h t 1 1 x i r * *4-  x 

, ■ 

7 0 0 

IF  • NO!  ITPI  IT.  FC. 

£) RO  TO  1 

7 1 0 

PRINT  H80 

780 

1 1 7 0 

FOPMRT  ."IV,  ■’ 

T I ME 

73  0 

74  0 

750 

7 6 0 
7'  7 o 
78  0 
790 

8 fi  0 

8 1 o 
88  0 
830 

84  0 
8 5 0 
86  0 
87  0 
8 8 0 
8 9 0 
9 1 j o 
9 1 0 


VFH  VEL 
- MO- 


VE N DIOR  ONES 

( H " S -'  '''FRS'« 


L PS  ) 


VFH  6 •••  ' 

: ■ T p P CMR  FDPCF  CPP  F 0 P C E " 1 , 1 

f;,  <;  I N C H E S'  > >'  I N C H F S "'  ( L E S ) 

PRINT  1170 
PRINT  1150 
DD  1881  0=1 , NRD 

1 88 1 PP I NT  1 1 8 0 , T <"0  '•  - <>;  0 «:  J , 0 ) , J=  1 , 6 
PRINT  1180 
1883  FORME t cik, - 

f:  RCF  FEM  RNR 
8 < L PS) 

PRINT  1883 
PRINT  H50 
DO  1830  0=1 i NRD 

1 63  0 PP  I NT  11  3 0 i T 1 0 1 i 'VI  1 .1  , 0 1 , J=  1 ■»  6 •' 

PRINT  1180 

185  0 FOPMRT  1 1 '/ 1 " TIMF  TORSO  DISP  TORSO  RNR  TORSO  VEL  TORSO 

:•  RCC  TORSO  R.D.  TORSO  R.  V.  “.-"ISN 


TIME  H-R  R.D. 

TIE;  R N R ••  1 , •• 

CL  PS)  'DEG'-' 


H-P  VEL  SERT  FR.  FEM  FO 
1 M s 1 1 I NI...  HE  S ' 1 FF  s : < 

1 D E 6 1 ' 1 1 


• MS;  :• 


8 > DER '■ 

PRINT  1850 


1 D -"SEC-' 


D -''S  EC  ♦♦8>  •:  INCHES) 


INCHES  :* 
(FPS)  “ ) 


98  0 

PRINT  1150 

9 3 0 

DO  1310  0=1 , NRD 

94  0 

1 31 0 PRINT  1 1 3 0 , T 1 0 ) , • 

■ ! ci  *'  .J  * f-  1 •» 

J=t , 6) 

95  0 

PRINT  1180 

96  0 

133  0 FnpMRT  >:  is;,  ” 

T I ME 

HERD  DISP 

HERD  RNR 

HERD  VEL  HERD 

97  0 

C RCC  HERD  R.D. 

HERD  F 

. RNR ” ■ 1 V , ' 

(MS) 

( INCHES) 

98  0 

S •:  DER) 

D SEC ) 

CD  ' S E C++8) 

•INCHES) 

1 DER)  " 

9 9 0 

PRINT  1330 

1 0 0 0 

PRINT  1150 

1 0 1 0 

DO  1380  0=1 , NRD 

1 08  0 

138  0 PRINT  1 1 3 0 , T . 0 ) , 

i ; S ( ! , (■•;  .. 

, ,i=  l , 6) 

1 0 3 0 

PRINT  1180 

1 i j 4 1 j 

1400  FOPMRT  IV,  " 

T I ME 

E'  ERR  RCC 

REV  UP  RND 

REV  MR  CST  E F V 

1 05  0 

;•  1,1  R DSN  RED  MR 

end  red 

mr  d s h " i ;• 

" CMS) 

G y 3 1 

1 06  0 

8 (FPS) 

(FPS) 

CFPS) 

■INCHES) 

• INCHES'  " 

1 07  0 

PRINT  1400 

1 08  0 

PRINT  1150 

1 08  0 

DO  1460  0=1, NOUT 

D - 2 


PRC 

08-'' 3 0 "88 

1 1 0 0 

1 4 ft  0 

PRINT  1 130- TOUT ck 

1110 

PRINT  1120 

1 1 £ 0 

1 4ft  5 

FGPMRT • 1 S, " 

1 1 3 ft 

c!  ! i IP  D l H Li F D i u It* 

114  0 

S:  CPF'S) 

1150 

PRINT  14ft 5 

1 1 ft  0 

PRINT  1150 

117  0 

DO  1 4ft 8 f"  = l - NDHT 

1 1 8 0 

14ft  8 

PRINT  1130- TOUT '4: 

| 1 Q ! j 

PRINT  1120 

\ £ H ! j 

1471 

FGPMRT  Cl S-  " 

1 £ 1 0 

fi  MR  STN  STD  UP 

1 ££0 

t CLBS) 

1 £ 3 0 

PRINT  1471 

1 =24  0 

PRINT  1150 

1 25  0 

DO  1474  F = i - NOI.iT 

1 2 ft  0 

1 474 

PRINT  1 1 30- TOUT (k 

1 2 7 0 

PRINT  1120 

1 2 8 0 

1 4 8 0 

FGPMRT ' 1 S- " 

1 O Q i'j 

J.  U J 

?'■  FORCE  HP  FORCE 

i 2 o n 

S CPS'  I E)  C 

1 3 1 0 

PRINT  1480 

1 32  0 

PRINT  1150 

1 33  0 

DO  154  0 K = 1 - N P D 

134  0 

1 54  0 

PRINT  1 1 30  - T CL)  - >: 

1 35  0 

PRINT  1120 

1 36  0 

1 5ft  0 

FORMRT Cl  X - " 

1370 

C HERD  SI". 'IX- 

1380 

C.;  CE"'  S)  CE  S' 

1 39  0 

!l:  CE  " S' ) " ) 

1 4 0 o 

PRINT  1 5 ft.  0 

1 4 1 0 

PRINT  1150 

1 42  0 

DO  1 ft. 2 0 K= 1 - NPD 

1 4 3 0 

1 ft80 

PRINT  1 130-  T Cft)  , 

144  0 

PRINT  1120 

1 45  0 

1 6£5 

FORMRT CIS- " 

1 4 6 0 

L VC  B CTP 

1 47  0 

L CHE  S'"1  CINCHES 

148  0 

PRINT  1ft 25 

1 4 9 0 

PRINT  1150 

1 5 0 0 

DO  1ft £8  K= 1 - NPD 

1 5 1 0 

1 ft28 

PRINT  1 130-  T CL)  , .:' 

1520 

GO  TO  1ft 37 

1 5 3 0 

1ft-  £9 

PRINT  1120 

154  0 

1ft  30 

FORMRT c l X - " 

1 55  0 

& FEM  FORCE  PB 

15ft  0 

%•  ( It  :» 

1570 

PRINT  1630 

158  0 

PRINT  1150 

1 59  0 

DO  1ft. 3£  L = 1 - NPD 

TIME  U BRE  mil:  IJBV  MR  END  UBV  WF  C ST  IJBV 

N D IJB'B  l.i P D S N " ' l y , M CM 3 > ■:  E •"  g 1 

C F P S'  ) C F p S > C I N CHE  S C I N CHE  S ) ' ' ) 


TIME  CST  F BSP  STN  F BSP  ST'v'  MR  LSI  RL.B 


CEPS  ) 


CINCHES) 


CMS:1 

c.i  nt  hes  ;• 


CL  ESC) 

I NCHES  ) " ' 


T I ME 


HERD  BP, 

■ n l " is-  1 

CL  BS  > 


B R E V DL.  B R E PRESS' . HI.'.-- ft 

C M S ) 1 I N f H E 3 ) * C l. 

ecu. IN. > " 1 


TIME 


CHEST  ftp  CHEST  SI  HERD  ftp 

MS)  Cft  •"  S) 


TIME  S' T N RCC  POLL  PRD  SC  B CTP 

l.ilRB  I.IUPB"  •-'IX'.  " CMS'1 

INCHES'-  CLBS  ) CLBS)";. 


I T ' I 


I I 


TIME  CST  R ES  HD  R ES 

V I..I R c S T B ft E P R E S S ‘ ‘ 1 X - 1 ’ 

CE S)  CLBS)  ' EPS:;. 


STN  RCC 

MS)  CG"'S> 

CPS  IS) 


i.  IN.  1 


C I N 


D - 3 


4 


PRC 

08  •■■'30 -■•'88 

1 6 1 j i j 

1 6 8 

PRINT  1 130?  T CK j - 

C Y 1 Ij  i'  J , L 

j , ..J  =1,6  j 

1 6 1 0 

PRINT  1180 

1 68  0 

1634 

POPUP T Cl  7, " 

TIME 

TDP  !:•’ D H N H 

HERD  PNG 

FEMUR  PNG 

1 6 3 0 

TP  CST  H-P  P, D» 

CHES  T 

DEF'L 17,  M 

CMS  j 

CDEG) 

1 64 1'* 

T C fiEG  'i 

CEPS:' 

CINCHES':- 

CINCHES) 

" '*1 

1 65  0 

PRINT  1634 

1 8 6 1 j 

PRINT  1150 

1 6 7 1 j 

DO  1636  K= 1 ? NPD 

168  0 

1 636 

PRINT  1 1 3 0 ? T 1 K ''  ? 

* !**•  i i 1 . J •»  K 

1 * = 1 ^ A 

1 b 9 l’1 
1 7 0 0 
1 7 1 0 
1 780 


1637  PRINT  1 IPO 

PRINT, "ENTER  1 TD  CELT UL PTE  HI 
INPUT  , NE'E  'I 

I F N P E T ,,  N E . 1 •<  G 0 T 0 P 0 0 0 


1 730 
1 74  0 
175  0 
1 76  0 
1 77  n 
1 7 3 0 
1 79  0 
1 3 0 0 


164  0 


PERE=  ij  „ 

NT TDP = NED 
DO  1715  1 = 1?  NSTOF' 
DU  1716  J= 1 - I 
L = I + 1 
SUM=0. 

DO  1717  K = 1 ? J 
L=L- 1 


1 8 1 0 

1 8 3 0 
1 8 4 0 
185  0 
1 86  0 
1 87  0 
1 88  0 
139  0 
1 9 0 0 
1 9 1.  0 
1 98  0 
1 9?  0 
194  0 
1950 
1 9 6 0 


SUM==SUM+HPGS  c L > ♦PI NT 8 

1717  CONTINUE 
DEL  T = T H I C c k ) 

CHECK  =o  UM--DEL.T 

I F PERK-CHECK 1 7 1 8 ? 1 7 1 6 ? 1 7 1 6 

1718  PERK= CHECK 
TLni.!=  ■ L — 1 '«  ♦PINTS 
TH'IGH=  I ♦PINTS’ 

1716  CONTINUE 
1715  CONTINUE 

HTC=PEHL  ♦♦£■  „ 5 
PRINT, "THE  NIC  I S " ? H I C 
PRINT,  " T 1 = " , TL.OW 
PRINT? " T£= " , THIGH 
8000  8 TOP 
END 


1 98  OC 
1 9 9 0 C 
P 0 Li  0 
8 Li  1 ij 
8 0 8 0 
8 1 j 8 0 
8 04  ri 


THIS  SUBROUTINE  SETS  Ijp  THE  DIFFERENT  I PL  ERUPTIONS  THPT  DEI 
THE  PH S S E N G E P L I N E MPT I C S . 

S U B P D U T I N E DIF  E Q < T , V ? D j ' ) 

COMMON  ■ MPNDPT.--  ZL  -779  ZS  ? ZH  ? PT  ? PN  ? PH  • PTOPH  - ■■••87  ? Y8Z  ? MB  - B 0 - MH , t, 
DOUBLE  PRECISION  Y C8) 

DIMENSION  DY  C8> 

C P L L F 0 P i"  E T H C Y ? T N E C K > 


3 05  0 
8 1 j 6 1 j 
8 07  0 
8 1 j 8 0 
8 0 9 0 


CPI.  L DEC  EL  - T,  GC  j 

C H L L BP  G :S  U F < T ? Y , [i  Y - T N E CL  ? FTH?  E K - E T T , G S > 
SH=S  I N -:  Y C6)  > 

ST=S  I N CY  •:  7 • > 

C H=COS  CY  ■'  6 T ) 


TV  W 
-"DEG:- 


:c  FIFE 
. T ? L E 


D - 4 


PRC 


1 1 i 


£ 1 0 0 
£ 1 1 i J 
• £1£0 
£ 1 3 0 
£ 1 4 n 
• £150 
£ 1 b 0 
£ 1 7 0 

£ 1 '4  0 
£ £ 0 0 

££  1 0 
Z:  Z‘  c'i 

U..  I—  L_  V 

££3  0 
££4  0 
££5  0 

££  b or 

££7  0C 
£ £3  OC 
££9  0 
£3  0 0 
£ 3 1 0 


£34  0 
£35  0 
£3b  0 
£37  0 
£ 3 3 0 
£ 3 9 0 
£4  00 
£410 
£4£  0 
£430 
£44  0 
£45  0 
£4b  0 
£47  0 
£43  0 

£ 4 '4  0 

• £5  0 0 
£51  0 
£ 5 £ 0 

. £53  0 
£54  0 
£55  0 
£5b  0 
£57  0 
£58  0 
£59  0 


C T = C □ :S  V ' ' 7 1 '■< 

DY  1 1 1 = 1 FX—  ■:  ZT+FT +ZH-*F‘N : ♦ 1 C T ♦ PY  >'3  —ST ♦V  ♦V  ' 3 ' •' 

k. ■ — 7 H ♦ P H ♦ 1 L H + DY  1 £ 1 - i H+  V '■  £ •’  ♦ Y 1 £ ' 1 1 / 1 ZL  + 4.T  + d. H 1 

DV  ■:  £ :<  = - FTH-ZH^PH^DV  ’ l '1  ♦CH-ZH^PN^PH^  . C 7>CH+DY  <3>  +C7>S'HVy  - • 

ft  ♦'£  i : 3 •'  —ST ♦ CH+V  ■:  :£>  ♦V  > 3 ’*  +ST ♦SH^DY  '■  '3  ’• 

;>  CZH+PH^PH':1 

DV  ''£;>  — y'FTT-  1 ZT+PT+ZH+PN.’  ♦py  Cl-1  ♦CT-ZH+PN^FH*  < CT^CH^PY  • £'*  +Y  >:'£Y 
& ♦ Y £ '■>  ♦ *•'  S T ♦ C H— C'T  ♦ R y + S T ♦ S H ♦ 0 Y 
b . ( z T ♦ P ' T ♦ p T + £ H ♦ P N ♦ P N -i 
DY  1 4 ■'  =-9S 
PY  1 5 1 =Y  ''  1 > 

DY  <b>  =Y  '•'£'■• 

DY  »:'?>  =Y  >'  £'' 

DY  v'S'y  = Y 1 4 '' 

PETUPN 

END 


THIS  SUBROUTINE  PER DC  IN  THF  INPUT  DRTR,  SETS  UP  THE  INPUT  DRTR 
FOP  D I S P L R V R N D I N J T I R L I Z E S 0 E Y V R P I R P L E S . 

S I J p P 0 U T I N E S E T I.  J P X - V 1 
FILENRME  INFILE, OUTE ILF 
p E R L N E C k C H k , k P N 1 L T 1 L F 

C □ M M 0 N S E R T F P I C •••- N S F - S U N , S FUT  - PEL S F , S F N < £ - £ 4 
c o M M o N S T F □ P C E - N ST,  S T E C £ - £ 4 
COMMON  CFO p C E NC - CF < £ - £4 • 

C 0 M M 0 N k N F EPF  S T - N k p , 3 k P , P U T , PE L k P , k P N ' £ - £ 4 

COMMON  •-NRMF  I NF I L E - OUTF I L E 

COMMON  -T I ME  • STEP  , XS  TOP"  DFLT'RT , PINT  1 •*  F'  I N T £ 

COMMON  • MRNDRT  ■ ZL? ZT, ZS- ZH, PT, F'N, PH, PTOFH- Y£Z- Y£Z- UP, BO, UH, L 7',  LF 
COMMON --MMRNpRT  STDRMP,  C DRMP , CHESTT,  HFRDT,  THF 0 , THLO,  DELTTP 
COMMON  NEC'k  NPN,  FNECk  < £,  £4  • , DC  N 
C 0 M M 0 N V EH  N V , V EH  OS  •:  £ , 5 0 '• 

C 0 M M o N ■ R R s F L 0 • N p R , R F N <"  3 , £ 4 '■ 

C 0 M M o N c R DRT H T M o p , P R Z , R T Z - I 1 - P N 1 , p N £ , P N 3 

COMMON.- BRRDRT  •-  VC  1 , V C £ , R V , FSR-  FSC,  FSB,  RO,  UFPR,  DPOLLZj  THETRp,  XI  , Y1 
COMMON -MBRRDRT  DM' , DINE, USOCk , DUPBZ, VCRH, VCRV 
c , R , S B , 3 C , 3 R P 

c 0 M M 0 N M D R T ■ D l J P BRH , DU P B R V - D P B R H,  DP B R V , VIJ P PRH,  V P B R H - V U P P R '• , V P B R V , 
S RIJF'BH,  RI.JPBV,  RP PH , RPBV , VSRH 

COMMON  - prppm  pft , THETRTZ,  THF'TRHZ,  PR  1 , FKNEE?  RP P,  DPP,  ROI 
C OMMON  - MF'RPRM  • S F , F " TBS  - PR  BS , Rk iPB , DUPP , DPC , DPS  - DS  C 
C OMMON  •-'MM  T SC  FM£ , PP8 , PR5 , RT , FPN , VOL  0 , Rl.l 
DOUBLE  P PEC  I: SI  ON  Y • 3 ;• 

X=  0 . 

£07  0 FOF'MRT  CV> 

£ 08  0 F OF'MRT  1 Y , "INITI RL  VELOC  I TY : " , R 1 0 . 3 - 1 X , 

C " I N I T I RL  HERD  R N R l E : " , R 1 0 . 3 •- 1 X , "TNITI  RL  TOPSO  RNRL  E : 

CR10.3> 

£100  FOF'MRT  '1  Y, 

k " ML. EG  MTOF'SO  MSTEPN  MHERD  PT  PN 


D - 5 


p H C 


HI-!..-' :-!!"!  ■ 


£ £ I'!  ij 

RH 

RTOF'H 

" 1 X ’ 8b  1 0. 31 

£61  0 

£ 1 £ 0 

PDF MPT 

'-  1 X " 

at  MOP 

PGZ 

GTZ  U 

RN1 

£ 6 £ 0 

RNc 

PM:?" 

1 i 7 A 1 IJ  „ . i 

£63  0 

£ 1 3 0 

FT!  PM  AT 

> 1 ; , " 

VC:  1 

VC£ 

RV  F3h 

FSB 

£64  0 

fj. 

FSC 

X 1 

V 1 " 1 X , 

8G 1 0 . 3! 

£ 6 5 0 

£ 1 38 

FORMAT 

1 1 X i 

A 0 

THETAD 

FABWGT  ST DAME 

C DAMP 

£ 1 6 0 

&: 

DM3 

DINF 

W30C  L " 1 X , 8 

G1 0. 31 

£ 6 £ 0 

£134 

FORMAT 

'■  1 X •> 

1 1 j M 

DR OL LZ 

V p “7  Y p "t* 

£68  0 

ft.. 

WE 

LF 

DON  ".-TV, 

8G 1 0 . 31 

p 6 9 n 

£137 

FORMAT 

'■  i x , 

A A r FLOW  T I MR  •-  1 X . 1 00 1 o . 9 

i 

P 7 fl  fl 

£1  36 

F OF'MPT 

r 1 X , " 

GAS'  F'LOli*  - 

LB-TEE  ".  IX,  1 OG 

1 0.  31 

£71  0 

£ 1 4 0 

FORMAT' 

' 1 X h ■' 

MR CL  ANGLE " 

' 1 X , 1 OG  1 0 =.  3 j 

£7£  0 

£ 1 5 0 

c nt mmt 

i | V » 

ME  r t TO POME 

- FT -LBS"  TO, 

1 OG  1 0.  31 

"7  f“j 

£ 1 6 0 

FORMAT 

' t , 

MR T£  NEC'L  NETS  kp 

NETS'  VEH  NETS  SEAT 

NRT'S  G 

£74  0 

6 El  . 5 

T 

SL  .IF-  . l x - 7 

A 1 Ij  . T 

£750 

£ 1 7 0 

FORMAT 

'•  l x - " 

MPT 3 3 

ter  MRT3  C'HS'T " 

• 1 X , 

£ 7'  6 0 

p£G i n. 

-1  ! 

*7  *T*  l”j 

p £5  m 

FORMAT 

* t X i 

VE'H.  PULSE 

- TIME  ".-  IX,  1 OG 

1 0 . 31 

£78  0 

££6  0 

FORMAT 

• i x:  •• 

VEH „ PULSE 

- DECELERATION 

“ -•  l X , 1 OG  1 0 » 3.1 

£790 

3 £6  7 

format 

'•  i x i 

SEAT  FP  T C T I ON  PI  S'  R L A L E M E NT"  1 X * 1 0 G 1 0 . 3 ■' 

£8  0 0 

£ £ 6 7 

F OF  MAT- 

'•  i ■ ’ 

SEAT  FPICTIOM  FORCE  - LBS 

" 1 X , 1 Obi  0. 31 

£81  0 

£ £ 6 8 

FORM  AT 

''  1 X - " 

L M EE  D I S' P l.  A C E M E M T " - 1 X - 1 0 G 

1 0 . 3 1 

£8£0 

P £ 6 9 

format 

•"  .1  X - 

KNEE  force; 

- LBS"  1 X, 1 Obi  0. 31 

£ 8 £ 0 

££70 

format 

''  1 X " " 

STERNUM  DTS 

PI.  AO  EM  ENT". TX, 

1 OG 1 0 . Si 

P R 4 i"! 

££71 

FORMAT 

i | x » 

STERNUM  FORCE  - LBS  " ■"  1 X - 1 OG 1 0 . 81 

885  n 

l Z C..  •'  c 

FORMAT 

1 1 X " 

CHEST  D I S' RI- 

ACEMEMT" IX,  1 OG  1 0.  3"' 

£ 9 6 0 

;!•  p p -Ij 

format 

' 1 X • " 

CHES' 7 FORCE 

- LBS " 1 X , 1 OG 

1 0 . 31 

£87  0 

c c!  «■''  4 

format 

■ 1 X , 

THFO 

THLO 

XS  TOR  STEP 

PINT  1 

P I N 


Tc’  " •••■’  1 X ' Sb  1 U . 1 


£88  0 

READ  ■ 

INFILE!. 

lSL 

0 7 1 j 

1 V 1 41  , V 1 61  , 

Y >''71 

£89  0 

READ  < 

I ME  I LE , 

u. 

0 7 0 

1 ZL  ’ 7 T , 7 S' , 7 

A,  FT,  F'N , RH,  RTDF'H 

p A fl  fl 

RE  AD  ■ 

I NEIL E, 

07  0 

1 NE'N  , NL  F' , N'v 

, NSF, NRG, SUN, SLR 

£91  0 

READ  >' 

INFILE, 

CZ. 

0 7 ij 

•1  NS  T,  ME 

£ 9 £ 0 

READ  <' 

infile. 

C 

0 7 0 

■'  < G E N 1 1 , L 1 , 

k -1 , NRG! 

P 9 3 fi 

READ  < 

INFILE, 

0 7 1 j 

1 c GEN  L‘£  , k 1 , 

k = 1 , NRG"' 

£ 9 4 0 

READ  • 

INFILE, 

i7. 

07  0 

1 ATMOR, RG7, 

GTZ,  IJ,  F'N  1 , RN£,  RN3 

£95  0 

READ  •: 

INFTLE, 

CZ. 

07  0 

1 VC  1 - VC  S' , A'v 

•»  F" ' h * F N'B 11  F 3 L 1 X 1 ^ Y 1 

£96  0 

READ  1 

INFILE, 

CZ. 

1 j 7 1 j 

*>  h U ^ THF.ThD  ? 

FABWGT , ST DAMP , C DAMP , DM3 

£97  0 

READ  < 

I NF I LE , 

CZ 

070 

l UM,  DPnLL.7 , 

Xc‘7 , Yc'7 , MB , L F , BLN 

£99  0 

READ  •' 

INFILE, 

l! 

07  0 

1 THRO, THLD, 

XS TDP , S TER , P I NT  1 , P I NT£ 

£99  0 

READ  < 

INFILE, 

£ 

07  0 

1 lSFN  1 , K 1 , 

k = 1 , NSF! 

:';!  0 0 0 

READ  ( 

INFILE, 

CZ 

07  0 

'■  ( S'  F N 1 £ , K 1 , 

k = 1 , NSF) 

8 0 1 0 

READ  <' 

INFILE, 

CZ. 

1 j 7 ij 

1 FNECK  < 1 , K 

.*  «i  fc  = t 1 N F'f ‘-j  * 

3 0 £ 0 

READ  '• 

INFILE, 

CZ 

07  0 

1 CFNECK  -:£,  k 

1 , k'=  1 , NPN"' 

3 0 3 0 

READ  •: 

I NF I LE , 

CZ 

07  0 

1 *■  V E H ^.7  £ 1 !»  f: 

1 , k = 1 , N 

3 04  0 

READ 

I NF I L E , 

CZ 

07  0 

j L V E H b S’  1 £ , L 

1 , k = 1 , NV1 

3 ij  5 0 

READ  t 

INFILE, 

CZ 

07  0 

1 . H F*  N 1 1 k • 1 

K=1 , NL  R 1 

3 06  0 

READ 

INFILE, 

CZ 

0 7 ij 

.'  '.'  k R N ' £ , k 1 , 

k = 1 , NLP! 

3 07  0 

READ  1 

INFILE, 

CZ. 

07  0 

1 • S TF  ('  1 , L I , 

L = 1 , NST1 

3 08  0 

READ  -' 

INFILE , 

CZ 

07  0 

l >:  S TF  '£,  k l , 

L = 1 , NST! 

3 0 A Ij 

READ 

INFILE, 

CZ. 

07  0 

J a ft,  k 1 , k 

=1 , NCI 

DINF • M5DCK 


0-6 


1 0 0 


READ ' INF! L F - £ 07 0 - 0CF C£ - k > ,K=1j NC > 


110  PRINT  £490 


1 £ 0 

PRINT 

" I NF 

IJT  VALUES 

— INPUT  UNTTSO  MSEC,  MEN , DEGREES , 

1 3 0 

8 INC HF 3-  L E8 " FT-LFC 

, 6 •"  0 > ■' 

1 4 0 

£48  0 

FORMAT 

- 1 O'  • 1 

061  0.  3'1 

150 

£49  0 

FORMAT 

• 1H-"- 

1 6 0 

£5£  0 

PRINT 

£ 08  0 1 

V >'  4 1 , V ' 6 ■>  , 

V < 7 1 

170 

PRINT 

£ 1 0 o , 

ZL ■ ZT ’ ZS - Z 

H'  FT,  RN i PH i PTOF'H 

1 8 0 

PRINT 

£ 1 6 0 ' 

NF’Ni  NO  R,  NV 

, NOP i NPP, SUN, Sk P 

1 9 0 

PRINT 

£ 17  0, 

NO T , NC 

£ 0 0 

PRINT 

£1  35, 

1 PEN  1 1 , k 'i  , 

k=l , NPPJ 

£ 1 0 

PRINT 

£ 1 3 8 , 

i PEN Op, k j , 

k = 1 , NR P :- 

£ £ 0 

PRINT 

id  d b 5 ■» 

1 z F N 'li  k 1 i 

k =1  < NSP:- 

£3  op 

PRINT 

q r A {•'  * 

’ i F H 1 d *»  '*  i 

k =1 , NSF> 

£ 4 0 

PRINT 

£ 14  0, 

OFNECk  Cl  , k 

:.  i , k = l - NRN'j 

£5  0 

PRINT 

£ 15  0, 

op  NECK  >:£,  k 

> , k = 1 , NF  N -, 

£6  0 

PRINT 

£85  o - 

OVE'HPS ' 1 , k 

;•<  , k = i , nv> 

£7  0 

PRINT 

i d.  d b U * 

OVEHPO  op,  k 

, k — 1 , N V 1 

£80 

PRINT 

£ £ 6 8 ■* 

o k F'N  1 1 , k > , 

k = 1 , Nk  P> 

£90 

PRINT 

££69  ’ 

< k R N '•  c ' k i , 

k = t , Nk  R > 

£ 0 0 

PRINT 

££70, 

os TF  '■  1 ■.  K>  , 

k:  = 1 - NC  T> 

3 1 0 

PRINT 

£ £71  ’ 

1 0 T F , 

k = 1 , NC  T> 

3£  0 

PRINT 

c c.  i' d * 

■ C F o l , k > , k 

= 1 - NC  :■ 

33  0 

PRINT 

d d 7 d •» 

1 ! _ F 1 d *•  1 i k 

= 1 , NCI 

34  0 

PRINT 

£ 1 £ 0 - 

ATMOF' * F’PZ ' 

P T Z , IJ,  F'N  1 , F'N£  , F'N 3 

35  0 

PRINT 

£ 1 3 0 - 

VC  1 , VC  £ , A‘y 

i F C A ' PSP , F'SC  , X 1 , Y1 

36  0 

PRINT 

d 1 :•  id.  i 

PO, THFTAD- 

RApl.IRT  - r T'DAMP , CRAMP , RMS  - D I NF , 1.10  DC  k 

3 7 0 

PRINT 

£ 1 3 4 , 

I.JH , DROL.LZ ' 

!:<£Z ? V £ Z ’ l.iP , L F , DUN 

3 p !"l 

PRINT 

d id  7 4 1 

THRU! , THL0 1 

XO TOP < STEP i PINT 1 , FT  NT  £ 

3 9 H 

V o £ •>  = 0 . 

4 I"!  0 

V C3>  =0. 

4 1 0 

V • 4 ■'  = Y ' 4 1 + 1 a 4F  b b F b 7 

4 £ 0 

V i c,  •.  = n m 

4 3 0 

Y o f,  = y i p > ♦ . 0 1 7 4 5 3 £ 9 

44  0 

THETAHZ=Y C6> 

450 

V C 7 ;i  = Y 1 7 J ♦ . ij  1 7 4 5 3 £ 9 

460 

THETATZ=Y 07> 

470 

Y 08' =0. 

48  o 

Y 0 1 :■  =Y  04  > 

4 9 0 

THPn=THPn^, 01745389 

5 0 o 

THLD=THLO+. 01745388 

5 1 0 

LT  = o Y£Z+LF^S I N - THFO  ■>  1 ■•-'£ I N 0 THl.  0> 

580 

DELTTB= 0 X £ Z - L T ♦ C 0 S 0 THLO>- L F ♦ C □ 0: OTHFO > > ♦ C 

530 

c~  L — 4 L £ u £ 

54  0 

ZT=ZT.--3£.  £ 

son 

77  “ — "7  \ . *■  p p 

56  0 

ZH=ZHP3£. £ 

57  o 

PT=PT-T£. 

58  0 

PN=PN.-T£. 

58  0 

F’H=F'H.y  1 £ . 

D - 7 


PRC 


HR  0 • 


3 S 0 0 
361  0 


364  r 


6 6 '' 


36  7 ( 


6 3 !' 


■3  6 R 


PTdPH=PTC!PH • 1 3 . 

SUM=SUM*1£. 

3:k  P=Sk  PM  3. 

THETRD=THETRD* . 01  746839 
1 = \ + , f ( f-i 3 - p t n p g . > + f; [I:”  'THFT RD 
V 1 =Y  1 + < DM  3-D  I ME  - 3 . > ♦S I N ''THETRD-1 
PPo=  - 3 . CRM  1 + 1 . .-PHI  , RN 1-1 . 

F M £ = V C 3 ♦ 3 Q P T 'PPP++  , PN1  > -PPR  + + . .-pNl  +1  . > - PNl  > 

P R = p G 7 + R T M □ P 

PR6=PR 


1 ! M 


331  f 


GT=GTZ 

FPM=PN3 


333  0 


334  r 


B 0=WB 
CHE3 TT=B 0 
HERDT=  1 . 3+I.IH 


3 3 6 0 3 R P = 6 » 3 8 3 1 8 ♦ ( 3,  ♦R  0 ♦ 3 9 F' T <•  <' F 3 H * ♦ 3 + F 3 C ♦♦£  :>  .-•••  3 . > + F 3 B ♦ F 3: 1"  +FS  R+F 3 B + F 3 R ♦ F 3: i_  :> 

3 3 6 3 l:J BRb  = '':  R P ♦ P R BW G T r‘  1 1 3 6 „ 


D U P B Z =3.  ♦ ' I'  M 3 - 1<  I N f - 


“■  9 f j 


9 8 0 f 


394  f 


DUPE=DUPBZ 

DPE=DIJPBZ 

D U P F R H - - D U P B 3 ♦ C n 3 < T H E T R D > 
DU  p B G V V = D I J P E ? ♦ 3 1 N ••  T H E T R D > 
DPBGH=DUPEGH 

dpegv=dup BGV 

',,'l.iPEGH-V  . 4",  ♦!  £, 

V'UPBGV-  0 . 

VP BGH=Y  f 4 •'  ♦ 13. 

VPEGV=0. 

GUP B= 0 . 

GIJPBH=0. 


3 9 0 C 


GUPEV=0. 


39 1 r 


G F1  E — 0 » 


3930  GRBH=0„ 

3930  gp FV=0. 


394  M 

7 i . G H — y 1 4 1 ♦ 

396  0 

V'1  ''j.  G 3 — 't  1 4 1 ♦ 

3;  9 6 0 

3C=DUPEZ  '3 

3 9 7 0 

3:R=3  C-  FSC  ♦ 

3 9 '3  0 

7E=3C  ■■  F3  C ♦ 

3 9 9 fi 

VdL Z=4.1 88 

4 0 0 0 

vdl  □=vnL.z 

401  0 

G 1.1.1=  '■  PR+VCIL 

4 030 

DC!  3010  J= 

4 1 j 3 0 

GEM  1 1 ’ ...1  = G 

4 04  0 

FNECK  - 1 ’ J':* 

4 06  0 

3 EM  l , J =7 

4 06  0 

301  0 \ PM  1 ’ J':-  =E 

4 070 

DC!  3 060  J = 

4 08  0 

VEHGS C 1 ’ J> 

4 09  0 

306  0 VEHG3  •:  £,  J:> 

13. 

13. 

F3R 

F3B 

7 ♦ 3 R ♦ 9 E ♦ 3 C + 6 . 3 8 3 1 ♦ 3 R ♦ 

3 j ( 1 j ♦GT  Z ') 

1 ’34 

EM  '1  ’ J'1  .-'1  000. 

=FMECK  1 1 ’ J:-  ♦.  01746339 
EH  r'  1 ’ J:>  13. 

PM  c 1 , J 5 1 £ . 

1 ’ NV 

=VEHG3 '1 ’ J'  1 000. 
=VEHG3 1 3’  j 1 ♦33 . 3 


C ♦ h 0 + 1 1 1 3: 0 C k ♦3R^  1 DM: 


C ' 


D - 8 


— ‘.-4  — 


PRC 


c 


4 1 0 0 
4110 
4 1 3 0 
* 4 1 3 0 

414  0 

4 1 5 or 
4 1 p or 
417  OC 
4 1 0'  0 
.1 1 9 0 
4 3 0 0 
431  0 
433  0 

433  0 

434  0 
43"?  0 
4?P  0 

437  or 

438  OC 

439  or 

4 3 0 0 
431  0 
4330 

433  0 

434  0 
4350 
4360 

437  0 

438  0 
4390 

44  0 0 
441  0 
4430 

443  OC 

444  OC 
4450 

446  0 

447  0 
4480 
4490 

45  0 OC 
- 451 OC 

453  OC 

453  OC 

454  0 C 

455  OC 

456  OC 

457  0 
4580 
459  0 


SFUT = 0 . 

PUT = 0 . 

CONTINUE 

RETURN 

END 

THIS  SUBROUTINE  IS  A GENERALIZED  TABLE  LOOKUP  AND  INTERPOLATION 
ROUTINE  WHICH  IS  CALLED  BY  OTHER  ROUTINES. 

SUBROUTINE  LOOKUP 'A, FUN  - NRTS h B > 

D I MENS  ION  FUN  CP,  5 0- 
DO  3190  J=1-NPTS 
3 1 9 0 I R e F IJN  >'  1 i . J :>  . AT,  A'1  GOT 03 8 0 0 
33  0 0 IP  i-'.j.  RCu  1 J---3 

K=  J - 1 

B=  >•  A-FIJN  i t - k > :■  ♦ ('FUN  ('3  - j ■ -FUN  CP . kJ  CFIJN  C 1 - ,J -FUN  1 1 - K>  +FUN  >'  3 - k > 

RETURN 

END 

THIS'  SUBROUTINE  CALCULATES  TMR  NECK • TORQUE  AS  A FUNCTION  OR  THE 
NECK  ANGLE  AND  RELATIVE  VELOCITY. 

S U B P 0 U TINE  FOR C E T H C Y - T N F f k > 

COMMON  NF  Ck  NRN,  FNRfk  >:p,  34>  - DC  N 
DOUBLE  PRECISION  Y (PC 
TNECK  = 0 . 

TDhMF-  — DUN*  1 V (.3  1 — v 1 3 1 1 
VREL=Y  CP':-  -V  -'  S '- 
TREL=  Y < 6 ■-  — Y ( 7 

f A L L.  LOO  k IJ  P T R E L - F N E C k'  - N P N - T N F f k 
IF  CTREL.  GT.  0.  0.  AND.  VRF'L  . LE.  0. TNECK =-TNFCk> 0 . 4 
- IFCTPEL.LT.  0 . 0 . A N D . V PEL  . GE.  0 > T N E C k = - T N E C k ♦0.5 
T N E C k = T N E C k + T D A M F' 

RETURN 

END 


THIS  SUBROUTINE  OBTAINS  THE  CRASH  PULSE  G-'S  AS  A FUNCTION  OF  TIME. 
SUB R O'  IT  I ME  DEC  EL  1 T < GS> 

C 0 M M 0 N V EH-  N V - V E H G S - 3 - c'  0 > 

C AL  L LOO k IJ p c T ^ V EHGS  ■>  N V . G S 

RETURN 

END 

THIS  SUBROUTINE  COMPRISES  THE  MAJOR  PART  OF  THE  ,,PAC"  PROGRAM.  IT 
EVALUATES  THE  KINEMATICS  OF  THE  DEPLOYING  AIRBAG-  CALCULATES 
THE  FORCES  THE  BAG  APPLIES  TO  THE  PASSENGER-  CALCULATES  THE 
BAG  VOLUME  AND  PRESSURE-  DETERMINES  THE  GAS  GENERATOR  FLOW 
CHARACTERISTICS:  AND  CALCULATES  THE  IMPULSIVE  LOADS  ON  THE 
BODY  DUE  TO  BAGS  LAP. 

S'  IJ  B R 0 IJ  TINE  BA  G S'  I J B - X - Y - D Y - T N E C k - F T H - F X - F T T - G S > 

REAL  N E C k C H k - K P N - L T - L F 

C D M M 0 N •- S EATFRI C •- N S' F - S U N - SFUT  - PEL S F - S F N C 3 , 3 4 > 


D - 9 


- 10  - 

PRC  OP-'?  0.-'8£ 

4 6 0 0 C □ M M □ N , - S T F 0 F'  C E N S T < S T F ■'  £ ' £ 4 

461  0 L □ M M □ N L F C3 F’ U E N C ’ i..  F 1 c'  i c 4 ■' 

4 6 2 0 C n M M n N ••••' L N EEF’E S T N K F'  - S' L P ' P I.J T , p E L E P - K P N •:  £ - £ 4 

4 £ 3 0 C n M M 0 N M H N P R T , ■ 7 L * Z T « Z S'  - Z H < P T - P N - P H , P T 0 P H « ‘ < £ Z « V £ Z - W E ? E 0 < MH  i LT  ? L. F 

464  0 COMMON  M M R N D R T S TER M F' n C D R M P 'CHE S T T ' HERDT-  THFCN  THE 0 ' DELTTE 

4 4. c 0 C D M M □ N G R S'  FLO  N P G - G E N •"  £ « £ 4 > 

4 P £ 0 C D M M □ N G R S' DRT RTMDPt  P G Z " GTZ'  U - P N 1 * P N £ - P N ? 

467  n C □ M M □ N F R G D R T VC  1 * VC  £ * RV  • FS'R  n FSC  - FS'E'  ■>  R 0 * WBRG ' DROLL Z - THETRD  ' X 1 - V 1 

46?  0 COMMON, -MERGDRT..-  D MS  ■>  D I NF  * l.iS  OCR  , DIJPBZ  - VCGH < VC  GV 

4 £ Q i S-R  i SB'  S-L  ' S'  R F‘ 

4 7 0 0 C 0 M M 0 N M D R T - DU P E G H - D U P E G V - D P BGH-.  P P B G V ' V U P B G H ' V P E G H - V UPB G V ' V P B G V ' 

4 7 1 0 & G U PER-  G UP B V * G F F N , G P P V - V S: G M 

4 7£ 0 COMMON  MM  I S C ,-  FMi? ' PPP  - PP5  - GT  - FPN  - VOL  0 ' GW 

47? 0 COMMON "-PRPRM - PFT,  TNFTRTZ-  THFTRMZ , PG1  , FENCE'  GPB-  DPB-  ROT 

474  0 COMMON  MFRPRM  ■ S F ' F S' TBS'  ’ FCBS  < GUPB ' DUPE , DEC ' DPS  ' DS  C ' GST 

4 7 5 0 C 0 M M 0 N - T I M F S T E P - X S TOP'  PE l-  TRTi.PI N T 1 , p I N T £ 

4 7 £ 0 C 0 M M 0 N - M M PRPR M p p , V D L ' F F T C ' F p C , V S C ' l,J PE'  DROLL-  P P DLL  * X P T ' DTD P S 0 

477  0 COMMON  - MPP I NT  FFTH-  FPH,  THFF  ' THEL  ' EPH'  VOT  ' WIJPB'  FOOT 

4 7 6 0 D 0 U P L E PPE C I S 1 0 N V -'  8 - E ' R ' P ' E ' R 1 ' E 1 ' C 1 - X £ ' V £ 

4 - 8 l.i  2.  F 0 F’  M R T • 1 M — < 

4800  DIMENSION  DV.FU 

4810  T HE T R T = Y <:  7"i 

4880  TRET RH=Y  >■£"< 

488  OC  CHECK  TO  SEE  IF  PR'S  ENG  EP  S'UEMRPINJNG. 

484  0 IF  ORBS'  OTHETPD-THETRT > . GT.  1 . 57>  GO  TO  500 

485 OC  CRLCULRTE  THE  SI. OPE  OF  THE  PRSSENGEP  TOPS'O. 

4860  OSL DPE  = TRN  <?.  14159  "'8.  + THFTRT"* 

4 8 7 0 I F 0 h B S 0 0 S'  L 0 P E . GT.  1 0 0 0 . 1 0 S'  LOPE  =1  0 0 0 . ♦ R E S'  0 0 S'  L 0 P E > 0 S L 0 P F 

486 OC  CRLCULRTE  THE  NEU  H-POINT  COOPDINRTES. 

4 8 9 0 X £ = X £Z-  ''  V 0 5 > - Y •'  8 > ♦18. 

4 Hi.  in  Yc‘=YcZ 

491  OC  CRLCULRTE  THE  WEIGHT  OF  THE  !.  INPESTPRINEP  PORTION  OF  THE  ERG. 

4 9 8 0 l.i IJ PB=l.ipR G ♦ oflfi+F S B - E 0 1 opt.lP B - D l J P B Z "•  ■•-'  0 8.  ♦F S C 1 > 0 R n + F S P > 

49  " 'J  WUPB  = I.IIJPB+  1 l.ipHb  — l.ii  (F'E  1 ^4  „ ♦ H ♦ 1 hi  i + Fs  E'  — F*  II  •’  •- s HF* 

494  OC  CRLCULRTE  THE  WEIGHT  OF  THE  PE'STPRINED  PORTION  OF  THE  ERG. 

4 9 5 0 i, IPB= I... i B R G ♦ B 0 ♦•:'!-  0 p P p-pi j P E Z I*  • - 0 £.  ♦F S C "« o R i;i + F S E > 

4 9 £ IJ  f >i F' E “ l.i F' 9 I.' h b — lilF'Fi '!  *4  „ ♦ ':■  H ♦ F^  U " R ^ 

49? OC  CRLCULRTE  THE  DISTPNCE  THE  UNPES'TPp T NFD  ERG  FPONT  HRS  MOVED 

498 OC  PELRTIVE  TO  THE  GROUND  IN  THE  HDPIZONTRL  DIRECTION. 

4 -I • i FiljF: F’ b H = DUP F'b H + V U PBbH^DE  LT R T + 1 b . *i „ c' ♦ b U R F' H ♦Fi E L T h T ♦ ♦ — 1 

5 0 0 OC  CRLCULRTE  THE  VELOCITY  OF  THE  UNPESTPRTNED  ERG  FRONT  PELRTIVE 

501 OC  TO  THE  GROUND. 

5 0 8 0 V U P BGH= V U P E G H + ? 8 6 . 4 ♦ G U P E H ♦ D E L T R T 

5 03  OC  DISTRNCE  THE  UNPES'TPP  I NED  ERG  HRS  MOVED  IN  THE  VEPTICRL  DIRECTION. 

5 0 4. 0 D IJ P E G V = D U P E G V + V IJ P P G V ♦DELTRT-1 93. 8^G U P P V ♦ D E L T R T ♦ ♦ 8 

505 OC  VELOCITY  OF  THE  UNPESTPRTNED  ERG  IN  THE  VEPTICRL  DIRECTION. 

5 06  0 V U P B G V = V IJ F' B G V - 3 8 6 . 44-GUPBV+DFL  TRT 

5 07 OC  COMPUTE  THE  GROWTH  OF  the  UNPESTPRTNED  PORTION  OF  THE  ER'T- 

5 08 OC  WITH  RESPECT  TO  THE  VEHICLE  IN  THE  HOPIZON.TRL  DIRECTION. 

5 1. 1 9 0 DUPEH= Y '■  8 1 ♦ 1 c' . — D U F F: i- H 


D 


10 


-11- 


p fir- 


51  one  Nfll'l  IN  THE  VERTICRL  DIRECT  TON’ 

51  1 0 pijR  PV=DUEEGV 

5 1 £'  OC  TOTRL  IJ.R.  ERG  GPOWTH  i,l-R  TP  THE  VEHICLE  IN  THE  DIRECTION  THETRD . 

5 1 3 0 D IJ  P p = 3 Q R T f D I.  J R P R ♦ ♦ P + D M P P V ♦ ♦ 2 "> 

514 OC  C'RLCUL RTE  THE  INSTfiNTPNFOUS  LENGTH  OF  THF  "C"  PfiG  fiXIS. 

5 1 5 n SCO=SC 


5 1 6 0 
5 1 7 nr 
5 1 8 0 
c 1 9 nr 
5 3 DO 
531  0 
533  DC 
■3  0 C 

4 D 

5 or 

6 nr 

7 n 
p n r 
9 r 
n n r 
1 0 
3 or 

3 OC 

4 0 

5 or 

6 0 
7 Of- 
O i'i 

539  or 
54  0 0 
541  n 
543  OC 


SC=DUPP.-"8. 

DD  THE  5RME  RPR  THE  ‘'fi"  fiXIS. 

5R=F3R-FCC^SC 

RND  NOW  THE  "P"  RX IS. 

0'  E = E 3 P F S C ♦ 3 C 

XC— — >:  3C-  r DM  3-D  INF  > > 

CRLCULRTE  THE  D I 3T fiNCE  THE  PESTER I NED  PfiG  FRONT  HRS  MOVED 
RELfiTI VF  TO  THE  GROUND  IN  THE  HOP  I ZONTfiL  DIRECTION. 

D R P G H-D P P G fi + V p P G H * D E L . T fi T 4 1 9 ? . 3 ♦ G R P R ♦ D E L T fi  T ♦ ♦ 3' 

CRLCULRTE  THE  VELOCITY  OF  THF  PESTER I NED  PfiG  FRONT  RELRT I VF  TO 
THF  GROUND. 

V P BGH= V P P G H + 3 8 6 . 4 ♦ G R P H ♦ D F L T fi  T 

DISTfiNCE  THE  PESTER I NED  PfiG  HRS  MOVED  IN  THE  VERTICfiL  DIRECTION. 

D P p G V = D P p G V + V P p G v ♦PFLTfiT-19 3 „ 3 ♦ G P P V' ♦ D E L T fi T ♦ ♦ 3 

'V E L 0 C I T Y OF  THE  RE S TPfil  N F D P fi G I N THF  V FF’TI  C fi L DIR E C T 1 0 N . 

V P P G V = ' ' P p G V - 3 -3 6 . 4# G P P V ♦DFLTRT 

COMPUTE  THE  GROUTH  OF  THE  PESTPfiTNFD  PORTION  OF  THE  PfiG 
MITH  RESPECT  TO  THE  VEHICLE  IN  THE  HOP I ZONTfiL  DIRECTION. 

DR PH=Y ♦13. -DRPGH 

NOW  IN  THE  VERTICfiL  DIRECTION* 

DP PV=PEPGV 

TOTfiL  R.  PfiG  GROWTH  U-R  TO  THE  VEHICLE  IN  THE  DIRECTION  THFTfiD. 

D R p = 3 0 R T D R P H ♦ ♦ 3' + D P P V ♦ ♦ 3'  > 

d;  RND  Y COORD.  DF  H- POINT  IN  XRRIME*  YF’P  I ME  COORD.  SYSTEM 
XH=  •-  >"  Y3—Y 1 :>  ♦S  I N THFTfiD)  + < X3  -X 1 > OS  ' THFTfiD  ■>  > 

YH=  •:  (Y3-Y1  - CX3-X1  > ♦TfiN  - THFTfiD)  ) ♦ L+COS:  < THFTfiD)  ) 

Cfil. CULfiTE  THE  X fiND  Y COOF’D I NfiTES'  OF  THF  NECK  PIVOT  IN  THE 


543  nr 

544  0 

545  n 

546  or 

547  or 

548  OC 

549  i'i 
S'  5 0 0 
551  0 


554  DC 
'*  5550 


556  0 


'■'/-'■r  SYSTEM. 

XNF C K = X H - R N ♦ 1 3 . ♦ S I N - T H F T fi T - T H E T fi D ) 

Y N F C K = Y H + R N ♦ 1 3 . ♦COS < THF  T RT-THET  fi D ) 

IN  THE  NEXT  SEVERfiL  STfiTEMFNTS  THE  INTERCEPTS  OF  THE  filPPfiG 
RND  THE  TORSO  RND  HERD  fiRE  CfiLCUL.fi TED.  THESE  INTERCEPTS  fiRE 
D E S I G NR TE D P Y X fi  * Y fi  fi N p X P * Y P IN  THE  fi I P P fi G COO P D INRTF  S Y S TEM, 

P = Y 3 - Y 1 - 0 S L 0 P E ♦ '■  X 3'  - X 1 ) 

fi=Sfi+#£+ LCDS ' THETRD  ' +0S LORENS  IN CTHFTfiD)  ) ♦♦£ 

+ S' C ♦ ♦ 3 ♦ ' SIN  - THETRD"'  -OSLOPE+COS  'THETRD)  ) ♦♦£ 

D=£.  ♦P^SC  + ^3^  '"OS  LOPE  ♦COS  0 THFTfiD)  -SIN  'THE  TfiD  > • 

F = P^^3t  '-SC  ♦♦3)  — S fi  ♦ ♦ 3 ♦ ('SO**3'  ♦ 'COS  'THETRD  1 +OSLOPE>S  IN  •:  THF  TfiD"'  '■  ♦♦£ 
TEST  FOR  SIGN  OF  D I SCR  I M I NRT'E 
I E 1 D ♦ ♦ 3' — 4 . ♦ R ♦ E ■'  3*4*4 
3 T VO 1=0. 


5570  GO  TO  14 

5 5 8 0 C F'Ffi L D 1 5 T I N C T POOT S'  ■'  TOP S 0 L INF  P fi S S' E S'  THP 0 U G H PfiG ) . 

5590  4 D I S C= > D*^3-4 . ♦R^E ) ♦♦ . 5 


D 


11 


- 18  - 


PHL 

08-' "3  0 -"88 

5 6 ij  0 f 

56 1 0 
5 6 HO 

B HO  INTER C E P T POINTS*  X ft * X E ft N D V ft * V E 
Xft=  '"-D-D I SO  ••••'  Oc'.  ♦ft'1 
X E=  (-E+E I SO:'.'1  .••-  <8.  ♦ftj 

563  0 

564  0 

565  0 
5 6'  6 f i 
567  0 

V ft  = '"E+Xft*  ."OS  L OPE  ♦ O' OS' ' THE  THE  :■  -3 I N < THETftE  '•  > > (COS  •:  THETftD:'  +DSLOPE  ♦ 
6 SIN  ' THETftE'  "' 

I F >:  S' ft , LE . HE'S  (VO'1  '>  V ft = HE  S'  •:  V ft '« .-"VO*  •:  S ft- . 0 0 1 > 

V E = ■" B + X E ♦ C 0 S LOPE^C OS '"THETftD > - S I NOTH E T ft E > > 1 *" 0 0 S ■'  T H E T ft E > + 0 S L. 0 P E ♦ 

S' IN.:  THETftD':-  > 

f-.  {•-■!  n 

I F •'  s ft  „ LE  . ftE  S ’ VP  1 1 VE=ftES  1 V E > •-■' YE'*  ( S ft— . 0 0 1 :• 

569  0 
5 7 ij  0 
571  Of 

573  0 
5730 

574  0 
5750 

IE  •: THETftT-THETPE . GT.  0.  . PiNE.Vft.LT.  0 „ ':>  '30  TO  518 
I E ■'  THFTftT-THETftE . L T „ 0 . . ftND . Vft . 9T . 0 . 80!  TO  51  0 
ftEST =D I STfiNO'E  FPDM  POINT  ft  TO  POINT  p. 

ft  J-  T — \I:.PT  '*  Xft — : E:  '*  1 V ft — V B 1 ^ ♦ 0 1 

IF  '•  THFTftT-THETftE  1 6*  6*  7 
6 VP=VE 
VN=Vft 

5 7 6 1 j 

X P = X E 

577  0 
579  0 
5790 
58  0 0 

XN=Xft 
GO  TO  8 
7 V E’ = V ft 
VN=VE 

5310 
588  0 

583  Of 

584  Of 

585  0 

XP=Xft 

XN=XE 

CftLCULftTE  THE  EISTftNOE  FROM  THE  H-POINT  TO  THE  EftG  SUPEftCE 
ft LONG  THE  TORSO  LINE. 

8 P p ft G = 3'QPT  ( ' ■' N - X H j ♦♦  8 + •" V N - V H 

586: 0 
587  0 
590Of 

I E f V r j . L T . V H '■  P B ft  G = - P E ft  G 
I E (P EftG . GT . PN+ 1 8 . '■ GO  TO  3 

COMPUTE  THE  X ftNE  V CO  OPE  I NOTES'  OP  THE  POINT  OP  FORCE  ftPPL  I OPTION 

533  00 

ON  THE  TORSO  IN  THE  X"'-Y"'  SYSTEM. 

5 3 0 1 j 

I F HE'S  T +PEHG- 1 c . ♦ P N :>  9 * 3 . i 1 1 

591  0 
5 98  o 

•4  YUT=YP 
XUT=XP 

593  0 

594  0 

595  0 

596  0 

597  0 

VLT =YN 
XL.T =XN 
GO  TO  11 
10  VUT=VNECK 
XUT=XNECK 

598  0 

VLT=VN 

5 9 9 1 j 

XLT =XN 

6 1 j 0 1 j 
6 Ij  1 Ij 

6 08  OC 
6 03  0 
6 04  00 
6 05  00 

1 1 VFT=  '"VUT+YL.  T '<  ..-'0  „ 

XFT=  '"XUT+XLT  :>  .-"0. 

F'FT =D I STft'NC E FROM  H-POINT  TO  POINT  OF  FORCE  ftPPL ICftTION  ON  TORSO 
PFT=SOPT  c'XH-XFT':.  ♦♦8+  •'  VH-VFT"'  ♦ ♦£':• 

PSLOPE=THF  SLOPE  OF  ft  LINE  PFPPFNEICUL.ftP  TO  THE  TORSO.  THIS  IS 
USED  TO  FIND  THE  X"'  ftNE  V"'  COOPE I NOTES'  OF  THE  POINT  MHEPE  THIS 

6 06  00 

LINE  INTERSECTS  THE  ftlPEOG  SO  THOT  THE  DEPTH  OF  EftG  PENETPftT I ON 

6 07  00 

IS  ON  EE  COMPUTED. 

6 08  0 

PSLOPE=-  (COS  (THETftE:'  +OS  LORENS  IN  (THETftD.'-'  > (DSLOPE+COS  (THETftD':' 

6 09  0 

C-SIN  '"THETftD > 

~ 
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"3  iV' 8 3 


6 1 0 0 ft  t = 1 . - S C ♦♦£+ P S L □ F'  E ♦ ♦ 8 S ft  ♦ ♦ 8 

6 1 1 0 B 1 = 8 . ♦PC  LOPE  S'  ft  ♦ ♦ 8 ♦ V FT-P  S'  L □ P E ♦ X F T > 

ft- 1 8 0 C 1 = V FT-P S L □ P E ♦ X. ET>  ♦ ♦ 8 S ft ♦ ♦ 8 - 1 . 

El  3 OC  THE  X"'  FIND  Y"'  CO  OFT I NftTES  of  this  point  ftPE  xg  and  yg. 

El  4 IJ  XU  — 1 — B 1 I 'FT  1 B 1 ♦ ♦ c — 4 M ♦ft  1 ♦!,_.  1 * • .*•*  L~'  u ♦ft  1 

E.  1 3 0 Y ft = P s LOPE  ♦ X G - X F T "<  + Y F T 

El  6 Of  THE  BAG  PENETRATION  IS  "BP". 

E.  1 7 IJ  F'  P ~ M P T ft  V '■■■’  E 7 — X I t *'  ♦ ♦ i~’  *+■  Y E T — Y ft  1 ♦ ♦ c' . 1 

ft 1 8 OC  COMPUTE  THE  TORSO  ftpEft  OF  INTERCEPT  IN  THE  FOLLOWING  STATEMENTS, 

ft  190  ftCHORDR=0. 

ft 3 0 0 AC  HOP'D N = 0 . 

ft 8! 0 ftCHOR DHP= 0 „ 

E88  0C  SOLVE  FOP  THE  ARC  SIN  OF  YUT.--SR  AND  YLT-SA. 

ftc'  7'  i.1  1 .7  ft  _ YUT Sft=HT hN  1 ; YUT  ft  1 1 ’ I,|RT  ft  1 „ — 1 y UT •••'  •::  ft  1 ♦♦c!  •'  ■'  1 

ftc'4  m ft  ■ y'LT  _ H = ft T H N ! 1 YL7  ft  1 ■ i *F  T < 1 „ — • y L T ■-  s ft  1 ♦♦c'.1  •'  1 

ft 8 5 OC  BAG-  TLINE  AND  ft CHORD  ARE  INTERMEDIATE  VALUES  REQUIRED  EOF1  THE 

ft 8 ft  OC  AREA  OF  INTERCEPT  CALCULATION. 

ft  8 7 0 IP  ft  X P.  LT.  0.  . A N D . Y N E C k'.  GT.  S'  A A C H 0 P D P = S C S A ♦ <3 . 1 4 1 3 A ♦ S A ♦ ♦ 8 

ft  8 8 0 S ■ 8 . - ft  Y IJ  T ♦ S Q ft  T >"S  ft  ♦ ♦ 8 - Y IJ  T ♦ ♦ 9 :<  + s ft  ♦ ♦ 8 ♦ ft  S Y I J T 'S  A > > 

ft  S'  9 0 B A G = S'  C ft  8 . ♦ S'  ft  > ♦ ft  1 Y U T ♦ S Q F'  T ft  S ft  ♦ ♦ 8 - Y I J T ♦ ♦ 8 > + S ft  ♦ ♦ 8 ♦ ft  S Y l J T 3 ft  :>  - 

ft  3 0 0 j < Y L T ♦ S Q R T S'  ft  * ♦ 8 - Y L T ♦ ♦ 8 + S ft  ♦ ♦ 8 ♦ ft  'r  Y L T S ft  •> 

ft 3 1 0 T'L  I NE=  (XL  T-  ( SN-YNE C K>  •-  >:'YN-YNECE>  ♦YLT>  ♦YUT  + ( YN-YNEC k > (. 8 . ♦ C YN-'y  NEC  k 

ft  3 3 0 '•  > ♦ Y IJ T ♦♦8-  f’ X L. T 

ft  ? 3 0 C - ft  X.  N - X N E C k > •-  ft  Y N - Y N F C k ♦ Y L T ♦ Y L T - ft  X N - X N E C k :■  ft  3 . ♦ < Y N - Y N F C k > > ♦ Y L.  T ♦ ♦ 8 

ft  34  m I F i'  YN . LT.  0.  > AC HDRDN=SC -"S  A+  ftYL T ♦SORT  •" t A ♦ ♦ 8 - Y L T ♦ ♦ 8 > +S  A^^3^ 

ft  3 3 0 j AS  Y L T S ft + 3 . 1413Q*Sft+^ 8 2 . '■ 

ft 3ft Of  SOLVE  FOR  THE  TORSO  AREA  AND  VOLUME  OF  INTERCEPT, 

ft 3 7 0 TAD  I =BftG-TL I NE +ftCHDRDP+ftCHDRDN 

ft 38  0 TVO I =7  AO  I ♦MB 

ft. 3 9 OC  CALCULATE  THE  HEAD  INTERCEPT  I.IITH  THF  BAG. 

ft4 0 Or  FIRST  CALCULATE  SLOPE  OF  HEAD. 

ft  410  14  IP  CABS  '' THE7  AH-THETftD '<  . GE.  3 .141 393-' 8 . • GO  TO  3 1 8 

ft 480  HSLOPE=  + TAN  ft'?.  14138  "8.  + ftTHETAH-THETAD-'  > 

ft.  4 3 0 ft  H = H S'  L 0 F'  E ♦ ♦ 8 S ft  ♦ ♦ 8 + 1 . ■•-'  S C ♦ ♦ 8 

ft 4 4 0 b H = 8 . ♦ H S L 0 P E ♦ ft  Y N E f k'-H S L. 0 R E ♦ X N E C k S A ♦ ♦ 8 

ft  4 3 0 C H = ft  H S L 0 F"  E ♦ X N E C k ♦ ■ H S L 0 P E ♦ X N E C K - 8 . ♦ Y N E C k + Y N E C k ♦ ♦ 8 > S A ♦ ♦ 8 - 1 . 

ft  4 ft. 0 IF  ft Ei H ♦ ♦ 8 — 4 . ♦ ft H ♦ C H •>  13-  17-  17 

ft. 4 7 0 13  HftOI  = 0. 

ft 48  0 GO  TO  3 0 

649  0 1 7 D I SCH=SQRT  ftBiH^P-4 . ♦ftH^C H> 

ftSOOC  THE  COORDINATES  OF  THE  HEAD  INTERCEPTS  I.IITH  THF  BAG  ARE: 

ft  3 1 0 X Hft-  ft-BH-DI  S C H ■•-'  <'8.  ♦AH 

6 3 c'  Li  '■  HF  ~ 1 — BH+B I !;  L H 1 1 . ♦ftH  1 

ft-  3 0 y H ft — H L 0 F"  E ♦ 1 Hft  — : : N E L k 1 ♦ V N E i..'  I7 

634  IJ  I P ft S H • L.E  ■ V Hft  > YHft  — ft  F' S ft  V Hft V Hft ♦ ft  N ft—  „ l.i  i.i  1 j 

633 0 Y HB=H S L 0 P E ♦ ft X HB- X N E C k > + V N E C k 

ft  3 ft  0 I F ft'SA . L E . Y H B > Y H B = A P S ft  Y H B > Y H B ♦ •'  S ft  - . 0 0 1 1 

6370  IE  ftTHETAH-THETAD)  19-19-80 

ft  3 8 0 C D I ST I N G U I S H BETWEEN  U P PER  ft N D L OWER  V A L U E S OF  HEAD  INTER C E PT  FT S . 

6390  19  YHF-YHB 
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66  00 

Y H N — Y H R 

6 6- 1 0 

XHP=XHB 

663  0 

XHN=X HR 

h?  P?  -J  J J 

GO  TO  £ 

6-  6 4 0 

3 IJ  Y H P = Y H R 

6 6 7 0 

Y H N ~ ''i'1  H B 

6 6-  6 0 

XHP=XHR 

6 6, 7 M 

XHN=XHB 

668 0C  DEFINE  X RND  Y COORD I NRTFS  DP  TOP  OF  HFRD  IN  X'-Y"'  SYSTEM. 

68  9 0 81  ' TnPH-VNFCP -1  £,  ♦ < PTDPH-PN^  ♦SIN  1 THETAH-THETRD> 

67  00  YT0PH=YNECK+13.  * >:  RTOF'H-RN:.  ♦CDS  'THETAH-THETAD  •' 

671 OC  SEE  IP  BOTH  NE  CP  PIVOT  RND  TOP  OP  HERD  RESIDE  WITHIN  BAR. 

8 7 £ 0 TOPCHP  = •••■  T 0 P M ♦ ♦ P s c ♦ ♦ 8 + V T 0 p H ♦ ♦ £ ..  •-  v R ♦ ♦ 2 

8 7 3 0 N E C K C H 1-  = ' N P C P ♦ ♦ 2 • ' '*  C ♦ ♦ 2 + V N E C P ♦ ♦ 2 S R ♦ ♦ 2 

674  0 I P • TOPCHP  RND . NEC  PC  HP  , L T . 1 , > GO  TO  37 

677m:  SEE  IF  THE  NFCP  PIVOT  PT.  RESIDES  WITHIN  BAG. 

676 0 I F ' NEC P CHP.LT. 1. > GO  TO  3 0 


677  or 

SEE  IF  ANY 

F'RPT  OF  HERD  RESIDES  WT  THIN 

B H N . 

673  0 

IF  > TOF'CHP  . i 

GE . 1 . > GO  TO  17 

6 79  OC 

NFCP  PIVOT 

NOT  WITHIN  BAG  BUT  HFRD  TGF 

IS. 

6 8 0 0 

YUH=YTORH 

681  0 

XUH=XTOPH 

6-88  0 

YLH=YHN 

6 8 3 0 

XLH=XHN 

6. 8 4 0 

GO  TO  38 

687  OC 

NFCP  PIVOT 

WITHIN  BAG-  BLIT  HFRD  TOP  IS 

NOT. 

636  0 3 0 YIJH=YHP 

6870  X I J H = H F' 

6 o 8 0 Y L H — Y NF LB 

6890  XL.  H=XNECB' 

690 0 NO  TO  38 

691 OC  BOTH  NFCP  PIVOT  RND  HERD  TOP  WITHIN  ERG. 

693  0 37  YIJH=YTOPH 

6930  v;uh=YTOPH 

694  0 YL.H=YNFCP' 

697  0 :-'LH=YNFCP 

696 OC  NOW  CRLCULRTE  THE  V RND  Y"  COOPD I NRTES  OF  THE  POINT  WHERE  THE 

697 OC  FORCE  IS  REPLIED  TO  THE  HERD . 

6 9 8 0 8 Y F H = 1 Y U H + V L H 1 ■'  2 . 

6 9 9 0 7 FH=  1 X UH+ X L H :«  - 2 . 

7 0 0 OC  CRLCULRTE  THE  PENETRATION  OF  THE  AIRBAG  BY  THE  HERD. 

7 0 1 0 H P S L □ P E = - 1 H S L 0 P E 

7 030  H R 1 = 1 . S C ♦♦3+HP S L 0 P E ♦ ♦ 3 S R ♦ ♦ 2 

7 0 3 0 H B 1 = £ . ♦ H P S LOPE S R ♦ ♦ 2 ♦ •" V FH-H P S L OPE ♦ X F H • 

7 04  0 Hi?  1 - - YFH-HPSLOPE^P'FH  ' ♦♦p.--SR*^3- 1 . 

7 0 7 0 X N H — : — HB  1 + S wpT  *'  H B 1 ♦ ♦ — 4 . ♦ H R 1 ♦ H i._  1 ■ .*  • L~'  B ♦ H R 1 1 

7 0 6 0 V G H = H P S L.  0 P E ♦ ' X G H - X F H > + Y F H 

i U 7 U E'  F'  H — ;■  i ; 1 F*  T 1 : F H — N H 1 ♦ ♦ 6 + 1 V F H — V N H ■'  ♦ ♦ c'  1 

7 08  OC  C OMPUTE  ARCS I NS  YLH  SA  RND  YUH  - SR. 

7 0 9 0 4 0 R '? V L H S R = R.T R N ' V L H S R "•  •'  S 0 PT  *T  . - Y L H S R ♦♦£  :•  > > 
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14 


15 


PRC 


iy  U-yy 


7 1 0 0 

R S Y U H S'  ft = ft  T ft  N •:  C Y Li  H 

7 1 1 0 

BRbH  = SL-"  1 £ . ♦ : R i ♦ 1 

- 1 8 0 

o C YLH^SUF'T  C S R++C— 

1 3 0 

HL I NE= CYL H- CYNECK 

714  0 

f!  EOF  -YHP>  > ♦Y(  !(-!♦♦£: 

7 1 5 'J 

C CYLH-  NECK- VHP 

~ 1 6 0 

L K — v H F'  i ♦ Y L H ♦ ♦ c' 

7 1 7 0 

IF  .""'HP „ LT.  0.  . RNP. 

7 1 8 0 

f;,-  . V’UH* S OPT  f'SP*4g 

7 1 R 'j 

IF  CRCHDPPP.  NE.  0.  -■ 

►v  I..IH+  .>:he  CK-YHP:.  /' 


one  -SOLVE  POP  THE  HERD  PEER  RHP  VOLUME  OF  I NTEPCEF'T . 

.1  0 HAD  I =BRGH-HL  I NE+ACHDF'PHP 

SO  5 0 HVO I =HRO  T ♦lilH 

8 0 C " V □ i ■ ■ - v o i i i m e OF  BR G I N T F P C EPT. 

4 PC  CRLCULRTE  THF  TDTRL  VOLUME  OF  BODY  INTERCEPT  MITH  THE  ERG. 

50  VOI=TVni+HVOT 

6 0 V C = S R ♦ S 0 PT  < 1 . - V C ♦ ♦ ? 0 C ♦ ♦ 3 T 

□ iL_ = 4 „ l 5 E:  ? ♦ S R ♦ :v  p ♦ " f + 6 , £ 8 3 1 ♦ S ft  ♦ S'  C ♦ ft  0 + M S 0 C L * V C ♦ C P M S'  - C S'  C + S'  C :<  - 
I F ■:>: -GEN  c 1 i 8 ft  0 , ft  0 - 1 0 1 
ftp  RR=ftTMnP 


7 3 0 0 
731  or 


733  0 

734  OC 

735  0 

736  0 


Ml 


7 4 0 0 
741  OC 
74-30 

743  OC 

744  0 
7450 
7 4 6 0 
7470 


’48  0 


74ft  OC 
7 5 0 0 


GD  TO  113 

COMPUTE  GRS  FL.OM  INTO  ERG 
1 0 1 c R L L Gft  S'  I N < Y , p I N * 

GUI  =GU+ Q I N ♦ PEL  T ft  T 

COMPUTE  NEI-I  TEMP.  RNp  PRESS.  DUE  TO  NET  GRS  GRIN  IN  ERG 
b T 7 = 1 b l.i!  ♦ b T + p I N ♦ b T 7 ♦ D E L.  T H T 1 b U 1 
PNIJM=U*GT7  + GU  1 
P El  7 = P n | j m v OLD 

COMPUTE  NEM  GRS  PRESS.  RNP  TEMP.  DUE  TO  POLY TROPIC  COMP.  OF 
Pftft=  '"PNUM.--VOL > ♦♦PPN-" PR7  + + CFPN-1  „ > 

GT ft=GT7*  cp ei ft  - PR7>  ♦♦  C FF'N- 1 . > FF'N  :• 

ERG  VENTING  COMF'UT  RT I ON'S  ’ FIRST  CRLC.  PRESS.  RATIO  RCPOSS  V 
P p 7 = r T M 0 P P R 8 

TEST  FOP  CHOKED  FLOWS  ALSO'  IF  PP7>  1 , EYPRSS  QEYH.&-  SET  Gl.i=G 
IF  CPR7.  L.T.  PRGO  GO  TO  108 
IF  CPF'S.  GE.  1 . 1 GO  TO  110 

FM 1 =VC  1 ♦S  OPT  fPp?++  (d . -•■'  P N 1 > — p P 7 ♦ ♦ •"  -'  PN1  + 1 . "<  ,'PN  1 ) 

GO  TQ  1 Oft 
1 08  FM1=FM8 

COMPUTE  EYHRIjST  FLOU  RNp  PES I PURL  GRS  HEIGHT 
1 0 ft  0 E Y H = S P P T c 778.  ♦ P N 1 > <" P N 1 - 1 „ } ">  ♦ R V ♦ P R ft ♦ F M 1 S 0 P T C IJ ♦GT ft 


5 


0 


754  r 


I Hi  I 


i7C 


I ftf 


GU=GUl-PE7H*DELTftT 
IF  CGU.  LT.  0.  ■'  GU=0. 

GO  TO  1 1 1 

110  GI4=GM  1 

COMPUTE  PRESS.  RNP  TEMP,  OF  GRS  RFTEP  VENTING 

111  PRTI0=GI.VGU1 
F’ft=p  Rft^pftT  I 0^+F‘N  1 
GT=GT 8+F'RT  1 0++  c PN 1 - 1 :< 


D - 15 


♦ <.  YN 


VO  I 


EYPRNS . 

ENT 
1.1 1 


- IS  - 


PRC 

oy. ■■-8  !.!,"•  83 

76  o nr 

COMPUTE  PRESS 

76 1 nr 

ESPPNDFD  THIS 

7 6 0 n 

PP6=PR8  PR6 

7 £ P 0 

IF  'PP6.LT.1. 

7 6 4 0 

FpN=pMP 

7 6 5 i*i 

GO  TD  l 13 

766  0 

1 1 0 

FPN=PN'3 

767  nr 

COMPUTE  FR 8 p 

7 p 9 n 

1 1 3 

PG1 =Ph-RTMOP 

7 h 9 ! - 

DROLL =EP0L LZ+ 

7 7 1 j n 

KPT COS  LOP 

7?i  n 

DTDPSD=SC+KPT 

776  0 

IF  CDUPE.  GE.  0 . 

773  n 

P ROLL =7' OPT  -I  •'  K 

"■74  n 

f;  -THPTRT':< 

776  '1 

IF  CPPOLL. GT. P 

7 7 6 1 - 

D . j b H=D  w i ? H 4*  y c j -i 

7 7 7 n 

V S G H = V S 0 N + 3 8 6 

7 7 P,  0 

flS  R V=DS  G V + r R 

7*  7 Q n 

VSGV=VS8V+386 

7 8 Ci  o 

DUbH=V 1 ^ 1 l „ 

781  n 

VC'GH=V  ■'  1 •'  ♦ 1 0 . 

7880 

1,867= P POL  L+  ■ r 

7 8 3 0 

VCG'V=-PP0L.  L>V 

7 8 4 n 

IF ' DUPE . LT.0. 

786  n 

SGH= VC  GH 

786  0 

VS GV=VCGV 

7870 

riPRH=DCGH 

788  n 

DSGV=DCGV 

7 p Q i'| 

1 14 

D CH=D  C GH-D  S'  G 

7 9 0 n 

PSC'V=DC  GV-DSG 

791  n 

D S C = D S C H ♦ C 0 S'  < 

7 R 8 Ci 

v s r h = v r g h - v s g 

793  c 

fJ  “■  p if  i>  fi  j~  j~  fi if  (i  *r-  j “ 

'-■i  4 ! ! 

VS  C= 'VS  CM  ♦CO  S'  •" 

796  0 

C RL  L CFOP  CDSC 

786  n 

FC  ES'=F  C ES  +CDP 

797  0 

IF CDPE-DTOF SO 

798  n 

115 

DP S = 0 . 

798  n 

DPC=  0 „ 

3 n n ri 

FSTES-O,. 

3 0 1 0 

PPOLL-SOPT  <:'K 

3 0 8 n 

GO  TO  118 

3 0 3 Ci 

1 17 

DP 8=  ' DROLL .-"0. 

3 04  0 

IF  < DROLL 70.  . L. 

3 06  n 

DPS-DPC -DS 8 

3 Cl  6 n 

CALL  ST FOP CDP 

3 ci  7 ci 

IF  ('DUPE.  GF.  0. 

3 08  o 

FSTES-FSTES+S 

3 0 9 Cl 

LETRT  i 

J0O 


V 0 v - V P F G V '>  ♦ S IN  1 T H 


D 


16 


17 


PRO 


I..IH  M 8 


S 1 o ;'i 

8 1 1 n 
8 1 £ 0 
1 !'i  r 

3 1 J-  r'i 
8 1 5 0 
8 1 8 nr 

81708 

P i p 0 


1 1 8 G 7 --  '■  FL  E i “7 T E ''  ■<  1 Zi  ♦38.  87 

GSTH=G5  T ♦CDS  •:  THE  TRT  ■ 

G 3 T V = - G T ♦ 3 T N '•  T'HFTRT 

IF  THE  ERG  PPESSUPF  18  MFGRTIVEG  CRLL  IT  0. 

PG=PG1 

T F '"PG  . L T „ 0 . PR-  0 „ 

IF  THE  TOP 80  18  NOT  IN  CONTRCT  1.1 1 TH  THE  ERG 
CRLCULRT I ON. 

I F >■  B ♦ ♦ 8 ~ 4 . ♦R#FI>  1 8 0 - 1 85  - 1 85 


180  F F T C = 0 . 


FPC=0. 


881  0 


BP=0. 


'?  8 '?  0 I 
884  0 
p p c 0 

88 E 0 


'I 


;;  ■■■’  0 1 i 

831  0 
838  0 


P 3 4 0 


835  n 

R3F  fi 


84  0 n 
8 4 1 0 
848  0 
8 4 3 0 


844  0 
8450 

8 4 1~  i i 

8470 
848  0 

3 4 P !'l 


I O i i 


"■1 


;g  n 


P 5 8 0 
P54  o 


GO  TO  188 

TMP  ERG  FORCES  RPF  CRLCULRTFD  IN  the  next  8f 
1 8 5 I F *•' P B R G . G t ,j  p N ♦18.  . OP.  X . LT.  G E N < 1 - 8 > GO  TO  1 i 
ENPH  I = F' G ♦ S R * •'  r > 8 » ♦ 8 I"' E T I-'  r R ♦ ♦ 8 + f ♦ ♦ p > 8 . ■'  '' 

C H E S T L. = S 0 P T < X N - •-i.it ■>  ♦ ♦ P + 1 V N - V IJ T "•  ♦ ♦ ? s 
FFT C = ;? . ♦ F N P h I ♦ c H E T L. 

IF  'IBP.  LT.  8.  ♦CHESTT-'  FFT C=FFTC ♦ CEP  - >:'8.  ♦CMFSTT) 
FPC=PG+l>!E  ♦ '“  HESTL 
1 88  FCHEST=FPC+FFTC 

IP  <;bh+*8-4 . ♦ r h ♦ r R ' 140>  1 35"  135 
1 35  I F •:  TO  PC  HE  . GT.  1 . . RND . NEC  EC  HE  . GT . 1 . GO  TO  1 4 1 
I F X . L T . G F N '•  1 " 8 > :■  GO  TO  14  0 
HEHDL=  ::--i.:.'PT  1 ' '■  UH  — XLH  ♦♦c  + 1 v !J H — v* I.H  1 ♦ ♦87 
E F'  H = P G ♦ l • ♦ H F R D 1. 

FFTH=8. ♦ENPH I ♦HE PEL 

IF  c'BPH.  L T 8.  ♦HFROT'i  FFTH=PFTH^  ■:  BPH  '8.  ♦HERE'T’ 
GO  TO  141 

140  FFTH=0. 

FF'H~  0 „ 

EPH=  0 n 

141  FHERB=FFTH+FPH 

F'HERE'  = mF'7  ' f y FH—  NFi  E j ♦♦,-'+  ; y‘F  H — YNE  E E 7 ♦♦8  1 
FTH=TNE  CE  -FHE  RE'*E  HFRE  ■ 1 8 . 

TPRN8T0P=— TNE C E 

F T T = - F H E R D ♦ P N ♦ C □ V < G 7 - V < 7 "•  ■ +TPR N 3 7 0 P - ■ P P 0 L L 
I F '■  PG  1 . GT.  0 „ 7 GO  TO  148 
G P E>  H = — G 3 3 8 . 8 
GPBV=0. 

G p B = G P E H ♦ C 0 S THETRB 
VPEGH=V  ■: 4"1  ♦18. 

VPEGV= 0 . 

•? U P E H = — i ? 3' '-! 8 . 8 
GUPEV  = 0 . 

G I J P Ei = G U P F H ♦ C 0 S 7 H E T p p > 
l,.,lUPEGH=V  '*  4 *'  ♦ 1 8 . 

VljPEGV=0. 
bO  TO  146 

1 48  I F ' POPE . GT . BUF  F 7 + 1 . GO  TO  1 43 


POP  ERG  FOPCE  CRL.CS . 
i 3'E  IE'  THE  ERE-'  FOPCE 

IVEPRL  3 TRTEMEN7 7 . 

0 


♦ F r E*  +PF  7 ♦F  CHEST'S  ■ 1 8 . 


O 


17 


- is  - 


PRC 

uy 

86  0 0 

gpp=  cfstbt  — PR  l ♦MS nr  o ♦fipni  l > .•-i 

SSI  0 

C-i  u PB=-PG  1 ♦ I.  ( 3 n C K ♦ D P 0 L L • • I.IBRG 

36  £0 

GO  TO  144 

8 6 3 0 

143 

GPP=  •: FSTBS^COS  C THFTRD-THFTRT: 

8 6 4 0 

H IJ  P'  p “ — F'b  j 4 V *4-  „ H4  1 H M — f!  1 J 1 4*  J q t * 

865  0 

144 

GPBH=GPF*C  03 CTHFTRP  1 

SS  6 0 

GPBV=GPE*S  TN  f.  THFTRD  •- 

S 6 7 n 

G IJ F pH-  G U P P ♦ f 0 ' c T HPT Rp  - 

S6S  n 

G IJ p p V = G U F' B ♦ 8 1 N . THFT'R D "< 

f,  Q n 

1 46 

I P DU F' B . GF . £ . ♦F 3 C GO  TO  15 8 

S 7 0 0 

148 

I P ' DPP . GF . £ . ♦ F:" C > GO  TO  1 55 

S71  0 

I P CBUPB „ GT . DUPBZ+ 1 . - GO  TO  1 1 

8 7'  £ 0 

IF  C D P PGM- D U P B Z ♦ C 0 3 <:  THFTRD-  . Gl 

8730 

PC  03  C THFTRD"' 

874  0 

GO  TO  165 

875  Cl 

158 

It !.J P' It H = — b ‘ " t r‘  a c‘ 

y 7 6 ? f 

GUPPV= 0 „ 

877  0 

G l j F'  p = G U P P H ♦ r 0 S'  C T H F T R D 

878  0 

VUPPbH  = y •'  4 ■*  * 1 c . 

879  o 

VUPPGV=0. 

88  0 o 

IF -PG1 . LT.O. > PG1=0. 

381  0 

GO  TO  149 

38?  o 

155 

G 3 T = '-  P C P 3 - F 3 T p 8 c Z 3 ♦ 3 8 . 8 ;■ 

88  8 0 

GS'TH=G3T ♦COS  CTHFTRTO 

334  0 

G 9 j V = - G 3 T ♦ 8 IN  CTHETRT 

8 8 5 i*i 

GPP=F8TP8^C  08  THFTRD-THFTRT 

386  0 

IFCGPp.LT.  - G S 3 £ . £ > GO  TO  15 6 

8870 

GPBH=GPP*C03  r THFTRD'1 

888  0 

G F'  B V = G P B ♦ 3 1 N -THFTRD  - 

389  fi 

VPEGH=Y  C4> ♦ 1 £ . 

39  0 0 

VPBGV  = 0 . 

8 8 1 f| 

IF  ''DPB.  GF.  DTOPSO:-  GO  TO  165 

898  0 

156 

bPBH=-G3  •••':.;!£,  £ 

893  0 

gpbv=o. 

8 q 4 0 

GPP= GPP  HUGOS'  -■  THFTRD'1 

895  0 

V P'  ?'  b H = V * 4-  > ♦ 1 ii  a 

396  0 

VPBGV=0. 

897  0 

GO  TO  165 

88  8 0 

158 

I F CDUPP . GF . £ . ♦FS  C GO  to  1 59 

p Q Q f | 

GO  TO  161 

9 1 j 0 o 

159 

G J T =DY  - 1 - :.:!£ . £♦!  OY  1 THE  T RT +P! 

901  0 

b S' TH=D Y 1 ■'  '3 £ o £ +PPOLL  ♦ - D Y - 

9 0 8 1 j 

G 8 T V = - P P 0 L L ♦ CDY  C3">  ♦SIN  CTHFTR 

9 0 3 0 

GP BH=G3TH 

8 04  0 

GPBV=G3TV 

9 05  0 

GFB=GST*COS  CTHFT RD-THFTRT  "■ 

9 06  0 

VPBGH=VSGH 

9 0 7 0 

VPBGV=V8GV 

9 08  0 

GO  TO  165 

9 09  0 

161 

IF CPG1 . GT. 0. > GO  TO  165 

3 y 6. 4 
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18 


19 


PRC 

08.-- 3 0888 

91  00 

IF  LFSTBS.  EQ.  0.  •'  GO  TO  167 

9 1 1 0 

G R B = F S' T E 3 ♦ C D S CTHETRP-THETRT:.  /l,IRB 

9 1 £ 0 

GRBH=GEP^COS  - T METRO  - 

9 ] 3 o 

G R B V = G E B ♦ 8 1 N • TMETR D ' 

9 1 4 0 

167 

v 1 = X 1 + ' ' S C - S C □ > ♦ C □ S < THETRB  1 

9 1 7 0 

v 1 = v 1 + ' : s c - s c n > ♦sin  <thetrp > 

3 1 6 Of 

IJPDRTE  OLP  VRLUES. 

9 1 7 0 

volo=vol. 

9 1 8 0 

Ph7=RR 

q i 9 0 

PEL  SF=Y  >:  7 -Y  0 8 - 

9 £ 0 0 

C R L L S PR  I N G <: F N - S F 1 J T - RE L S F - S IJ N 0 . 0 - S F 1 N S F > 

981  0 

HPRD=  1 V «:' 7 > -V  ' 81  > *18. 

988  0 

I F r HPRP-BEL.TTE -CO'S  CTHLD>  170-170,  173 

98  3 0 

1 7 0 

THEF=THFn 

984  0 

THEL.=THL  0 

9 8 7 0 

RELKE-Y  <'7':>  -Y  C8> 

986  0 

GO  TO  £70 

987  0 

178 

I P 'PELT tp . l T . 0 . 1 GO  TO  7 1 4 

9 8 8 0 

NLOOP=l 

Q P Q fj 

EPSLON=. 0001 

Q “ f|  f i 

ROOT G=l  T ♦'*  OS  CTHFl  ■' 

9 8 1 0 

POOT=ROOTG 

9880 

177 

FROOT =RDDT ♦X8+YS^S0RT  *'1  T ♦♦8-RODT ♦♦£>  - ''LT^E- 

Q 3 ! j 

D F R 0 D T = X 8 - Y 8 ♦ROOT S' 0 R T a T ♦ ♦ 8 - R □ 0 T ♦ ♦ £ > 

934  0 

DELPT=-FROOT  -PEROOT 

9 3 7 0 

root =root+ del.rt 

9 9 *9 0 

I F '•  R E S •:  P E L R T . L E.  EPS L 0 N > GO  TO  1 8 0 

9 '3  7 0 

NLOOP=NLDDP+t 

9 8 8 0 

I F -NL OOP . GT . 8 0 . > GO  TO  7 1 6 

9 3 9 !i 

GO  TO  177 

84  0 0 

1 8 0 

V ROOT=S  1? R T ■- L T ♦ ♦ 3 —ROOT ♦♦  8 

9 4 1 0 

THEF=RTRN  <YRD0T-Y2>  - LX 8 -ROOT  - > 

9 4 j1  i'i 

THEL=RTRN  - YROOT.-- ROOT 

9 4 3 o 

IF  '-ROOT,  LT.  0.  1 THEL  = 1 . 77 07963-RT RN  (ROOT  - YROOT 

944  0 

P E L K R = L T •- 18.  ♦ f C 0 O'  c THLO - C 0 S ( T M F L V:-  + P E L TIB •--  < 1 

9470 

870 

1; R L L SPR I N G < L R N , PUT  - R E L K R - S K F - 0 . 0 - F Y N E E - N Y F' 

946  0 

FX~-  ('SF +FKNE.E+  ( FC  HEST+FCBS-'  ♦COS  ' V •:'?>  > -t-FHERP^ 

947  0 

H=  1 

948  0 

GO  TO  740 

949  0 

7 0 1 j 

XSTOF-X 

9 7 0 0 

N=N+ 1 

97 1 0 

PRINT  £ 

952  0 

r-  " 

I E ( N „ E 0 . £ > P PINT  -"PR  0-  0 E N G E F'  0 U P M R R I N I N G - N 0 T 

jj  * 

953  0 

GO  TO  740 

8 7 4 1? 

7 1 0 

XSTOP=X 

O c c fi 

' .J  J v 

N=N+ 1 

8 7 E 0 

PRINT  8 

9770 

I E CN . EO . £> PR I NT  - " LOWER  CHEST  I MPRC T W I TH  PR 

9780 

GO  TO  740 

979 1"1 

718 

XS TOF-X 

D - 19 


STDPF’F 


£0 


Phi;  08  ■■■'30. •■■'8? 

9 £00  N=N+1 

96 1 0 PRINT  £ 

9680  IF  CM.  EG'.  £}  PRINT  - " UPPER  CHEST  IMPACT  WITH  DASH*  RUN  STOPPED.  " 

-M_<  GO  TO  5 9 !i 

964  0 SI  4 SSTDP=S 
96 5 0 n=N+ l 

9 6 6 1 j p P I N T £ 

96 t 0 IF  c'M.  EG.  £>  PRINT'  "X£  INPIJT  NOT  CONSISTENT  M I TH  OTHER  LEG  INPUT.  " 

9 6 8 0 GO  TO  54 0 

9690  516  YSTOP=X 

9700  N=N+1 

9710  PR  I NT  £ 

97  c!  0 I E ( N . EG . £':<  PR  I NT.  ■■  k NEE  RNGLE  SOLUT I ON  NOT  CONVERG I NG . " 

9730  GO  TO  540 

974  0 518  '■■  0 TOP= 7 

9750  N-N+l 

9760  PRINT  £ 

977 0 IE >N. EG. £} PRINT, "HERD  POTATION  > 90  DEG. - RUN  STOPPED." 

o y 9 0 9 n TO  54  n 

9 7 9 o 5 £ 0 ':■•; STDP= X 

9800  N-N+i 

9810  PRINT  £ 

988  0 I E •:  N . EG . £>  PR  I NT  - " ERG  POLL.  DEPLOYING  HIGHER  THRN  HERD*  STOP  RUN.” 

9830  540  RETURN 

9840  END 

9:85  Of 

986 OC  THIS  SUBROUTINE  CRLCULRTES  THE  STERNUM  BRGSLRP  FORCE  RS  R FUNCTION 

987 OC  OF  THE  RIRPflG  POLL  INTO  THE  STERNUM. 

9 8 8 0 S U B R OUT  I NE  STFO R ■" D R - F S T F S > 

9 3 9 0 C 0 M M 0 N ■-  S TFOR C E N S T - S TF<£*  £ 4 "• 

99  0 0 CRLL  LOOKUP  -"DPS < S TF , NS T * EG  JBS> 

9910  RETURN 

9980  END 

9 9 9 ft  f ■ 

994 OC  THIS  SUBROUTINE  CRLCULRTES  THE  CHEST  BRGSLRP  FORCE  RS  R FUNCTION 

995 OC  OF  THE  DEFLECTION  OF  THE  STERNUM  RELATIVE  TO  THE  CHEST. 

9 9 6 0 S U B R 0 U T I N E C EOF'  C D S C * FCB S > 

9 9 7 0 C 0 M M 0 N CFOPCE N C * C F (. 8 - £ 4 j 

9 9 8 0 C RLL  LOO E U P 1 D S C * CF*  N C * F C B S 

9990  RETURN 

1 0 0 0 0 END 

1 0 0 1 OC 

1 0 0c! OC  THIS  SUBROUTINE  COMPUTES  THE  RRTE  THRT  GRS  ENTERS  THE  ERG. 

1 0 0 3 0 S U B R 0 U TINE  G R S I N •"  y.  „ G I N 

1 IJ  0 4 U L 0 M M 0 N •■■  R R !;;  FLO  N F‘  G ’ G E N 1 8 ’ 8 4 '< 

1 0 0 5 0 C R l..  L L 0 0 E U F'  C Y.  - G E N < N F'  G , G I N > 

10060  RETURN 

10070  END 

1008 OC 

1 0 0 9 0 C T H I S S U F R 0 U T I N E PLACE S CERTRI N V R L. IJ E S I N M RTF' I Y FORMRT  FOR  PRINTING 
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81 


PRC 

08-'  30  ■••'8  8 

1 0 1 0 0 
lull  o 
1 0 1 8 0 
1 0 1 3 o 
1 0 1 4 0 
i o 1 s n 
1 0 ! 6 0 
1 0 1 7 0 
1 0 1 8 0 
1 0 1 9 0 
1 030  0 
1 08 1 0 
1 088  0 
1 083  0 
1 084  0 
1 085  0 
1 086  0 
1 087  0 
1 0880 
1 0 8 9 0 
1 0 3 0 0 
1 031  0 
1 038  0 
1 0 33  0 
1 034  0 
1 0350 
1 0360 
1 037  0 
1 0380 
1 0 9 9 0 
1 04  0 0 
1 041  0 
1 048  0 
1 043  0 
1 044  0 
1 045  0 
1 046  0 
1 047  0 
1 048  0 
1 049  or 
1 05  0 0C 
1051 OC 

.1  058  OC 
1 053  OC 
1 054 OC 
1 055  0 
1 0560 
1 057  0 
1 058  0 
1 059  0 

3 ij  B P D U TINE  P P T N 7 1 C X - V - D V ) 

L □MMnN-"'  OUT  1 •'  N t T 175."'  ' X4  -:  6 , 1 75  1 ? X5  6 * 175)  ’ >*,6  ••  S ? 1 <’  5 - 
C 0 M r-1 0 N .-■■■  s TFOP C E ■ " N S T 7 STF  < 8 - 8 4 ) 

C □ r-1  M □ N •••■■  C FOP C E •••••  N C 7 C F >: ' 8 - 8 4 > 

C 0 M N 0 N IT N R N D R T 3 T D P M P k C D R M P 7 C H E S T T 7 HERDT - THFCU  THLO-  DELTTB 
C 0 M M 0 N B R G DRT V C 1 7 V C 8 - R V - E 3 R - F 8 C ? FSEb  ROj  MBPG-*  DPOLLZ?  THETRDj  XI  « V 1 
COMMON  ■MERRDHT  .-r'N  - „ BINF  7 MSOCK  * DUPBZj  VCGH-  VCRV 
':!  ? SR  , v B i YL  7 YRp 

c o M M 0 N ••••  M n H T ri u P BOH,  D U P B 6 V - D F B 6 H 7 DF' B G V 7 V U P B G H 7 V P B G H 7 V 1 1 P B G V - V P B G V 7 
CGUPBH  7 bljPBV ' bPBH ' b F1  E'  V * Vo  bH 

COMMON -PRPRM  PET, THETRT7- THETRH7 7 PG 1 7 FENEE’ GFB, DPB' ROT 
U 0 M M 0 N M P R P R M .■■■' F 7 p y T E! ' E CB'Y • bUPE: 7 BUPP 7 DF‘L-7  DF'  o 7 Di-L 
C0MM0N.""MMPRPRM.-BP7  VOL  - F ETC  7 E PC  7 VSC'7  MPEb  DP  DLL  7 FF'DLL  7 XPT  7 D TOP  80 
DOUBLE  F'PEC  1 8 1 ON  Y CR> 

DIMENSION  DV C8> 

N-N+l 

DRT R F.""l  £.  "" 

IF  CM.  GT.  175:-  PETUF'N 
X 4 1 7 N = b F'  E: 

X 4 3 ’ N ■'  = V F'  B b H • ' F 

X4  -:'3 7 N>  = t;  c'VpBGH-VC GH;-  ♦COS’  -"THE T RT--  + <'VPBGV-VCGV>  ♦SIN  - THETRT".  :« ,-F 
X4  4 7 N - = t'VPBGH-Y  <‘4>  ♦18.  > ♦C03  THETRD':-  +VP BGV+8IN  >'  THFTRD  - - F 
X4  -"57  N"-  =DPBGH 
X 4 6 7 N = — D P B 
X5  1 7 N 1 =GUPB 

X5  >:  £‘7  N - =VUPBGH  "F 

X5  3 7 N>  = ,:'VUPBGH--VCGH ■-  ♦COS  -"THETRT:-  + -" VIJPBGV-VC GV :•  ♦SIN  - THETRT)  ) .-"F 

X5  -"4.7  N-  = •"  -'  VUPBGH-V  4 - ♦18.  > ♦COS  -"THETRD)  +VUP  EG '/♦SIN  - THETRD)  ) ■ F 

X5  -"5  7 N)  =DUPBGH 

X5  '"67  N"-  =-DUF‘B 

X6  1 7 N - =FCBS 

X6  ■:  8 7 N)  =FSTBS 

X6  -"3  7 N)  =-VSC-"'F 

X6  -:  4 7 N)  =-DPS 

X6  -"57  N - =-DSC 

X6 <6- N) =DTOPSO 

T -N) = X ♦ 1 00  0. 

PETUF'N 

END 

THIS  SUBROUTINE  SOLVES  THE  D I FEE  PENT  I RL  EC'URTIONS  THRT  DETERMINE 
THE  PRSSENGEF'  L I NEMRT I CS  RND  VEHICLE  MOTION.  THE.  FOURTH  ORDER  FUNGE 
KUTTR  METHOD  IS  USED  TO  STRPT  THE  I NTEGPRT I ON 7 BUT  ONCE  THE  FIRST 
FOUR  POINTS  RF'E  OBTRINED  WE  SWITCH  TO  THE  MOPE  ECDNDMICRL  FOURTH 
ORDER  RDRM3 -MOULTON  F'F'ED I CTDP-C OF'PECTOP  METHOD. 

SUBROUTINE  SOLVE  rN) 

COMMON  T I ME  ■••"STEP  7 XSTOF'  7 DELTRT , F'lNTl  7 PI  NTS 
C 0 M M 0 N .••"  G R S E l.  0 N p G 7 G E N 8 ' 8 4 ) 

D 0 U T L E P P E C I S 1 0 N Y C 3: ) 7 V T < 3; ) 

DOUBL E F'PEC  I S I ON  B87  0 7 B 1 9 7 BS5 1 

D - 21 


PHU  0 8 v-' 3 iV' HE' 

9600  n=N+1 

9610  PRINT  P 

9 6 SO  IF  CM.  EC.  S>  PRINT  , " UPPER  CHEST  IMF'RCT  WITH  DRSH*  RIJN  STOPPED.  " 

9 6 3 0 60  TO  54 0 

90  4 0 514  95  TOP =9 

9650  N-N+l 

9660  PRINT  £ 

90  TO  I P >:  N . EC . 8 > PRINT"  " PE:  INPIJT  NOT  CONSISTENT  I.IITH  OTHER  LEG  INPUT." 

90,80  GO  TO  540 

9690  516  PS TOP =9 

9700  N~N+ 1 

9710  PRINT'  2 

9 7 S 0 I F 1 N . E >0 . £ ">  PRINT."  E N F E R N G L E S 0 L U T 1 0 N N 0 T C 0 N V E R G I N G.  " 

9730  GO  TO  540 

9740  518  PS TOP =9 

9750  N-N+l 

9760  PRINT  2 

977  0 I F ON . EP . £> PP I NT , " HERD  POTRTION  > 90  DEG.  - PUN  STOPPED." 

9 7 9 0 G 0 TO  54  0 

9790  530  PSTOP=P 

9800  N-N+l 

9310  PRINT  £ 

988  0 IE -:N.  EQ.  £i  PRINT,  "ERG  POLL  DEPLOYING  HIGHER  THRN  HERD-  STOP  PUN.  " 

9830  540  RETURN 

9840  END 

835  OC 

986 OC  THIS  SUBROUTINE  ERL CULRTE S THE  STERNUM  BRGSLRP  FORCE  RS  R FUNCTION 

987 OC  OF  THE  RIRBRG  POLL  INTO  THE  STERNUM. 

9 8 8 0 "7 IJ  B R OUT  I N E ST  FOP  D P S - F S T B S > 

9 8 9 0 C 0 M M 0 N S TFOPCE N S T - S' T F C £ - £ 4 > 

9 9 0 0 C R L. L L 00 KU P >" D P S'  - S T F - N ST  * F S T B S' 

9910  RETURN 

9980  END 

9 9 '3  0 C 

994 OC  THIS  SUBROUTINE  CRLCUL.RTFS  THE  CHEST  BRGSLRP  FORCE  RS  R FUNCTION 

995 OC  OF  THE  DEFLECTION  OF  THE  STERNUM  EEL RT I VE  TO  THE  CHEST. 

9 9 6 0 S U B R 0 L l T I N E CFO  R < D S'  C - E C E S ■' 

9 9 7 0 C 0 M M 0 N C EOF' C E N C - C F < 8 - £ 4 

9 9 8 0 CRLL  LOO L U P >'  D S C -OF,  N C - E C B S •' 

9990  RETURN 

10000  END 

1 0 0 1 OC 

1 0 0 8 0 C T H I S S U B P 0 U TINE  C OMPIJTES  THE  PRTE  THRT  G R S ENTERS  T H E ERG. 

1 0 0 3 0 S U B P 0 U TINE  G R S I N X - n I N "■ 

1 0 U 4 ij  U 0 M M 0 N b R S FLO N P G - G E N (. c'  - c 4 1 

1 0 0 5 0 C RLL  LOO K U P Y - G E N - N P G - 0 I N > 

10060  RETURN 

1 U IJ  7 ij  E N D 

10080C 

1 0 0 9 0 C T H T S S U B POUT  I N E PL.RCE  S CEF'TRI  N V R L U E S I N M RTPI X FOPMRT  FOP  PRINTING 
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E'l 


PRC 

08-"90-'8i. 

1 0 1 0 0 

SUBROUTINE  PRINT! < 

X-  V-  DV’-1 

1011  0 

COMMON. 

OUT1.--N,  TC17 

5 i X 4 6 - 175  * X 5 (6,  * 175'-  - X 6 - 6 - 175’- 

1 0 1 8 0 

COMMON. 

-'STFORCE-  - NST 

- STF  US  - 84> 

1 0 1 3 0 

COMMON 

■"CFOPCE  - NC  - C 

E -S- 84' 

1 0 1 4 0 

COMMON. 

••MMRNDRT--STDRMP-  CDRMP-  CHESTT  - HERDT-  THFO-  THLD-  DELTTB 

1015  0 

COMMON- 

BRGDRT  -'VCl  - 

VC 8 ? RV  - F SR  - FSC  ? FSB  - R 0 - MpJRb  - DROL  LS  - THF  TRD- 

1 0 1 6 0 

COMMON. 

•'MBRGDRT-'DMS 

- DINF- MSOCK- DURBZ - VCGH- VCGV 

1 0 1 7 0 

!v  i ''Hi 

3 - I i X R R 

1 0 1 8 0 

COMMON. 

'MDRT  •••'  DURPGH 

- D U R B G V - D R B G H-  DR  S' G V - V U R B G H - V R BGH-  V U R B G V - V 

1 0 1 9 0 

& GUP EH 

H 1 J P p V i H P PH  i 

Gppv- VSGH 

1 0 8 0 0 

C OMMDN 

••  PRPRM  PFT.  THETRT7-  THETRHZ-  F'Gl  - FKNEE- GEB-  DEP - RD I 

1 081  0 

COMMON. 

■ MPRRRM..-  SE-  R 

STBS  - FCES - CURB- DUPB- DEC- DPS- DSC 

1 088  0 

f OMMDN. 

•'MMRRRRM.--' BP  - 

VOL  - FETC-  FPC-  VSC  - WEB  - DROLL  - RROLL-  XPT  - DTOR 

I 083  0 

DOUBLE 

PRECISION  V 

r j-j  j 

1 084  0 

D I MENS 

ION  DV  C3> 

1 085  0 

N=N  + 1 

1 086  0 

DRTR  E. 

• 1 8 . 

1 087  0 

IF - N. GT. 175- RETURN 

1 0 8 8 0 

X4  (1  - N' 

- =GRB 

1 089  0 

X4  as  - N 

- =VP  BGH  E 

.1  0 9 0 0 

H4  - "3 , n: 

. = >•  .rVRBGH-VC 

GH':>  ♦COS  -THE T RT>  + VEBGV-VCGV  - ♦SIN  ''THETRT'-  > 

1 031  0 

V 4 -:'4  - N 

, = , i:  VRBGH-V 

4>  ♦18.  > ♦COS  ( THETRD  ■-  +V  P BG '•.■'♦  SIN  THE  TRD:-  - -F 

1 038  0 

X4  -:'5  - N: 

- = DP BGH 

1 0 3 3 0 

X *4  ■*  k i N 

- =-DPB 

1 094  0 

V5  1 - N 

- =GUP B 

1 0 350 

X5  - 8-  N' 

- = VUE BGH  "E 

1 0 36  0 

xs  c9 , n: 

, = ,•  CVURBGH-V 

C G H ♦ C 0 S -•  THETRT  > + -'  V U R B G V - V C b V > ♦ S I N -THETRT 

1 037  0 

V5  -:4  - N 

• = »'  ' V U P P H — y 

4 > ♦18.  > ♦COS  - THETRD"-  +VURBGV^SIN  -THETRD')  > 

1 033  0 

X5  <5  - N' 

- =DUPBGH 

1 0 39  0 

H 5 S'  - N 

- =-DUPB 

1 0 4 0 0 

X6  >:  l , N! 

- =FC BS 

1 041  0 

X 6 -•  8 - N 

• =ESTBS 

1 048  0 

X6  C3-  N 

- =-VSC--'  E 

1 043  0 

X6  -'4 - N 

- =-DPS 

1 0 4 4 o 

X6  -:'5-  N 

- =-DSC 

1 0450 

X6  6'  - N 

- =DTDPSO 

1 046  0 
1 047  0 
1 048  0 
1 049 0C 
1 06  0 OC 
1 05! OC 
1 058  OC 
1 053 OC 
1 054  or 
1 055  0 
1 056  0 
1 057  0 
1 058  0 
1 059  0 


T =>:>  1 f 

return 

END 

THIS  SUBPniJT I NE  SOLVES'  THE  DIFFEPENTIRL  ERUPTIONS  THRT  DETERMINE 
THE  PH S' SENGEP  CINE MR TICS  RND  VEHICLE  MOTION.  THE  FOURTH  ORDER  PUNGE 
KUTTR  METHOD  I S USED  TO  STRRT  THE  I NTEGRRT I ON?  BUT  ONCE  THE  FIRST 
FOUR  POINTS:  RRE  OPT  R I NED  ME  SMIT'CH  TO  THE  MORE  ECONOMICRL  FOURTH 
ORDER  RDRMS-MOULTON  PREDICTOR-CORRECTOR  METHOD. 

SUBROUTINE  SOLVE  -TO 

COMMON  -TIME ••-STEP  - XSTOP-  DELTRT , PINT  1 - PINTS 
C 0 M M 0 N b p s E L 0 N p G-  GE N 8.  84 > 

D 0 U P L E PPE IS  I S 1 0 N V < 8 > - V T 8 > 

DOUBL  E PPEC I S I ON  B87  0 - B 1 9 - B85 1 
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PRC 

i' 

:i9..-30--83 

1 0 6 1 j i j 

D I MEN 8 1 ON  DY  CS  '<  , F •:"?  - 8>  * J < 3':-  • 

1 OR  1 0 

B85 1=851 . 

t 06?  fi 

ESS 0=870. 

l fiRpn 

E 1 9= 1 9 „ 

1 064  0 

PPI 1=0. 

1 065  0 

PP 1 8=  0 . 

1 066  0 

DFLTRT = o „ 

1 i j 6 7 ''i 

CRLL  SETUP  y. 

1 068  0 

r RL  L I PEL1 1 ' n y i PY  1 

1 0 6 8 o 

DELTRT = STEP  "4 . 

1 0 7 0 0 

CRLL  PPI  Nil  - - '9  DY:' 

1071  0 

CRLL  PPI  NTS  >8,  y,  PY  •' 

1 0 7 8 0 

IE  ' /'STOP.  hi.  , 85-  7ST0P=.  85 

i n 7 7 o f 

STRPT  OP  INTERPRTIDN  POUT INI 

1 074  or 

PUNR'E  LIJTTR  STRPT  UP 

1 1'!  7 7 0 

J -■  1 ' = 1 

1 076  o 

J ' S3  '1  =8 

1 077  0 

! r 8’  = 8 

1 078  0 

DO  7870  0=1*  N 

1 0 7 8 0 

787  0 

P I-  3-  L =DY  ' O'- 

1 'i  8 0 o 

DO  755  0 JO  = 1. , 3 

1 081  0 

D 0 7 3 i'i  0 L = 1 - N 

1 088  0 

7 8 0 o 

- 1 -!  =py  ."0:  ♦ STEP 

1 088  0 

YN=7+STFP 

1 084  0 

DO  7330  0=1 i N 

1 085  0 

7* f 1 

Y t i-'  !•••  'i  — y •:  k '■  + S ( O'  'i  8 . 

1 086  o 

CRLL  PI  FEU  *:'7N,  YT-  PY:- 

1 087  0 

DO  7360  L=1 ’ N 

1 088  0 

796  0 

8 i' 6 ''  =■"’ 1 0 '<  +8.  ♦DY  ' 0 ♦'“  TEP 

1 089  0 

DO  738  0 0 =1  * N 

1 0 8 n o 

738  it 

Y T ' O'  > = Y ( 0 ' + Y T E F‘  ♦ P Y 1 0 1 8 • 

1 08  i o 

L’RLL  P I EEU  1 ' N * YT  * pY  '< 

1 08  8 0 

PO  741  0 0 =1*  fi 

1 098  0 

741  0 

S • 0 '<  =S  i'  O'  ■'  +8 . ♦PY  ''  O'  ♦ST EP 

1 084  0 

PD  743  0 O'  = 1 * N 

1 095  o 

743  0 

YT  '"O'  > = Y 1 O'"1  +PY  •'  0 "•  ♦ " TEP 

1 fi  96  0 

Y=Y+S'TEP 

1 097  0 

L R L.  L.  PIPE  0 1 ''  K * Y T * P Y 1 

1 0980 

PO  747  0 0"  = 1 , N 

1 099  0 

747  0 

Y '"0" ''  = Y '"0  ' + '" S ’'O  ' +PY  ' 0 '■  ♦STEP'1 

1 1 0 o o 

L’RLL  PIPE  U '■  " * Y * P Y 1 

1101  0 

COTO  >'  75  0 0 * 753  O'  * 755  0 '■  * jo 

1 1 08  0 

75  0 0 

PO  7510  0=1* N 

1 1 0 80 

751  0 

6 r 3 * O'  ''  = P Y '"  0'  •' 

1 1 04  o 

GO  TO  7550 

1 1 07  0 

753  0 

PO  754  0 0=1  * N 

1 .1  06  0 

754  0 

F '"  1 * 03'  =PY  ■:  O'  ' 

1 1 07  0 

755  0 

PPI 1=7 

1 1 08  0 

PPI 8=7 

1 1 09 OC 

P P E P I IS  T 0 P - C 0 P P E C T 0 P :S  E C T 1 0 N 
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PRC 


i.i  y 0 ■•'b'c 


1 l l nor 
1 1 1 1 0 
1 1 1 £:  o 
1 1 1 ? 0 
1 1 1 4 HP 
1 1 1 *5  0 
1116  0 
1 1 17  nr 
1 1 1 8 0 

I 1 1 9 h 

I I ? n n 

1 let  oc 
1 1 88  0 
1 1830 

I 1 £4  0 
1 18  “5  0 

I I £ 6 0 r 
1 1£7  0 

1 1 £• ! j 

1 1 £9  0 
1 1 3 0 OC 
1 1 3 1 0 
1 1 3£  0 
1 133 OC 
11340 
1 1350 
1 1 36 OC 
11370 
11380 
1 1 39  0 


PREDICTDP 

7590  DO  760  0 0=0  N 

76  0 0 YT  Ck’>  = Y CK>  + S'TEF'  + ' 55.  ♦ BY  - k '-  —59.  ♦F  C J Cl  ’>  . 0 +37.  ♦F  C J <£?  •>  k 1 
H . + F l I ' ;■  *t  . I-*  1 1 L-'4  . 

S'Rp'F  DV  ••••": 

DO  764  0 0 = 1 i N 
764  0 P C J 3 > . 0 :>  =DY  -'  k 
EVRLURTF  STEP 
Y=Y+STEF' 

DFLTRT ='c  TEP  ■••'£. 

7 7 0 1 j C h L L DIFE 0 C Y , Y T •>  D Y > 

POthTF  vector  pointer 

IJT=.J  r 3 :« 

J 1- 3 = J f-p) 

J (’£'■'  - J < 1 > 

J i-  J =1  IT 
CORRECTOR 
DO  775m  0=1, N 

775!'!  Y 1 O'  j = Y '•  O'  1 + S' T E F‘ ♦ •’  9 . ♦D Y 1 0 '<  +19.  + F ( J < 1 > , O’  1 — 5 ..  ♦F  <.  J •’  c > » k 1 + F >.  J t :>  . 
6 k > > £4  . 

RDDITIDN  OF  ERROR  TERM 
DO  78  00  0 =1  - N 

78  0 0 Y CO:  > = CP£5 1 ♦Y  1 k +E 1 9 + VT  C0">  '•  •■•'  P£ 7 0 
S ECOND  EVPUJRT I ON  STEP 
L ' H L L DIFE  y c Y , Y * T'  Y ’’ 
ijhLL  UPDHTE  t Y ' Y i D r 1 
PRINTING  SECTION 
F'R  I 1 =PR  1 1 +STEF' 

RRI£=RRI£+STER 

I F CRR II  . L T . R I NT  1 GO  TO  789  0 


114  0 0 

RR 1 1 =RR 1 1 -P I NT 1 

11410 

CALL  PR  INTI  CY,  Y.  DY':- 

1 148  0 

789 0 I P CPE I8.LT.PI NT£> G0T07R8 0 

1143  0 

F'R  I £=PE  I £-R  I NT  £ 

1 144  0 

CRLL  PRINTS'  - Y-  Y - DY':- 

11450 

7 9 £ 0 I F C Y . LT.  Y S' T 0 P > G O' T 0 7 5 9 0 

1 1 46  0 

RETURN 

1 1 47  0 

END 

1 1 43 OC 
1 1 49 OC 

THIS  SUBROUTINE  COMPUTES  THE  O'NEE  RESTRAINT 

CRUSH  FORCE  RND 

1 1 5 0 OC 

THE  SERT  FRICTION  FORCE.  HYSTERESIS  EFFECTS 

MRY  BE  INCLUDED 

1151  0 

S'  U E P 0 U T I N E S'  P R I N G F , D F L ' T R , D I S T - S L 0 P F 1 < S L 0 P F £ 

, FORCE  ’ NF'TS '- 

1 158  0 

DIMENSION  F -'  £ • £4  - 

11530 

E=  0 . 

11540 

I F CD I ST.  GF . DFL.TR;-  GO  TO  834 0 

11550 

M=£ 

11560 

UT=DIST 

11570 

GO  TO  8360 

11530 

8110  M=3 

1 1 59  0 

UT=DFLTR 

D - 23 


£4 


PRC  0 8 ■■•"3  0 ■••'8  c' 


! 1 6 o o 

GO  TO  8360 

1161  0 

8 1 4 0 

F 1 , k>  =DELTA 

1 1 6cO 

F £ i k>  = FORCE 

1 1 6 Cm 

IF  OC.  GT.  3>  GOTHIC  1 8 0 

1 1 6 4 i"i 

G □ TO  8 8 4 0 1 8 8 4 0 » 8 3 0 0 '>  - k' 

11650 

8 1 8 0 

k k=k-3 

1 1 66  0 

R=1 . 

1 1 67  0 

DD  8810  L=8 - NF'TS 

1 1 68  m 

F 1 1 - L>  =F  r 1 , L +k'K':> 

1 1 6*0 

881  0 

F C8,  L>  =F  Cc'n  L+Kk  > 

1 1 7 0 0 

N F T '3  - N P T 8 — k k 

1171  0 

GO  TO  8300 

1 1 78  0 

884  0 

k k-3-k 

1 178  0 

P = 8 _ 

1 1 74  0 

DO  888  0 L.L=  1 ’ NRTS 

1 1 75  0 

L-MPT8+1-LL 

1176  0 

F ' l i L +k  k 1 =F  >'  1 • L> 

1 1 77  0 

888  0 

F 1 c.  * L + K • — F *.  i~' i [_) 

1 178  0 

npts=npts +k k 

1 1 790 

8 3 0 0 

F 1 , £ > =F  e 1 , 3 '■  -F  > 3 , 3 > • C LOPE  1 

1 1 8 0 0 

F 1 c *•  c'  •'  = H „ 

1181  0 

6 can  l>  =-SL 0PE8+F  .-1  , £> 

1 188  0 

FCln  1 = 0 „ 

1 188  0 

884  0 

M=1 

1 184  0 

U T = D I '7  T 

11850 

836  0 

DD  837  0 .1=8,  NPTS 

1 1 86  0 

837  0 

I F - F e i , J-.  . RT . IJT>  GDTD838 0 

11870 

8 9 8 n 

k=,J-l 

.1  1 88  0 

839  0 

FORCE  — 1 LIT  —F  • 1 , k 1 ♦ 1 F 1 c' , .J  1 - F 1 c , k ' 1 ' 1 F ■:  1 , _l  1 —F  '•  1 , k ' +F  '•  8 , k ■' 

1 189  0 

IF  ep.  EQ.  1 . :■  N F‘ T S' = N F* 7 3 + k k' 

1 1 9 0 0 

IF  ep.  EC.  8.  NPTS=NPTS-kk 

1 1 9 1 0 

GO  ID  C845  0 1 84 1 0 ■>814  O '1  ■>  M 

1 l 98  0 

841  0 

I F > FORCE . LE . 0 . G0T0845  0 

1 1 98  0 

I F 'HOC  e e F CS , J > -F  ca , k :>  > • -:  F e | , J > -F  e 1 ■>  k > > -8  LOPE  1 > . LT.  . 0 1 > 

1 1 9 6 0 

GOTO 8 450 

1 195  0 

GO  TO  8110 

1 1 96  0 

845  0 

RETURN 

1 1 97  0 

END 

1 1 98 0C 
1 1 99  Of 

THIS:  SUBROUTINE  PLACES  CERTAIN  VALUES  IN  MATRIX  FORMAT  FOR  PRINTING 

1 8 0 0 0 

8: U FP OUT  INF  p p I N T 8 X , V ■>  D V > 

1 8 0 1 0 

L 0 M M 0 DU  TV  fS  i T * | i ■ S • n t‘\  m *.  F«  i t i‘  5 1 ^ t F»  * 1 75  •'  * •■311'  '*  F*  n 1 7’,5.1  •»  s'  '•  Ft  •»  1 -1  q v-s  »■  •’  F.  i IT*1 

1 £ 08  0 

!j.: X o ' F.  ^ | X S > * !•••!  'H  • F>  n 1 7 5 ’*  * !•••!  1 i.i  * F.  i 1 , # 15  1 i 11  v F.  i ] 5 .-1 

1 8 0 8 0 

DOUBLE  PRECISION  V - 8 ' 

18  040 

DIMENSION  DY ' 8' 

18  050 

C OMMON  - MANDAT  - ZL  - ZT , ZS , ZHi-RT-  RN , RH  - PTOPR  - X8Z  - Y8Z - WB  - B 0 , MH , L T , LF 

1 8 06  0 

COMMON  MMANDAT  -'STDAMP,  CDAMP,  CHESTT » HFADT  - THFO-  THL 0->  DFLTTB 

1 8 07  0 

r 

COMMON  ■■■BAG  DAT  VC  1 , VC8  > AV  ■>  F S A ■>  FSC , FSB , A 0 , ME  AG  * DROLL  Z , THE  T AD , k ONF  . YON 

t: 

1 £ 08  0 

COMMON  - MB  AGDAT  DM  ~ , DTNFn  MS'CCk',  DUPBZ-  VCGH  ■>  VCGV 

1 8 09  0 

’!• , S' A , SB » SC  i SAP 

D - 24 
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1 £ 1 0 0 
1211  0 
1 £ 1 £ 0 
. 1 £ 1 3 0 

1 £ 1 4 0 
1 £ 1 5 
1 £ 1 6 0 
1 £ 1 7 0 
1 £ 1 ft  ij 
1 £ 1 '7  0 
1 ££  0 0 
1 ££  1 0 
1 £££  0 
i p p n 
1 ££4  0 


COMMON  PRPRM,-PFT>  THETRTZ,  THFTRH7,  PR1  , FKNEE3  GPP  < DF'B , ROI 

COMMON-  MPEPEM.--  3F,  FS TPS  - FCB3"  RIJPP , DUPEN  DPC  - DP  3 - DSC  * GST 

C □ M M 0 N ••••' M M P R P R M ••••' BF‘«  VOL  * FFTC>  FF'O  VSf  , MFp , DROLL ? PPDLL  , LPT,  DT0P3 0 

C □ M M □ N M F'  P I N T F F T H , FPM , THEE  - THEL  BPH-  VO  I - liJUPP-  ROOT 

COMMON  MDRT  • DUPPRH-  DUPBG'A  D PpRN  - DPBRV,  VUPBGH,  VPFGH , VIJPpRV,  VPpRV 

'CRiJPpR , R ij P p v , R P B H - R P P V , V 3 R H 

COMMON  - H I C TH I C < 1 75>  ’ HPRS  <.  1 75-'  - CPG3  ■'  1 75"' 

C □ M M 0 N T I M E 3 T E P , V £ TOP,  BEL.TRT , PINT! , PI  N T £ 

DRTR  F •••■  1 £ > , R :•  £ . £ •••'  * D 57 . £ 9 5 7 y U 
N = N+  t 

IF  f'N.  RT  , t 75'  PE  THEN 
C H=C  0 7 1 V 1 6 ’■* 

CT= C 0 S 1 V r 7 .< 

S H=3  I N V '■  6 > > 

ST=£IN  T V- 


1 ££5  0 

T ' N'i  ='•••♦  i fnjO. 

1 ££6  0 

v 0 <"  1 , N > =D  V <:  4 j R 

1 ££7  0 

7 0 ( £ , N 1 - V ( 4 > 

1 £££0 

7 ij  f 3 , N j = V 1£  > ♦ E 

1 ?PR  0 

7 0 >'  4 , N 1 = P'  E’ 

1 p p f!  fi 

7 0 rO,  r-f*  =FFTC 

1 £3 1 0 

• . U 1 *z-  *»  N 1 ~ f-  P‘ ! 

1 £300 

v l >:  l , M = c y f 5 '>  - V >' 

3 ' • ♦F 

t £33  o 

7 1 1 £ , N 1 = V 1 1 '' 

1 £34  0 

71  C 3 " N =3F 

1 £35  0 

LI  f 4 , N':>  =F  KNEE  ♦CO 

3 . THEE:'  •••- 

C.  - 

1 R n 

71  5 - N •'  =7  HP  F*P 

1 £37  0 

7 1 1-6,  N , =THEL  ♦D 

1 £30  0 

.y-t  ci  1 1 n N 1 — * V >.  *5  j 4- F' T ♦ 1 N J N V V I 

7>  -3 IN  ■:  THET RTZ>  > > ♦E 

1 £ 3 9 0 

7 £ ( £ , N 1 — V 7 j ♦ 7i 

1 £4  on 

7 £ '•  3 , N ''  ~ V 3 ’■  ♦ D 

1 £ 4 ! 0 

7 El  '■  4 , N =DY  63  j ♦D 

1 £ 4 £ 0 

7 c ''  5 , N •'  ~ 7 c!  '■  1 , N > — 

Y <8> ♦F 

1 £ 4 3 ij 

7 c! 1 6 ’ N 1 ~ 1 P T ♦ V ' 3 1 

♦CT +Y  1 1 •' 

-y  '■  4 i i 

1 £44  0 

73  '■  1 , N ''  =F+  6y‘  >■  5 j +E'N+  i:'3 1 N < 

Y ' 7':.  -3  IN  'THE TRTZ ■'  > +PH 

1 £45  0 

73 IN  •' THETRHZj  > > 

1 £46  0 

73  C£j  N>  =Y  '''6>  ♦D 

1 £47  0 

7 3 C3,  M j =y  63 ) -*-ri 

1 £40  0 

7 3 i 4 , N j = D V £ j ♦ D 

1 £4Q  n 

.j  1 ^ i N 1 ~ f 1 i N j — 

Y 3 j ♦F 

1 £5  0 0 

7 ';!  i'  6 , N 'i  =:  i V t E'  j — Y 1 

7>  > *0 

1 £51  0 

77  r i , m > — pF'H 

1£5£0 

77  ■:'£,  N :•  =VOL 

1 £530 

77  «:'3  N j = P R 1 

1 £54  0 

77 <4 , N 1 =PFTH 

1 £550 

77  >:  5,  N>  = F P H 

1 £56  0 

77  -'  6 - N > =V0 1 

1 £57  0 

7 '3  ''  1 , N j = '•  D Y ( 1 ♦ C 

T+PT^DY  '• 

3 1 :•  •■■■' 3£ . £ 

1 £58  0 

78  '-£,  N ' = f-DV  < 1 j ♦ 

7 T + E‘  T ♦ Y 

3 j ♦ Y 63  > > ..■•'30  . £ 

1 £59  0 

7 '3  1 ’ ■! , N 1 = 1 E'  H ♦ Ci  Y '■  £ 

j +DY  ''  1 > ♦ 

L H + P Cl  ♦ (.  Cl  Y < .!•  ♦ '•  C H ♦ C T + 3 H • 
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1 £6  0 0 

6 ♦V  03  > ♦ OSH+C  T-CN  + S T 

■- - 3£.  £ 

t £ 6 1 0 

ft  y ( 4 * N ~ 1 F H ♦ V * c.  V 1 

£"'  +RN^  OY  06  •-  ♦V  0 6"'  * '■  S H ♦ S T +CH+CT  •-  — DY  3? 

1 £6 9 0 

g 0 £ H ♦ C T - C H ♦ S T '*•  ">  - 11  V 

' 1 -i  ♦SH':-  ' 5 £ - £ 

1 £6  6 0 

69  ••  1 - N>  -COT 

I £64  0 

69  •:'£-.  N =EROL.  L. 

13650 

69  .'3 , -i  =;,:ohe 

1 £66  0 

£ 9 i 4 ' N 1 — V D N E 

1 £67  0 

ft  U * C|  *>  N •=  i » i F‘  p! 

1 £66  o 

C F.  n f '-j  f ~ 1 1 j ! J F'  p- 

1 f j 

rp 03  ON"'  = £PRT  f v '3  i- 1 , N 

■>  ♦ ♦£+ '•  •" -6  i'  £ , N •'  ♦ ♦ £ "' 

1 £ 7 0 0 

T H I C 1 N "'  = 0 N ♦P I N 7 £ ♦ ♦ 0 . 6 

1 £ 7 1 0 

HP 0 S ON  •'  CRT  ••  66  0 6 i N 

♦♦£+£6  04.,  N>  ♦ ♦£"■ 

1 £7£  0 

6 1 0 >•  1 1 N '•  =f  R6S  ON  •' 

1 £730 

ft  | M ».  c..  * N 1 = H F’  h *j’  * N .* 

1 £74  0 

1 U 1 * N 1 = U . T 

1 £7*50 

£1  0 04"  N"'  =EKNEE+COS  o 

THEE"-  ■•-£. 

i p 7 6 0 

6 1 0 f 5 ■ N ">  = i"  0 V P BCR -VC 

EH  1 ♦COS  0 THETRT + OVPE:RV~VCPV>  ♦S  I N ■ THETRT 

1 £ 7 7 0 

£1  0 06' N 1 =PG1 

1 £76  0 

£ j 1 0 1 , N "'  =V  • 7 • ♦B 

1 7'  7 Q 0 

£ 1 1 ''  £ ' N ■'  =V  ''  b ■'  ♦ O' 

1 £6'  0 0 

£1 1 06,  N > = JHFF ♦ D 

.1  £8 1 0 

61  i 04'  N ">  =-VSC  F 

1 £ 8 £ 0 

61  i 05,  N ~ <:  V 05>  -V  0 6". 

"•  ♦F 

1 £83  0 

ft  | 1 v F,  i N i = — pi  F*  i _ 

1 £64  0 

PFTUPN 

1 £86  0 

END 

1 £86 0C 

1 £87  0 

SUBPOUT I NF  UPBRTF  06 

f pi  V 1 

13880 

DOUBLE  PRECISION  yo 

J-!  j 

1 £ 6 6 0 

DIMENSION  DY  06  "- 

1 £6  0 0 

COMMON  -KNF EPF ST. -NL R 

, 3 p p , pi  it  , pel  0 P , k PN  o £ , £4> 

1.  £91  0 

C 0 M M 0 N S ENTER  I r N s F 

•»  UN  i PUT  * P F l..  .'sF  * ft F N i‘  j;-’  * y *4ft 

1 £9£  0 

3 EUT =P,ELSF 

1 £ 9 3 0 

RUT =RELKR 

1 £94  0 

RETURN 

1 £950 

END 
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Appendix  E 


PAC  Sample  Run 


. . , 


PAC  SAMPLE  RUN  - FULL  LIST  OF  OUTPUT 


READY 

RIJN 


PRC 


16:  04EDT  OR  -'  08 ••'8c 


INPUT  FILE  NRME'"'TRCST7P 


ENTER  1 IF  YOU  WANT  FULL  LIST  OF  OUTPUT : ENTER  £ IF  YOU  WANT  ABBREVIATED  LIST. ? 1 


INPUT  VRLUES  — INPUT  UNITS  < MSEC.  MPH.  DEGREES.  INCHES.  LES.  FT-LBS.  GyS> 
INITIAL  VELOCITY:  .34.0 

INITIAL  HERD  ANGLE:  -9.50 


INITIAL  TORSO 

ANGLE: 

-£7.  5 

MLEG 

M TORSO 

IX 

LU 

h- 

•y~ 

MHEAD 

RT 

EN 

RH 

RTOPH 

7 1 . 0 

58.  4 

£.50 

11.4 

14.  0 

£ 0 . 5 

4.75 

C O ■ f 

NPTS  NECK 

NPTS  KE 

NPTS  VEH 

NPTS  SEAT 

NPTS  GAS 

SL.ST 

SL . KE 

Z\ 

o 

18 

8 

14  0. 

500E+04 

0. £4  0E+  04 

NETS  STEP 

NPTS  CHS  T 

5 

cr 

,_i 

GAS  FLOW  TIME 

0 . 

15.  0 

1 8 . 0 

£ £ • U 

£5.  0 

30.  0 

35.  0 

40.  0 

45.  0 

50.  0 

GAS  FLOW  TIME 

55.  0 

60.  0 

75.5 

1 0 0 . 

GAS.  FLOW  - LB--' 

SEC 

0. 

0. 

3 . 4 0 

3 . 9 0 

7.54 

8.  09 

7 . 5 0 

cr  ~i 

_•  ■ i 

4.  0 0 

£.54 

GAS  FLOW  - LE-" 

SEC 

1 . 49 

0. 8£0 

0. 

0. 

r rOT  f"  r*.  t -•  tthp 


E 


1 


•j  cn  i 


1 1 j ■ 1 


DISPLRCEMENT 


1 4 . 0 


r-5  0.  0 

0. 

SERT 

FRICTIDN 

FDRCE 

0. 

0. 

NECK 

RNGLE 
-81 . 0 

1 8 . 0 

NECK 

TOR CUE  - 

FT-LBS 

117. 

0 . 

VEH. 

PULSE  - 

TIME 

1 . 0 0 

LBS 

35  0.  350. 

90.  0 
-87.  0 


15.0 


10.  0 


0 . 

7.  0 0 

1 7 . 0 

8 1 . 0 

86.  0 

33.  0 

45.  0 

58.  0 

54  . O 

58.  0 

VEH.  PULSE  - 

TIME 

65.  0 

70.  0 

76.  0 

80.  0 

9 0.  0 

93.  0 

1 07. 

150. 

VEH.  PULSE  - 

DECELERRT I DN 

0. 

0 . 

1 3 . 3 

18.  7 

15.  6 

10.7 

1 8 . 3 

8 9 . 8 

3 0 . 7 

89.8 

VEH.  PULSE  - 

DECELERRT I DN 

3 1 . 3 

L-ji  a 

17.8 

19.5 

11.7 

18.6 

0 . 

0. 

KNEE  DISPLRC 

EMENT 

0. 

8.  0 0 

8 . 5 0 

8.75 

3.  0 0 

3 . 5 0 

4.50 

15.0 

KNEE  FDRCE  - 

LES 

0. 

0. 

8 0 ij . 

7 Li  0 . 

0. 180E+04 

0 . 83  OE  + 04 

0. 850E+04 

0. 870E+04 

S TER  N U M P I S P L.  R C E M E N T 

-5  0.  0 

0. 

0. 85  0 

1 . 0 0 

1 0 . 0 

STERNUM  FDRC 

E - LES 

0. 

0. 

4 0 0 . 0 . 

1 60E  + 04 

0. S50E+05 

CHEST  P I SF'LRCEMENT 

-11.8 

-1 . 85 

0. 

1 . 85 

11.8 

CHEST  FDRCE 

- LES 

— ij  . 465E+04 

-15  0. 

0. 

1 5 0 . 

0. 465E+04 

RTMDP 

PGZ 

GTZ 

IJ 

PN1 

PN8 

PN3 

14.  7 

0.  0.  16 

6E+  04 

668 . 

1 . 4 ij 

1 . 4 0 

1 . 4 0 

VC1 

VC  8 

RV 

FSR 

FSE 

FSC 

x 1 

V 1 

0 . 7 0 0 

0 . 7 0 0 

4.  0 0 

1 3 . 0 

8.  0 0 

11.0 

13.  4 

84.  5 

RO 

THETRD 

FRBI.IGT 

ST DR MR  CDRMP 

DMS 

DINF 

MS  DC K 

9 . 5 0 

0. 

8 . 4 0 

0. 

3.  0 0 

6.  0 0 

4.  0 0 

1 3 . 0 

KIN 

DRDLLZ 

r\  C.'  L_ 

Y c Z 

1,1  B 

LF 

DCN 

6 . 1 0 

5.  0 0 

31.7 

y „ yy 

18.  4 

16.7 

ij . 

THFD 

THLD 

XSTDP 

STEP 

P I NT  1 

PINTS 

1 9 . 0 

44.  0 

0.  150  0. 

1 00E-08 

0. 500E-08 

0. 500E-08 

E 


2 


T I ME 

‘"'EH  G-'S 

VEH  VEL 

VEH 

DISP 

CHEST  BP 

CM  ft  FE1PCE 

C PR  FDPCE 

'Me"' 

( G **'  S •' 

''FPS''' 

CINE 

HES"' 

I NC  HES> 

'"LBS"' 

(LBS  -• 

0. 

0. 

43.  97 

r 

. 

0. 

Ci . 

0. 

5.  0 0 

0 . 

4 

9.  87 

Cl 

99 

0. 

0 . 

Ci . 

1 0 „ 0 0 

— 3 . 9 3 

4 

3.  68 

c 

. 98 

0„ 

0. 

0. 

15.  0 0 

-1 0. 64 

4 

3 „ 5 0 

c 

. 93 

0. 

o. 

0 . 

8 0.  0 0 

-13. 35 

46. 46 

1 1 

■ >’  O 

Ci  „ 

0. 

Ci . 

35 . 00 

-15. 03 

44.  37 

14 

.51 

o„ 

0. 

fl  B 

3 0.  0 0 

-13. 8 0 

41.9  4 

•17 

. 09 

0 . 

0. 

0. 

35.  0 0 

-11.97 

4 

0.  08 

18 

. 55 

0. 

0. 

Ci . 

4 0.  00 

-15.13 

■“ 

7 „ 9 0 

31 

09 

0 ,, 

0. 

0. 

45.  0 0 

-18. 3 0 

5.31 

34 

„ 09 

0 , 

0 , 

0. 

5 0.  0 0 

-36.51 

1 . 6 0 

36 

. 1 Cl 

1 „ 73 

6 a 8 8 

117. 53 

55.  0 0 

-30. 47 

C 

6. 85 

Cl.  r 

„ 86 

3 6 3 

35 . 89 

38  0 „ 98 

6 0.  0 0 

C.  f*  u -Z-  C 

*“ 

3 . 1 0 

39 

30 

4.  74 

118. 67 

7 0 9 . 8 Q 

65.  0 0 

-31. 3 0 

1 

8 . 1 8 

3 f 

. 53 

5.  96 

Li  I'  ' H a Cl  C 

1 1 3 3 . 6 0 

7 0.  0 0 

C Z‘  n c!  U 

1 

4 . 6 0 

31 

cr  1 

a L 

7 „ 3 0 

398. 04 

] 569. 1 5 

75.  0 0 

-18.7 0 

1 1 . 33 

iZ 

n lZ  C- 

8 „ 4 0 

5W  U.  *4 

1958. 47 

30.  0 0 

-19.50 

8.34 

a O l’ 

9.  47 

753.43 

Z 'Ll  "7  1 

U.L..  O a < 1 

35.  0 0 

-15.6 0 

5.41 

■“/ 

w lZ  l‘* 

1 0.  33 

y 9 7 a 6 5 

3459.84 

90.  0 0 

-11.70 

3 . 3 1 

c •-« 

n -.Z’ 

1 0,  38 

Q48, 33 

3 469. 6 6 

95.  0 0 

- 1 0 . 8 0 

1 . 39 

“i  •“ 

. E'E’ 

1 1 . 06 

886. 55 

6c69 . c 5 

1 0 0 . 0 0 

-6 . 3 0 

0.  09 

. 69 

1 0.  84 

737. 43 

1898. 16 

1 05. 00 

-1.8  0 

0.  74 

“i  •“ 

a 67 

1 0.  38 

517. 78 

1485. 70 

110. 0 0 

0 . 

0 „ 8 0 

“i  •“ 

. 6 c 

9.44 

307. 13 

930. 54 

1 1 5 „ 0 0 

0 . 

0 „ 3 0 

•“i  ■“ 

cr  —« 

n J 1 

8.  43 

131 . 59 

441. 83 

1 c'  0 . 0 0 

0. 

0 „ 3 0 

cr 

7.34 

1 0. 86 

41 . 86 

135..  00 

0. 

0 . 8 0 

•“i  T 

. 47 

6 . 35 

0. 

0 . 

1 3 0 . 0 0 

0. 

0 . 8 0 

“i  ~ 

. 43 

5.18 

0. 

0 . 

1 35. 00 

0. 

0 . 8 0 

*Z 

u .1*  C' 

4.15 

C' . 

0 . 

1 4 0 . 0 0 

0. 

0 . 8 1? 

o .r'  z' 

3 „ 1 3 

0 . 

0. 

145. 00 

0. 

0 . 8 0 

u e'e' 

3.  09 

0. 

Ci . 

15  0.  0 0 

0. 

•j . 8 0 

~i 

a Cl  .1 

1 . 05 

0. 

0 . 

- 3 


E 


I 


TIME 

H-P  E.D. 

H-P  *7 

EL 

CM3 

^1  „ 

CINCHES:;- 

CFPS 

> 

0. 

0. 

49. 

o r 

5. 

0 0 

0.  0 0 

49. 

c-  f 

1 0 . 

0 0 

0.  00 

49. 

i”i  i~i 

15. 

0 0 

0.  05 

49. 

o o 

£ 0 . 

0 0 

0 . 1 9 

49. 

c«c 

£5. 

0 0 

0. 45 

49. 

64 

3 0 . 

0 0 

0. 83 

49 . 

£6 

35. 

0 0 

1 . 31 

43. 

64 

40. 

0 0 

1 . 37 

48. 

0 0 

45. 

0 0 

£.  53 

47. 

c r' 

5 0. 

0 0 

j m _T-  -Z> 

45. 

98 

55. 

0 0 

4.  £1 

41  . 

36 

60. 

0 0 ' 

5.  10 

36 . 

!-!  hj» 

65. 

0 0 

5.9  0 

3 0 . 

£6 

70. 

0 0 

6 . 53 

£ 3 . 

65 

75. 

0 0 

6 . 9 7 

1 6 . 

CT  i“. 
._l  »“i 

8 0 . 

0 0 

7.16 

9. 

15 

C.ET 

0 0 

7.  07 

1 . 

r’  i‘* 

9 0 . 

0 0 

6 . 73 

-4. 

95. 

0 0 

6 . £ U 

- a 

6 0 

1 0 0 . 

0 0 

5.59 

- 1 0 . 

£9 

105. 

0 0 

4.93 

- 1 0 . 

91 

110. 

0 0 

4 . 36 

-11 . 

JL 

115. 

0 0 

-11. 

4 0 

1 £ 0 . 

0 0 

3.  09 

-11 . 

•j  o 

1 £5. 

0 0 

£ . 4 6 

-11. 

■Z1  C. 

1 3 0 . 

0 0 

1 . 83 

-1 1 . 

39 

135. 

0 0 

1.19 

-11. 

49 

14  0. 

0 0 

0.  54 

-11. 

61 

145. 

0 0 

-0.11 

-11. 

r*  -Z* 

15  0. 

0 0 

~~  0 ■ r r* 

-11 . 

y 5 

FP. 

FEM  FORCE 

FEM  RNG 

TIB  HNG 

S .’i 

CLBS:;* 

CBEG> 

cb  eg;- 

0. 

19.  00 

44.  0 0 

. 03 

0. 

19.  00 

44.  0 0 

. 0£ 

0. 

19.  0 0 

44.  0 0 

B OC 

0. 

19.  00 
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451 . 88 

CZ  C ■- 

— ...!  » -1 

- 0 . 88 

- ij . ij  1 

8 1 . 78 

50.  no 

883,.  78 

6 6 „ 1 0 

-4.  8 0 

- 0 „ 04 

-0.  43 

p 0 . 8 5 

5*5.  no 

0. 

0 . 

Ij. 

o „ 

0. 

1 Q . f.  0 

6 0.  0 0 

0 .. 

0 . 

Ij. 

1 j n 

0 . 

1 3 . 8 1 j 

65.  0 0 

Ij . 

0. 

0, 

o 

ij . 

16.6  3 

7 n „ 00 

'j . 

o . 

0. 

o . 

0 , 

15.  05 

75.  0 0 

o . 

Ij. 

ij  „ 

o „ 

0 . 

1 3.  5 ~ 

80.  00 

o „ 

o. 

0. 

ij . 

o . 

1 8 . 3 0 

85.  no 

Ij . 

0 . 

0. 

0 . 

Ij . 

1 1.31 

9 0 , 0 ij 

0 . 

0. 

0. 

0 , 

Ci , 

1 Li  • 6 8 

9S.  on 

0„ 

Ij  „ 

Ij  „ 

Ij . 

ij . 

1 0.  4 5 

1 00. 00 

o. 

Ij . 

o. 

0 . 

ij . 

1 0.  P4 

1 i j 5 . 0 o 

0. 

Ij . 

0. 

Ij . 

0. 

11.8  0 

110.  0 0 

tj„ 

Ij  „ 

0 . 

Ij . 

0. 

1 8 „ 0 7 

115.0  ij 

Ij  „ 

Ij , 

o . 

1 j u 

Ij . 

17.14 

1 8 0 . 0 0 

0. 

Ij  „ 

0 „ 

0 „ 

Ij . 

1 4. 34 

185. 00 

f"; 

l'l 

- 

- 

1 5 , 5 7 

l_!  „ 

' u 

. 

l_l  „ 

l.l  „ 

1 3 0 „ 0 0 

0. 

0. 

0. 

0 . 

l'l  „ 

16. 79 

135. 00 

0 „ 

o „ 

0. 

1 j . 

Ij. 

1 3 . ij  1 

.1  4 0 . 0 0 

Ij . 

0 . 

Ij. 

0 . 

0. 

1 7 . p p 

145. 00 

0. 

0. 

Ij. 

Ij . 

0. 

8 0.  43 

1 5 0 ..  0 0 

o „ 

0 „ 

0. 

o . 

0. 

8 1 . 6 1 

E 
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TIME 

HERD  BP. 

ERG  VEIL. 

ERG  PRESS. 

HI.I.-R  FORCE 

HP  FORCE 

I NT.  VOL 

c'MSV  . 

< I MCHES> 

<CU. IN. > 

<PSIG> 

<LBS> 

<LBS> 

CCU .IN. > 

Cl 

0. 

4 1 £ . £3 

0. 

0. 

0. 

0. 

5.  0 0 

0. 

4 1 £. £3 

0. 

0. 

0. 

0. 

1 0 . 0 0 

0. 

4 1 £ . 4 0 

0. 

0. 

0. 

0. 

15.  0 0 

0. 

413. 07 

-0.  03 

0. 

0. 

0. 

£0.  00 

0 . 

554. 3£ 

1 1 . 99 

0. 

0. 

0. 

£5.  0 0 

0. 

1719.61 

£ . 3 1 

0. 

0. 

0. 

3 0.  00 

0. 

3 097 . £3 

5.  34 

0. 

0. 

0. 

35.  00 

0. 

5£90. 75 

1 . 73 

0. 

o. 

0. 

4 0.  00 

0. 

6364 . 94 

1.40 

0. 

0. 

0. 

45.  00 

0. 

3149. 05 

0 . 9 0 

0. 

0. 

0. 

50.  00 

0. 

9 063 . £6 

0.  49 

0. 

0. 

£31 . 1 £ 

55.  00 

0.  65 

93 1 5 . 43 

0.  9£ 

1 .£1 

15.  £3 

843. 00 

6 0.  00 

1.57 

3935. 70 

c.  ■ ii  r - 

13. 35 

69.  15 

1 £36 . 1 3 

65.  00 

£ . 6 9 

3434. 60 

3 . 6 1 

49.  07 

147. 93 

1 698 . £7 

70.  0 0 

3. 93 

3006. 1 0 

4.  94 

115.45 

£38. 56 

£187. 07 

75.  00 

5.31 

7539. 76 

6.17 

£03. 0£ 

3 1 0 . £6 

£ 6 6 £.91 

30.  00 

6.  6£ 

7 1 37 . 53 

7.  14 

£9£ . 98 

359. 1 0 

3073. 68 

85.  0 0 

7.  71 

6335 . 65 

f’  a T* 

370. 1 1 

389. 06 

3383. 56 

90.  0 0 

8 . 54 

6664. 39 

7.74 

410. 49 

339. 63 

356 1 . 0 0 

95.  0 0 

9.  0 0 

6636 . 95 

7.  09 

396 . 5 0 

356. 9 £ 

er  q j“j  p 

1 0 0 , 0 0 

8.  95 

6739. 05 

5.  9£ 

3£9, 17 

£97. 93 

3495 . 63 

105. 00 

3.  35 

6949. 04 

4.45 

£3  0 . 58 

-?  CT 

1 1 ■ 1 

c o r ■ i 

110.  0 0 

f*  o o r’ 

7£45 . .43 

c!  ■ yy 

1 3 1.70 

144. 33 

£991 . 09 

115.0  0 

6.17 

7 6 U 3 . 6 £ 

1 . 39 

ET  “I  2 

69.  93 

£63£. 90 

1 £ 0 . 0 0 

4 . 96 

7991.9  0 

0 . 1 3 

4 . 10 

6. 69 

££44. 6 £ 

1£5. 00 

O C;  ;I| 

8383.11 

0. 

0. 

0. 

1853. 4£ 

1 3 0 . 0 0 

£ . 95 

8767. 41 

0. 

0. 

0. 

1469. 1 £ 

135. 00 

£ . 1 £ 

9140. 45 

0. 

0. 

0. 

1 096. 03 

14  0.  00 

1 . 33 

9497. 50 

0. 

0. 

0. 

739. 0£ 

1 45. 00 

0. 47 

9831.51 

0. 

0. 

0. 

405. 0£ 

1 5 0 . 0 0 

0. 

1 0 094 . 87 

0. 

0. 

0. 

141.65 

E - 10 


T I ME 

CHEST  HP 

CHEST  ST 

HERD  HP 

HERO  ST 

f'Mv'i 

1'  G 7 » 

•:  G"  S ’• 

('  G S '' 

G 3 'i 

0. 

-0.  00 

-0.  0 0 

-0.  0 0 

0.  00 

5,  no 

0.  22 

0 „ 1 0 

0.  48 

0.  08 

1 0 „ 0 0 

0 . 3 1 

0.  08 

0.  38 

-0.  05 

1 5 „ 0 0 

0 . 3 1 

0.  0 0 

o . 3 8 

- 0 . 1 4 

80.  00 

0.  £7 

-0.  35 

0.  45 

-0.  43 

85.  0 0 

0.  19 

— 0.  73 

0.  48 

-0.  98 

30.  0 0 

0.  OH 

-1 . 43 

. 0.  58 

- 1 . 78 

? 5 . 0 0 

0.  08 

. U f' 

0.  36 

-8.  15 

40.  00 

0.  1 0 

-1 . 68 

0 . 5 8 

-3.  1 0 

4 5 . 0 0 

-4.36 

-8.  08 

-0.  36 

-0.71 

50.  00 

-8.  38 

-6.81 

4 « c b 

-4.88 

5=;.  00 

— 6 . 17 

-13.  1 3 

- o „ 88 

-13. 3 1 

b 0 . 0 0 

-13. 56 

- 1 8 . 7 0 

-5.  08 

— q a 79 

65.  0 0 

-81 . 46 

-11. 08 

-13. 54 

-4.81 

70.  0 0 

-£H. 55 

-8,  6 0 

-87, 1 8 

l l a b or! 

75.  0 0 

-37. 71 

-5.71 

- 4 0 . 4 8 

C-  “| 

n 1 — 

30.  00 

-45. 07 

Cl  n i 

-51 „ 81 

R f 

05  „ 0 0 

-48  „ 8 7 

0 . 6’  7 

— 6-  0 ..  7 6’ 

9 „ 9 b 

90.  00 

-51 . 08 

3 . 1 6' 

— b 3 m b 9 

7.  34 

95 . 0 0 

-4H. 38 

1 . 6 0 

-53. 80 

1 . 6 

t 0 0 . 0 0 

-41 . 91 

0.  33 

-41. 84 

-3.  43 

105. 00 

-88. 91 

0 . 11 

-36. 00 

-i.o  o 

no.  o o 

— 8 0 . c r* 

• ! n l 9 

-16. 38 

1.31 

115.  o 0 

-10.1 3 

1 . 5 0 

-4,  79 

5 . 1 3 

1 8 0 . 0 0 

-1.41 

1 . 95 

1 . 98 

5.  86 

1 £ 5 . 00 

0.  33 

1.91 

0.  35 

5 . 3 1 

1 3 0 . 0 0 

0.  93 

1 . 78 

“K  • :« b 

4.34 

135. 00 

1 . 33 

1.61 

-8.  84 

3 . 3 1 

14  0.  0 0 

1 . 53 

1 . 43 

-4.  31 

1.40 

1 45. 00 

1 . 6 c 

1 u c'6 

-5.4  0 

o„  88 

1 5 0 . 0 0 

1 . 66 

1.11 

- 6' . 1 3 

-0.  45 

E 
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TIME 

STM  RCC 

POLL 

PhD 

7iS  B CTP 

VC  B 

TP 

1,1  pp 

l,H  iPE 

i-  Mr 

'? 

( 9 

**» 

f INCHES 7 

'."INCHES::1 

■"INC 

HE  S ' 

•:l.BS:> 

. CL  PS'  J 

ij „ 

0. 

ij  fi 

18. 

17 

17.  35 

84. 

c 

0 

ij.  8P 

0. 88 

cr 

0 0 

cr 

9 0 

1 . 

17 

17.  85 

84. 

cr 

.j 

fi 

0.  35 

'i.  88 

1 0. 

Ij  I j 

"7 

51 

1 8 . 

17 

17.  35 

84  „ 

cr 

0 

0.  85 

■ J « c!  ci. 

I 0 „ 

Ij  Ij 

-7. 

6 u 

1 8. 

17 

17.  35 

84. 

cr 

fl 

ij . 35 

0.  88 

c‘  ij  .= 

ft  f j 

7 

♦7-  7 

18. 

i 7 

1 » H »'  .7 

84. 

5 

fi 

0. 88 

0.  88 

~J  C. 

0 ij 

"7 

41 

1 8 . 

17 

1 P n y P 

84. 

cr. 

0 

|J . 6 8 

! I . 19 

' j 'j  „ 

Ij  I j 

— fr*  n 

96 

A i-i 

.1  O c 

1 5 

81 . 44 

84,. 

cr. 

ij 

0.  87 

‘j . 1 6 

90,, 

Ij  Ij 

— K~, 

54 

18. 

1 8 

83.  94 

84. 

cr 

. j 

ij 

0.  86 

i j . i 3 

4 I j . 

0 0 

~ A „ 

cr  a 

L 

18. 

07 

84.  4C 

84. 

cr 

ij 

0.  88 

' j ,,  1 1 

45  „ 

1 j Ij 

--  1 fl  fi . 

41 

1 8 . 

i j S' 

cr  p cr 

84. 

c 

ij 

0. 8 9 

o „ 1 0 

=;!'!_ 

fi  fi 

A , ~'z  a 

ij? 

1 8 „ 

1 ij 

cr  q ~i 

C'Ja  - 7 

84. 

cr 

o 

ij  . 3 3 

0.  09 

cr  cr 

1 j Ij 

1 . 

7 ij 

A ~ 1 

1 *'  n 

g 0 

8 P. . 37 

84. 

cr 

Ij 

0„  39 

i j . i j 8 

6 0 . 

Ij  Ij 

— A n 

7 A 

17. 

cr 

86. 38 

84. 

cr; 

o 

0 . 44 

0.  08 

65. 

Ij  1 j 

-16. 

ci  cr. 

17. 

C C* 

P A „ 3 8 

84. 

c 

j 

0 

0.  49 

0.  08 

7 ! j „ 

i j 0 

cr 

Q p 

1 p,  . 

9P 

86.  89 

84. 

cr- 

0 

0.  54 

0.  08 

“7  cr 

1 j Ij 

•”«  cr 

cr  ”? 

1 A n 

A A 

96.  99 

84. 

cr 

j 

Ij 

ij . 59 

0,  08 

9 0 . 

Ij  I"1 

-44, 

1 9 

1 p. . 

48 

86.  4 0 

84. 

cr 

.j 

Ij 

0 „ 6:8 

0,  08 

O c: 

1 j 1 j 

-5  0. 

59 

16. 

cr 

C J 

8p. . 4 ij 

84. 

cr 

j 

Ij 

0 „ 66 

0 . 03 

9 ij . 

ij  n 

cr  ~i 

A A 

I A B 

16 

86.  4 0 

84. 

cr 

Ij 

0 . P 8 

0.  08 

97  „ 

Ij  I j 

-54  „ 

fi  ij 

16. 

18 

86.  4 t 

84 . 

cr; 

Ij 

I I „ h;i  9 

0.  fi  9 

1 0 fi . 

Ij  Ij 

-49. 

H 7 

16. 

18 

86.  41 

84. 

cr. 

Ij 

0. 63 

0 . ij  9 

1 05, 

Ij  Ij 

-33. 

45 

16. 

18 

86 . 4 1 

84. 

cr; 

Ij 

0 „ p,  P 

0.  09 

110. 

Ij  Ij 

l“  7 ii 

17 

16. 

18 

86.41 

84. 

cr 

Ij 

ft . P'  4 

0.  09 

115. 

0 ij 

A 

t Cl  n 

IP. 

16. 

1 3 

86.  41 

c 4 n 

cr 

. t 

o 

0 „ p.  i j 

ij.  09 

1 9 ij . 

Ij  Ij 

-?  . 

08 

1 P< . 

16 

86.  41 

84. 

cr 

.j 

Ij 

fi . 5P 

i'i  „ 09 

i 1=  cr 

Ij  Ij 

fl . 

9 

16. 

81 

86. 41 

84. 

cr 

0 

0.  53 

0.  09 

1 3 ij .. 

Ij  Ij 

1 . 

14 

16. 

88 

86.41 

84. 

cr 

j 

o 

ij  „ A 9 

0 . 0 9 

135. 

Ij  1 j 

1 . 

P-5 

16. 

-.7  i 

86  „ 4 1 

84. 

cr 

fl 

0.  45 

0.  1 j 9 

1 4 ! j „ 

ij  fi 

1 . 

9 P 

1 6 „ 

48 

86.  4 1 

84. 

cr 

_i 

Ij 

0 . 4 1 

rC  09 

145. 

Ij  Ij 

c ■> 

05 

16. 

p.  ij 

86. 4 1 

84. 

cr; 

Ij 

lj„  97 

0 . 09 

1 O'  '■* . 

Ij  I'1 

cr.  * 

1 1 

1 6 . 

}'  7 

86  „ 4 1 

84. 

cr 

.j 

Ij 

0.  83 

09 

ENTER  t TO  CPLCUL.RTE  HTC‘7 1 

1 HP  H 1 1.  I Pi 1 948P  ■4  1 F + ! f i ~ ' 

t i = ? . o o ij  i'i  o ri  n f - o £■ 

T £:  = 1 1.  1 ij  i j n ij  i j o p - 1 j t 

P F‘  □ I •“  P h M ■“  J n P H T 1 9 0 1 j 


E 
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PAC  SAMPLE  RUN  - ABBREVIATED  OUTPUT 


RUN 


PRC 


16:£3EDT 


09.-'  fi  ftp 


INPUT  FILE  NAMESTRCST7P  • 

ENTER  1 IF  YOU  WANT  FULL  LIST  DF  OUTPUTS  ENTER  £ IF  YOU  WANT  RBEPEVIRTED  LIST.TS 


INPUT  VALUES  — INPUT  UNITS  < MSEC.  MF'H ■>  DEGREES.  INCHES.  LBS.  FT-LBS.  G'S> 
INITIAL  VELOCITY:  34.0 

INITIAL  HEAD  ANGLE:  -9.50 


INITIAL  TORSO 

ANGLE: 

-£7.  5 

MLEG 

M TOR SO 

MS TERN 

MHEAD 

RT 

RN 

PH 

PTOPH 

7 1 . 0 

58.  4 

£ . 5 0 

11.4 

14.0 

£ 0 . 5 

4.75 

C.O  ■ >' 

NPTS  NECK 

NPTS  KR 

NPTS  VEH 

NPTS  SEAT 

NPTS  GAS 

SL.  ST 

SL  . KR 

T; 

8 

18 

6 

14  0 . 

500E+04 

0. 840E+04 

NPTS  STEP 

NPTS  CHST 

5 

5 

GAS  FLOW  TIME 

0. 

15.0 

1 8 . 0 

££.  0 

. £5.  0 

30.  0 

35.  0 

40.  0 

45.  0 

5 0.  0 

GAS  FLOW  TIME 

55.  0 

60.  0 

75.5 

1 0 0 . 

GAS  FLOW  - LB-' 

SEC 

0. 

0. 

3 . 4 0 

3 . 9 0 

7.  54 

8.  09 

7.50 

5.87 

4 . 0 0 

£.54 

GAS  FLOW  - LB-'’ 

SEC 

1 . 49 

0 . 8?  0 

0. 

0. 

E - 13 


SEPT  FRICTION 

DI SPLPCEMENT 

-50.  0 

0. 

1 . M 

14.  0 

15.  0 

50.  0 

SEPT  FRICTION 

FORCE  - LBS 

p 

0. 

0. 

350. 

35  0. 

0. 

0. 

NECK  PNGLE 

-81 . 0 

18.  0 

90.  0 

NECK  TORQUE  - 

FT-LBS 

117. 

0. 

-87.  0 

VEH.  PULSE  - 

TIME 

0. 

7.  0 0 

17.0 

£ 1 . 0 

£6.  0 

33.  0 

45.  0 

5£.  0 

54.  0 

53.  0 

VEH.  PULSE  - 

TIME 

85.  0 

70.  0 

76.  0 

80.  0 

90.  0 

93.  0 

1 07. 

150. 

VEH.  PULSE  - 

DECELERATION 

0. 

0. 

13.  3 

1 £ . 7 

15.  6 

10.7 

18.  3 

c'9 . 8 

3 0 . 7 

£9.  8 

VEH.  PULSE  - 

DECELERATION 

£ 1 . 3 

C.  s'  ■ C 

17.  3 

19.  5 

11.7 

1 £ . 6 

0 . 

0. 

KNEE  DISPLP CEMENT 

0. 

£.  00 

£ . 5 0 

£.  75 

3.  0 0 

3 . 5 1 j 

4.50 

15.0 

KNEE  FORCE  - 

LBS 

0. 

0. 

£ 0 0 . 

7 0 0 . 0 . 

180E+04  0. 

£3  OF  + 04 

0. £50E+04 

0. £70E+04 

STERNUM  D I SPLPCEMENT 

-50.  0 

0. 

0.  £50 

1 . 0 0 

1 0 . 0 

STERNUM  FORCE 

- LBS 

0. 

0. 

4 0 0 . i J . 

160E+04  0. 

£5  0E+05 

CHEST  DTSPLPC 

EMENT 

- 1 1 . £ 

-1 . £5 

0. 

1 . £5 

1 1 . £ 

CHEST  FORCE  - 

LBS 

-0. 465E+04 

-150. 

0. 

15  0.  0 . 

465E+04 

PTMDP 

PGZ 

GTZ 

I_l 

PN1 

PN£ 

PN3 

14.  7 

0.  0.16 

6E  + 04 

66c . 

1.40 

1.40 

1.40 

VC1 

VC£ 

PV 

FSP 

FSB 

FSC 

ftl 

VI 

i j . 7 0 0 

0 . 7 0 0 

4.  00 

1 3 . 0 

£.  00 

11.0 

13.  4 

£4.5 

PO 

THETPD 

FPBMGT 

ST DP MR 

C DP  MR 

DMS 

DINF 

USDCK 

9.50 

0 . 

3 . 4 0 

0. 

3.  0 0 

6.  00 

4 . 0 0 

1 3 . 0 

I..JH 

DROLLZ 

C 

V £ Z 

I.JB 

LF 

DCN 

6 . 1 0 

5.  0 0 

31.7 

8 . 3 3 

18.4 

16.7 

0. 

THFO 

THLO 

ft  ST  OP 

STEP 

P I NT  1 

PINTS' 

19.  0 

44.  0 

0.  150  -0. 

100E-08  0. 

500E-08  0. 

5 0 OE - 0£ 

£ - 14 


TIME 

CST  R GS 

HD  E GS 

STM  ACC 

FEM  FORCE 

REV  WE  CST 

ERG  PRES 

CMS'1  ~ , 

05 S> 

05  -■  S> 

cg  S) 

CLES> 

ii 

ii 

LL  II 

Ll  II 
II 

CPS I G) 

0. 

0.  00 

0.  00 

0.  0 0 

0. 

0. 

0. 

5.  00 

0.  £4 

0. 4£ 

-5.  9 0 

0. 

0.  00 

0. 

1 0 . 0 0 

0 . o'  d 

0 . 38 

-7.51 

0. 

- 0 . £ 0 

0. 

1 5 . 0 0 

0 . 3 1 

0.41 

-7.  6 0 

0. 

-1 . 33 

-0.  03 

£0.  00 

0.  37 

0. 6£ 

-7.  67 

0. 

-45. 71 

11.99 

35.  00 

0.  75 

1 . 09 

-7.41 

0. 

-184. 94 

£.81 

30.  0 0 

1 . 4£ 

1 . 86 

-6.96 

0. 

-1 13.98 

5.  34 

35.  0 0 

1 . 73 

£.  19 

-6.54 

0. 

—9.  76 

1 . 73 

4 0.  00 

1 . 69 

£.  1 £ 

-6.51 

0. 

-11.95 

1.40 

45.  0 0 

4 . 83 

0.  79 

-100.41 

5£.  66 

-14. 49 

0 . 9 0 

50.  0 0 

1 0 . 8 0 

6.  43 

63.  07 

546. £5 

-4.64 

0.  49 

55.  00 

14.50 

13.31 

1.70 

1 039. ££ 

0. 

0.  9£ 

60.  00 

18.  58 

11.0  0 

—6.  76 

1 079. 40 

0. 

a • CL  »*’ 

65.  00 

£4.  16 

14.37 

~16.  c'c! 

1 096 . 3 1 

0. 

3 . 6 1 

70.  00 

30.  77 

£7.  19 

-£5. 9£ 

1 087. 19 

0. 

4. 94 

75.  00 

38.  14 

4 0.  84 

c*  -? 

1 03  0 . 9£ 

0. 

6.17 

80.  00 
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